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TTie combined use of optimi 2 ation and simulation tech- 
niques proves to be efficient for planning and designing of large 
water resource systems. Optimization technique selects the best 
out of many generated alternatives for a specified objective, 
while siimilatlon attempts to study the performance of a specific 
alternative. These techniques are applied, in this thesis, to the 
Narmada river basin in the central zone of India, The system 
consists of five serially linked reservoirs located on the main 
stem of river, A mathematical model is used to represent the 
main featxires of the system. Four periods in a year are selected 
for use in the discontinuous linear programming (LPD) model. 

ZX3LP (IMSL) computer package program available on DEC 1090 
computer system is used to run the LPD model. This model is 
implemented for different cases such as different configurations, 
different flow conditions, lirwsar and nonlinear objective 
functions, etc., to obtain 14 design variables namely, gross reser- 
voir capacities, seasmal water releases for irrigation and 



installed capacities of power houses. Deterministic inflows 
are used. LPD is used in the preliminary design stage. 

If the system analysis is based only on the LPD model, 
the design of system may be overoptimal or underoptimal. Hence 
the results of the LPD model are needed to be screened further. 

A range for each variable is assumed on the basis of the results 
of the LPD model for use in steepest ascent method of sampling. 

A set of design variables from this range is selected arbit- 
rarily as an initial base. The jump or step value is set to 
about one-fifth of its range for each variable. The step values 
are decreased in successive iterations. A nonlinear objective 
function developed for LPD model is used to represent the 

response surface, A net maximum annual benefit of Rs. 120,65 x 
7 

10 is obtained. The set of variables which gives maximum net 
benefits is selected for simulation study. Values of the 
selected variables namely gross storage in MAP, annual irriga- 
tion target in MAP and installed cajfecity in MW for each reser- 
voir are 6.25, 4.00, 9.00 (Bargi), 14,82, 10,30, 240,00 (Narmada- 
sagar), 1.60, 4,14, 108,00 ( Omkareshwar ) , 0.41, 0, 67.00 
(Maheshwar), 7.76, 9.52, 125,00 (Sardar Sarovar) respectively. 
Maheshwar reservoir is designed only for hydropower and hence 
irrlgaticn component is zero. 

A computer porogram developed by HEC (1966) namely 
'Reservoir Yield' is modified sioitafoly, to fit the present 
system for simulation. Simulation attempts to study the perfor- 
mance of a system for specified operating policies. The opera- 
ting policies inclucfe: (i) water release policy (li) the maximum 



and minimiOT permissible storage® for each period and for each 
reservoir. Simulation study requires long sequences of stream- 
flews. Since only 28 years of streamflow record is available 
at. gauging sites at Jamtara, Mortakka and Garudeshwar^ synthe- 
tic flows are generated using stochastic hydrology preserving 
important statistical parameters of the historical streamflows. 
HEC-4 computer program: ‘Maithly streamflow simulation* is used 
to generate 250 years of synthetic flows at each reservoir site, 
using transferred streamflow data available at each reservoir 
site , 

Two types of simulation studies are made: (i) moithly 
simulation and (ii) flood simulation. Four monthly simulation 
runs are taken for NWOT award with and without modification, 
using historical and synthetic flows (NWDT: Narmada Water 
Disputes Tribunal). Rigid rules are developed for each reser- 
voir using end of month storages. Shortage index and types of 
shorta^s for irrigation and power are evaluated at each reser- 
voir site and for each simulation run. Loss functions are 
necessary to take into account the contingency of not being able 
to meet irrigation and power targets. Hcaithly statistics of 
spills from the terminal reservoir nan«ly Sardar Sarovar are 
evaluated for each simulatiai run. Tine net average annual 
benefits from the system are Rs. 66.37 x lo'^ and Rs. 73.93 x 10 
for modified NWiyp. and NWDT award respectively for monthly simu- 
lation, In the case of flood simulation the net annual benefits 
are Rs, 71.53 x 10*^, and Rs. 80.13 x 10*^ for modified NWOT and 


NWDT award respectively 



Flood simulation study is carried out to find the 
maximum 3-hotjrly spill from the terminal reservoir namely 
Sardar Sarovar. The Lag method of flood routing is used for 
channel routing, A flood magnitude of 24,50,000 cusecs gets 
attenuated to 5,42,407 cusecs for NSnTOT award, while it gets 
attenuated to 11,00,840 for modified NWDT award. 



CHAPTER 1 


INTRODUCTION 


1.1. General 

Water Is an essential basic resource which is necessary 
for the existence of all life form. It is required in abun- 
dance to meet irrigation, domestic and industrial, etc., needs. 
Its availability is restricted in terms of quantity, quality, 
time and space or ccanbination of these factors. 

The water resource systems are complex hence assumptions 
are made to simplify it. The complexity of the problem may be 
due to the complicated inter-relationship between various compo- 
nents of the system, data limitations, uncertain hydrology, 
demand pattern, etc. 

The design of a water resource system consists of 
(A) determination of physical size to yield certain target out- 
puts and (ii) establishment of an optimal operating policy. 

This involves the interaction of political, legal, governmental, 
economic and engineering processes. In the development of any 
river basin, generally, the first priority is given to irriga- 
tion. The water requirements for power generation, domestic 
and industrial purpose are given second choice. The reason is 
that there is no other substitute available for the irrigation 
water, while power may be generated by the use of coal, oil, 
nuclear energy. It is observed that large habitations developed 
first in the delta reaches oE the river basin because of easy 
availability of water in large quantities. 



System analysis is used to search out 'best* solution 
in terms of objective functions by joining engineering and 
eccnomic aspects together- The data required for analyzing 
water resource problems include hydrologic and hydrometeorolog- 
ical records, cost and benefit functions, economic loss functions 
and demand patterns. In general benefits cannot be estimated 
with absolute reliability. 

The system components and their interrelationships are 
first identified. An idealized representaticn of the system is 
made through a mathematical model by making simplifications and 
assixmptions , Care must be taken to see that the model is still 
a valid representation of the problem during the process of 
simplifying the system. An objective ftmction is a very 
important part of mathematical model. It is described by a 
mathematical expressicn. It quantifies the effects of alter- 
native courses of action in the system model. 

All possible short -comings, errors or oversights in 
the model should be Identified in advance before its results are 
implemented, t>ata requirements should be scrutinized for 
compatibility with data availability. Solution for the model 
will ally be an approximation to the solution of the real world 
problem. 

1.2, Optimization Techniques Commonly Used 

It is essential that river basin model and selected 
optlmizaticai technique are compatible. Tte linear prograiwnlng 
(liP) and dynamic larogramming COP) or a joint use of these two 
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are the commaily used optimization techniques for solving the 
water resource problems. Their uses are restricted to preli- 
minary screening. DP has limited use if there are many state 
variables in each stage. Its use ie limited due to the so 
called "curse of dimensionality" (Bellman, et al,, 1962). 

Simulation model is ixsed to evaluate the performance 
of various cc»if iguration of reservoirs, px^wer plant cap>acity, 
water use allocation and operating policies, for given hydro- 
logical inputs. 

Optimization model is designed to select the best plan 
out of many alternatives for a specific objective. Simulation 
model is better suited for the detailed analysis of specific 
alternative. The final design and operating policies is 
selected csn the basis of results of these studies. 

1.3. Present Study 

The aim of the present research study is to determine 
optimal operation of five serially linked reservoirs located 
on the main stem of the river Narmada using optimization- 
simulation technique. Discontinuous linear programming model 
(LPD) with deterministic inflows is \ised to arrive at design 
variables such as reservoir size, installed capacity of power 
house and seascmal water releases for irrigation and pxDwer 
for each reservoir. Further screening is done by using the 
steepest ascent method. A set of design variables which 
produces roaxlmum net benefit is selected for the use of simu- 
lation study. A simulation computer program developed by HEC; 



"Reservoir yield" is modified and lased to test the efficiency 
and performance of the system for specified operating policies. 

Simulation study requires a long record of streamflows. 
Historical record of streamflows are too short to irclude all 
possible pattern of drought and floods. Generation of synthetic 
streamflows provides longer sequences, Streamflow data are 
generated using technique of stochastic hydrology. 

A monthly streamflcw simulation computer program has 
been developed by Hydrologic Engineering Centre (HEC-4), The 
monthly streamflw data are generated for a desired length 
using this program and with available record of historical 
streamflows as input data. A number of simulation studies are 
carried out using historical streamflcws and synthetic stream- 
flows as input data. The other input data are: design para- 
meters# release policies# and maximum and minimum permissible 
storage in each period and for each reservoir. 

1,4. Organisation of Jteport 

Chapter 2 consists of various methods and techniqi^s 
available and their ULmitations to analyze water resource 
systems. Chapter 3 deals with various techniqi^Kss of generation 
of synthetic streamflows. Chapter 4 consists of the use of LPD 
model for prelirainary screening. Chapter 5 includes the use of 
the steepest ascent method for further screening of design 
variables. The descriptim of simulation program is given in 
Chapter S, The jcesults obtained from LRD model# steepest ascent 
method and simulation studies are described in Chapter 6. In 
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the same chapter discussion of results are included. Surwnary 
and recommendations for further studies are given in Chapter 7. 

1,5. Units Used in the Present Study 

In order to keep the study as close to a real problem 
as possible, actual data from published literattore, and from 
the report of the Narmada Water Disputes Tribunal (NWDT), 
Volumes I to IV are used. These data are generally in PPS 
system and hence the sane units are retained in this study. 



CHAPTER 2 


REVIEW OF LITERATURE 


2.1. Introduction 

Mathematical programming and simulation techniques 
have been applied to water resource systems, during last few 
decades. Research has progressed frcxn analysis of a single 
reservoir to a system of multireservoirs, encompassing an 
entire river basin. Harvard Water Group did pioneering work 
to study the basic principles, concepts and methodology for 
designing and planning of water resource systems, 

A combination of linear programming, simulation and 
response surface methods was used by the Texas Water Develop- 
ment Board CTWDB), Their aim was to seek *near optimum* solu- 
tion rather than exact optima* The procedures adopted by them 
are as undter: 

1 . The components of the system and operating rules are 
determined by an optimization method. 

2. Initial screening is carried out by using the method of 
random sampling. 

3. A gradient search is used to reduce further range of 
variables , 

4. The most attractive schedules are improved by succe- 
ssive perturbations. 

5- A pattern search (Hooks and Jeeves, 1961) is used to 
reach the optimum in the vicinity of the crevices 
(Hlramelblau, 1974), 



2.2. Available Analytical Techniques for Sttjdy of Water 
Resource Systems 

2,2,1, Linear pr ogranroing 

Linear programming is widely used in water resource 
systems. Manne (1962) has shown that the problem of infinite 
planning horizon can be solved by using linear programming with 
the assumption that the objecti^re function and constraints are 
a linear function of decision variables. Thonas and Watermeyer 
(1962) used the linear programming forroulaticai of Manne for 
seqtjentlal decision model, Tlieir main objecti-^^ was to deter- 
mine the optimal draft rate for each period and optimal capacity 
of reservoir to maximize the benefits. Their decislcn regar- 
ding releases is based on the initial storage and random inflows 
during that period. 

There arc two types of linear programming models namely 
linear programming continuous (UPC) and linear prograitaning 
discontinuous (LPD). LPC obeys the continuity law for one or 
itKare periods while LPD violates the contlnxiity law for one or 
more periods. Loxicks (1981) advocates the xise of both models, 
Dorfman (1962) illustrated the use of mathematical programming 
techniques^ especially# linear programming and a special type of 
nonlinear programming# namely separable programmlr^, for solving 
water resource systeire, Martinez (1971) used deterministic 
LPD for solving a single multipurpose reservoir. Rogers (1969) 
used a linear programming formulation for the lower Granges and 
Brahmaputra rivers in India and Bangladesh, He considered 
seven resexrvoir sites in India 'and a single site in Bangladesh. 
The objective function includes benefits due to power and 



irrigation^ cost of surface and ground water storages^ and 
cost of eniban>aTnents, ail reduced to annual costs. 

2.2.2. Chance constrained programcning 

One of the techniques used to ta3ce into accoiant the 
stochastic nature of inflow is the chance constrained program- 
ming technique. This method was developed by Chames aol 
Cooper {1958). 

The linear decision rule (LDR) was first developed by 
Re-^mlle, et al. {1969). He applied this technicfue for finding 
out optimal capacity and operating rules for a single reseirvoir 
using stochastic inflows, Rewlle and Kirby (1970)^ Joeres, 
et al. (1971), Nayak, et al. (1971, 1974), Estman, et al. (1973), 
rnClerc, et al, {1973) and Revelle and Gunc3telach (1975) have modi- 
fied, extended, and/or applied this method to reservoir manage- 
ment problems. Eisel (1972), Loxieks (1970), Sobel (1974), have 
pointed out tlie drawbacks of lirraar decision rules. Loueks 
and Dorfman (1975) recccwnended that these rules may be used for 
initial screening study, as they are conseirvatlve for design 
purpfjse. Lane (1973), Askew (1974) have vaed this method in 
their study. Houeks and Datta (1981) developed multiple LDR 
model. It is superior to single IDR model for equivalent rest- 
rictions on reservoir operation and performance , The multiple 
LDR model gives a smaller capacity reservoir than a single LEtfl 
model, 

Dantzlng and Wolfe decomposition principle can be 
used to reduce riKJdeling effort as well as to save comptiter time, 
when the system is very large. In case of nonlinear objective 



ftincticm and a few ncailinear constraints/ a piecewise linear- 
ization technique is fovind suitable. The problem beccmes 
complicated and unmanageable as the number of nonlinear const- 
raints increases. 

2.2.3. Nonlinear programming 

Nonlinear programming techniques are not ccmmonly tised 
for the analysis of water resource systems. Few authors who 
have used this technique are Chu and Yeh (1978)/ Hiroroelblau 
(1977)/ and Bayer (1974), Bayer used nonlinear programming in 
river basin water quality model. He compared his results with 
the results obtained by other researchers using linear and 
dynamic programming to solve river basin water quality optimi- 
zation problems, Himmelblau used two nonlinear progranming 
techniques; (i) a Generalised Reduced Gradient (ii) a Conjugate 
Gradient Projection method in his water quality model, 

2.2.4, Dynamic programming 

Ttere is vast literature available for DP technique 
used by many researchers. This technique was used by them in 
their study of a single multipurpose reservoir. The reasons 
were (i) dimensionality of the problem was well within the 
realm of the computer memory and (ii) nonlinear objective 
function. Little (1955) used DP for optimizing the operatlcn 
of a hydropower system. Buras (1963) used DP for solving 
conjuctive operation of a reservoir-grourKi water system. Hall 
and his associates (1961/ 1964/ 1966, 1968) proposed a number 
of DP formulaticans to analyze a multipurpose reservoir system. 



Their aim was to find optimal schedule of releases such that 
the maximum benefits are obtained from the sale of water for 
irrigatican and for power. Young (196?) ccarblned deterministic 
DP and hydrologic slmulaticxi to determine the optimal draft 
from a single reservoir. Mobasheri, et al. (1970) xased modi- 
fied DP f or determining an optimal long term operating policy 
for a single reservoir. Parlkh (1966) used DP for a single 
reservoir with a deterministic hydrology and used a deccxnposi- 
tion principle of linear programming. Meir and Brighter (1967) 
used DP in their study to optimize a multistage water resource 
systems. Schwing and Cole (1968) used DP to stxKJy the optimal 
control of linked reservoirs. Lieu and Tedrcw (1973) used DP 
in their study of a multi lake river system to find operating 
rules, Harley and Chidley (1978) used deterministic DP to 
determine long term reservoir operating policies for a reser- 
voir system. Collins (1977), and Oprlcoblc, et al, (1976) used 
DP in ttelr study to find optimal long teinn control of a imlti- 
purpose reservoir, Bhaskar, et al, (1980) used DP for deri- 
ving mcnthly reservoir release policy. 

2.2.5. Stochastic dynamic programming 

The stochastic nature of inflows are taken into accomt 
in dynamic programming by using stochastic dynamic programming. 
Pew authors who have used this technique are Butcher (1971), 
Dudley and Burl: (1973), Su and Deininger (1974), Mawar and Thorn 
(1974) and Torahi and Mobasheri (1973). 

Torahi aid Mobasheri used stochastic dynamic program- 
ming for detejcmining the optimal operating policy of a single 



multipurpose reservoir. They developed a dynamic programming 
model with a physical equation and a stochastic recur sicai 
equation for deriving the optimum operating policy. Butcher 
used this technique for a multipurpose single reservoir. In 
his study, the optimal operating policy is stated in terms of 
the state of the reservoir (storage volume) and riverflcw JLn 
the preceding month. 

2,2.6. Incremental dynamic progrananing with successive 

appr oximat i ons 

This method has been used by Troat and Yeh (1973) and 
Nupmongcol and Askew (1976), Troat and Yeh used this roethiodi 
to determine the optimal design of a system of reservoirs with 
series and parallel connections. The retxjm from the system is 
determined for an optimal operating policy and a specified set 
of reservoir sizes. This policy is detezrmined by decomposingi 
the original problem by Bellman's method of sxiccessive appreaci— 
mations. The multistage problem is decomposed into a series 
of stjbprcblems of one state variable, A sequence of optimiza- 
tion over the subproblems converges to the solution of the 
original problem. The costs are assumed to be a function of 
reservoir size and are computed from storage capacity versus 
Cost curves, A modified gradient technique is used to deter- 
mine the set of reservoir sizes which maoclraizes the net benefits 
with the imposed constraints, Nopmangcol and Askew developed 
a multilevel Incremental dynamic programming. The basic caneepat 
is directly pendent on certain characteristics of incremental, 
dynandc programming (IDP). They found marked reductlca in 



computing time by using this technique. It increases the power 
of IDP. It can be used as a tool in the optimization of multi- 
dimensional deterministic systems, 

2.2.7, Incremental dynamic programming 

State increment dynamic programming technique was 
introduced by Larson (1968), A different version of an incre- 
ment concept for state variable has been xxsed by Hall^ et al. 
(1969) and is known as Incremental dynamic programming. ‘Ihe 
major difference between them is the^time Interval used in 
computation. Time interval Is a variable in the former while 
it is fixed in the latter. 

2.2.8, Multiobjective dynamic programming (MC30P) 

This has been \ised by Tauxi# et al. (1981) 

They developed mathematical models for determining release 
policies that maximize dump energy and minimize cumiilatlve 
evaporation . The major characteristics and advantages of MOOP 
are as trrxier: 

(a) NcjncommensurablQ objective functions can be handled 
quantitatively, 

(b) The entire ncnlnferlor solution set is available from 
one solution of the problem, 

(c) Trade-off ratio are available explicitly between objectives, 

2. 2 .9, Markov models 

A discrete probability distribution concept is used in 
Markov model* The following two apprcociroatlons are made to 



describe the probability distributicai of streamflows, reser- 
voir volumes^ and other hydrologic events. 

(a) The discrete probabilities can be approximated by 
continuous probabilities. 

(b) Discrete Interval of time can be approxiroated by 
continuous time interval. 

Many other approximations are also necessary in addi-tion 
to the above, to preserve mathematical tractability and essence 
of an actual water resource system (Shen, 1976). Gablinger, et 
al, (1970) used Morkov model for flow regulation. 

2.2.10, Monte Carlo technique 

Monte Carlo analysis is valid for sequential decision 
processes. It is used to test the policy, using equally likely 
sequences (Hall and Daracup, 1975), Askew, et al, (1970) used 
Mcnte Carlo technique in the design and operation of multi- 
purpose reservoir system. 

2-2-11, Queuing theory 

This can be applied to the problem of selecting the 
optimal design of a single multipurpose reservoir. Langbeln 
(1958) used this technique to find reservoir capacity, by main- 
taining a minimum regulated draft, providol that inflow is 
normally distributed aixi draft is a linear function of storage. 
Gani and Moran (1955) have given the outline for using queuing 
theory and Monte Carlo analysis, Thowias, et al, (1982) gave an 
elagant (queuing theory solution for selecting the optimal scale 
of development and optimal oierating policy for a single 



multipurpose reservoir. Fierlng (1961) applied this technique 
associated with Monte Carlo method to study the optimal size 
of a single multipurpose ireservoir using an economic efficiency 
objective . 

Mattematical progranroing technique has been extensively 
used by many iresearchers for planning and designing the water 
resoxjrce system and for finding operating policies. There is 
no general algorithm. Ihe choice of the technique depends 
mostly on the characteristic of reservoir system^ availability 
of data, objective and constraints. 

2.3, Limitations of Optimization Models 

The optimization model is \inable to give exact solution 
to the problem dx® to following llndtaticns, 

CD Many assumptions are made to simplify the rmal problem 
in order to make the problem solvable. The objective 
function (cost and benefit function), energy production 
function, evaporation function, flood damage function, 
loss functioi, etc*, axe generally nonlinear. One has 
to take recourse to nonlinear programming technique in 
such situations which takes lot of cat^uter time and 
thus proves to be costly . Nonlinear programming 

problems have no general algorithms fur solution. To 
overcome this, linearization may have to be done. 

(2) The second llmltaticai is mcace ccacepfcual. It stems 
from Inability to fix a criterion for evaluating each 
possible management alternative with regard to quanti- 
fications, : 


<3) The third is nonavailability of reliable and cc»tq?atible 
data. Hence the optimization models are restricted to 
preliminary screening designs (Biswas, 1976). 

Stochastic optimization model contains more variables 
and constraints than that of deterministic model of the same 
system. It requires huge computer memory, hence it is rest- 
ricted to small subbasins. Houok, et al. (1978) used stochastic 
linear programming irodel for design and management of multi- 
purpose reservoir system. 

2.4. Simulation Models 

Slmulaticn model is capable of evaluating the hydrolo- 
gical and econcmic performance of a large scale river basin. 
Ackoff (1961) defines simulatlcn as a process which "duplicates 
the essence of the system or activity without actually attaining 
reality itself”. Naylor (1973) defines simulation as a “numer- 
ical technique for conducting experiments with certain types of 
mathematical models which describe the behaviour of a complex 
system on a digital computer over extended periods of time”. 

He states fxjrthor that "the principle difference between simu- 
lation experiments and real world experiments is that with 
simulation the experiments are conducted with a model of real 
system rather than with the actual system itself”. 

Six major phases of a simulation have been listed by 
Chtirchman, Ackoff and Amoff (1957): (i) formulatjlsn of problem, 
(ii) constructing mathematical model to represent the system 
under study (iii) deriving solution from the study (iv) testing 
the model and solutlcn derived from It (v) establishing control 


over the solution and (vi) putting the solution to work ( Imple- 
mentation). Simulation model takes into account nonlineaur, 
dynamic, and stochastic responses of the system. 

2-5. Optimlzaticai-Simulation Technique 

Considering the complexity and Inherent characteristics 
of a reservoir system, a joint use of optimization and simula- 
tion technique is found more stiitable. Seeker and Yeh (1974) 
used this technique for their study of operation of multi- 
reservoir system, Jacctoy and Loueks (1972) used analytical and. 
optimization models to screen the set of design variable worthy 
for slnnalation analysis. The combination of LP and stochastic 
DP was developed by Takeuchi and Moreau (1974). Roefs and 
Bodin (1970) classified optimization into two approaches, the 
explicit stochastic optimization (ESO) and implicit stochastic 
optimization (ISO) . The ESO uses the probability distribution 
of inf lews instead of inflow sequences. In the ISO, the 
generated inflow sequences and deterministic optimization 
method are used to determine optimal releases, Hirsch, et al. 
(1977) Studied the gains from the joint operation of mixLti- 
reservoir system, Hufschmidt and Fiering (1966) carried out 
systematic study of a system of reservoirs by simulation. 
Fredrich and Beard (1972) used simulation to analyze the water 
resources problem Including (i) determination of the specified 
set of projects to be constructed, (li) timing of construction, 
and (ill) method of operation, so as to minimize the sum of 
capital and operaticn cost over the planning period. 



Reservoir zoning concept was given by Beard (1967), 
According to this concept^ each reservoir is divided into a 
niirriber of storage zones and during sirnulatiun all reservoirs 
are maintained in the same zone as far as possible, operating 
policies are based on priority concept or some prescribed inter- 
reservoir relaticnship. A generalised simulation program for 
tl» operation of reservoir system for conservation purposes 
such as water supply# recreation# lew flow augmentation# and 
hydroelectric power was developed by Hydrologic Engineering 
Centre# U.S. Army of Corps of Engineers (HEC-3# 1974). 

Ford and Fulkerson (1962)# Durbin and Kroenke (1967) 
developed a new approach to simulate a complex water resource 
system using out -of -kilter algorithm (OKA) to a network of 
reservoirs, A system can be represented by a series of nodes 
and arcs in a 'capacitated network' form analogous to an 
electric circtiit. Texas Water Development Board (TWDb) deve- 
loped a nurrdDer of computer programs (TWDB-1# 1970? TWDB-2# 19727 
TWDB-3# 1972) using this technique, Sigvaldascn (1976) used 
the concept of reservoir zoning and OKA to simulate tte opera- 
tion of a multipurpjose and multireservoir system. 


2,6. Streamflcw Generaticn Techniques 

Simulation study requires a long seqence of streamf lows. 
In most of the river basins of India# the historical records 
are not long enough to make predictlcais aboiit mean# high and 
low flows. Hence it is not pxjssible to assess reliable perfor- 
mance of the water resource system. A Ic^ig sequence of flow data 



is generated using technique of stochastic hydrology, preser- 
ving the important statistical parameters of historical flew 
data. 

After testing the system with possible alternative 
sequence of hydrologic input data, the planner may be able to 
formulate effectively and efficiently the development proposals. 
"Synthetic flews do not Improve poor records but merely improve 
the quality of design made with whatever records are available" 
{Piering and Jackson, 1971), A design based only on short 
historical sequences may lead to under or over design to an 
xonknewn extent. 

The data generating models can be classified into two 
groups- short memory models, and long memory models, Markov, 
and ARIMA (Auto-regressive-integrated moving average) models 
belong to the short memory group, while Fractional Gaussian 
Noise (PGN), and Broken line (BL) processes belong to th^ long 
memory group. The term memory is tjsed as a synonym for persis- 
tence. It means a tendency for high flows to follow high fltws 
and low flews to follow low flows at a station. The short 
memory models have been criticized for their inability to 
simulate a relatively high valixj of Hurst ccjefficient (around 
0.7) observed in long historical record* long neroory models 
have been developed to take care of Hurst coefficient (Singh,! 979)^, 
Fierlng (1967) and Fiering and Jackson (1971) suromtirised 
the use of short memory processes. Mandelbrcut and Wallis (1969) .and 
Matalas and Wallies (1971) have considered long term hydrologic 
persistence in generating synthetic flow sequences using FGN. 



FGN model is fotjind to be expensive, Lettenmalr Buries 
(1977) have proposed a ndxed ARIMA-Markov model. 


2.7. Time Series Components of Streamflowi 

Four components are associated with stream £lc3w 

(Kottegoda, 1970): trend periodic correlation and 

random e^. if they are combined together by an additive model 

then the expression becomes = t + SE + + %• ^ 

t t t t ^ 

series is a sequence of values arranged in ordsr of their ^occur** 

rence. It is staticnary if the statistical properties are time 
invariant. The monthly data can be standardized by the following 
equation ; 



where, 

* monthly flows; 

X^ = standardized monthly flows; 

Xj - mean flow for the jth month. 

In general^ it is found that annual flows obey normal 
distribution^ while the monthly flows obey Log-normal distribu- 
tion. Markov or an auto~regresslve (AR) model of order p can 
be written as. 
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where, 

* AR process at time t; 


( 2 . 2 ) 



£ = AR parameters; t , '2- , •- 
e = random component, b - h'^i 


_ |!> 
-■\ 


Tte Lag-me Markov model has been widely used in gener- 
ation of annual flows, using the equation 


X. 


X^ 


.1 <1 -i2'i> 


1/2 g 


(2.2a) 


Persistence is one of the characteristics of time series. It 
can be calculated in term of serial correlation. Lag-one 
serial correlaticn is serial correlation between an event and 
the immediately preceding event, Kottegoda (1970) has reported 
that Lag -one models are quite reliable. 


2,8. Synthetic Streamflow Generation Techniques 
2,8.1, Short memory models 

The Russian mathematician Markov introduced the concept 
of Markov process. The probability distribution of any trial 
depends directly on the preceding trial, ihe following Markov 
model is proposed by Brittan to represent actual stream flows; 

X + fix. - X) t,S (1 - 

X as I I p. I -r-^r r-.i-r 4 . (2 3 

1+1 deterministic randan * 

where, 

^i+l' ^i * anhual rtmoff for (1+1 )th and 1th years; 

X - mean annual historical flows; 

y s= annual Lag-one serial correlation coefficients? 

tj^ = normal random variate of sero mean and unit 


variance (McMohan and Mein# 1978) 



Thomas and Flerlng (1962) developed the following algorithm; 


X 


1+1 


where^ 


X.^, + b.{X. ~ X,) + t, S,^, (1 - t2j 1/2 
j+1 J 1 j 1 j+1 


(2.4) 


^1+1' ^1 ~ generated flows dxaring the (1+1 )th, ith seasons; 
^j+1' ~ flews dxaring (j+l)th/ jth seasons; 

bj = least square regression coefficient for esti- 

mating (j+l)th flows from the jth flows and 

b = T fi±i . 

t^ = normal random variate with zero mean and xonit 

variance ; 

Sj+i^r Sj = standard deviation of flews dxiring the (j+l)th, 
jth seasons. 


Thomas and Burden (1963) transformed the normal variate t^ to 
a skewed variate t^ with an approximate Gamma distribution,, losing 
Wilson and HUferby transformation to take care of skewness 
in data; 




j-1 
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36 
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(2.5) 


where, 

V , s co-efficient of skewness of "like Ganroa" variate; 

J 

t^ « N(0, 1); 

t„ ®sG(0,l,t j)? 

r r,j 

j » repetetlvB annual cycle of season. 

McMohan and Miller (1971) reported that the Wilson and Hllferty 
transformation failed in case of large skewness, Matalas (1967) 



presented monent trans format ioi equations. Hie main assunption 

£3 

made^that logaritlms of flews follow normal distribution. 

The algorithms for estimation of parameters for three parameter 
lugncarmal model are as follcws; 


^i+1 ®j^^i 


Xj) - r 2 , 1/2 


( 2 . 6 ) 


Aj + exp{0.5 


(2.7) 


exp [2CSj + Xj)] - exp(Sj + 2Xj) 


( 2 . 8 ) 


exp [3Sj ] ~ 3 exp(S^) -f 2 
[exp(Sj ) - 1]^^^ 


(2.9) 


oxp [S^ R^] - 1 

- ■"» I I 

.As. A// \ __ 


[exp(Sj) - 1] ' [exp(Sj^j^) - 1] 


( 2 . 10 ) 


B. 




( 2 . 11 ) 


Moreau and Pratt (1970) generated weekly flews with the 
assxmptiun that weekly flows obey lognormal distribution, A 
moving average process is used to reduce sampling errors^ i.e.^ 
Smoot hening of raw data. They developed a nonstationary model 
for Synthesizing weekly flows such that means and variances 
of weekly and monthly flew will be preserved. A scheme for 
geTOrating daily flows at a single station has been developed 
by Quimpo (1968). Yevjevich (1966) used the concepts of stati- 
onary time series and Fourier analysis to produce sequences at 
a single site. Young and Pisano (1968) developed an operational 



model based on a ccxnblnation of the model suggested by Matalas 
and transformation studied by Yevjevich. 

2.8,2. Long memory models 

(i) Fractional Gaussian Noise (FGN): Mondelbrot and Wallis 

(1969) introduced the Fractional Gaussian Noise (FGN) model to 
take care of large val\^ of Hurst's coefficient, say h >0.5. 
Approximations are necessary to compute FGN variate which 
involves an infinite number of operations. They have used 
type 1 and type 2 approximations. Type 1 approximation is 
expensive with respect to computer time. Type 2 approximation 
is deficient in high frequencies. Low frequency approximation 
is satisfactory when h is close to one. Matalas and Wallis 
(1971) proposed a filtering technique to account for high 
frequency behaviour of type 2 approximation. The drawback is 
that it takes large amount of computer time. Mandelbrot (1972) 
developed Past Fractional Gaussian Noise (FFGN) to overcome 
the difficulty met with FGN. 

(ii) Broken Line models (BL) : Mejia, et al, (1972) presented 
a long term memory model, namely. Broken Line (BL) model. It 
is a continuous time process which produces Hurst effect. Lag- 
one correlation of flows, and the second derivative of correlation 
functim at the origin- The second derivative is related to 
crossing properties of streamflc:w function. Lawerence and 
Kottegoda (1972) have <3tefined the second derivative as the 
number of times the flow rate crosses a given level and the 



expected time between such crossings. This model has several 
defects. It is slightly non-Gaussian arid is unable to 
prodxice actual crossing properties. The definition of crossing 
properties is nut well defined. 



CHAPTER 3 


STREAMPLCW GENERATION MODEL 


3.1. Introduction 

The present study deals with studies relating to the 
optimal developnrasnt of the Narmada river basin. This basin 
will have several reservoirs for meeting demands of irrigation, 
power, etc. The historical record of streamflows at a few 
staticffis on the river is not very long. Hence there is a 
need to generate streamflows by using streamflow generation 
techniques . 

It is a simple procedure to generate synthetic flow 
sequences at a single site. It becomes cumbersat^e when more 
sites are included. It is desirable that the data generating 
technique must preserve important statistical parairoters of 
the historical flews at stations and at the sane time be able to 
generate canpatible flow data. 

It is a major problem, with the generating models used 
for multisites, to preserve cross correlation among various 
sites. A number of models are available for generation of 
monthly streamflows for multisite situaticxi. Looking at the 
llmltatiai of comHiter time, it has been decided to use 
'HEC-4: monthly streamflow simulation model’ developed by 
Hydrologic Engineering Centre, U.S, Army Corps of Engineers, 
for the Narmada river basin under study. It has been found 
that in general, annual flows obey normal distribution, while 
monthly f lows obey lognormal distribution . 



zo 


3.2. ■HEC-4 Model 

3.2.1. Mode 1 opt ions 

The HEC-4 model has the following opticais (Reference 

HEC-4 Manual). 

(a) Analyze monthly streamflow at a number of interrelated 
stations (up to 10 stations can be analyzed simultane- 
ously) to determine their statistical characteristics. 

(b) Pill-in missing flows and test for consistency of 
matrices of correlation coefficients. 

(c) Generate a sequence of hypothetical streamflows of any 
desired length, having statistics obtained in item (a). 

(d) Obtain maximum, minimum and averages of recorded, 
reconstituted and generated flows for each month and 
for specified duoration. 

The model can be used for other purposes such as 

rainfall, evaporation and demands, 

3.2.2, Computational steps 

(a) Increment the monthly flows at each station by 1 % of 

average flew of the corresponding month. The incremen- 
ting of historical data is to avoid possibility of 
negative logarithms while transforming the flew data 
into logarithmic form in the initial stage of compu- 
taticMi , The increment is substracted f ran the 
generated sequences in the final stage of c^omputatiens . 



(b) The mean, standard deviation and skew coefficient for 


each station and for each month aore computed using the 
follcwing equaticnsi 



log j + qj^) 


( 3 . 1 ) 
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Where# 

X - logarithm of incremented monthly flew; 

Q = monthly recorded flow; 

q « small increment of flow used to prevent infinite 
logarithm for month of zero flew; 

X = ntean logarithms of incremented mcnthly flews; 

N a total years of record; 

S « unbiased estimate of population standard deviate; 
g = unbiased estimate of population skew coefficient; 
1 * subscript used for month; 

j = subscript used for year. 


(c) Using the above derived statistical parameters# histo- 
rical flow data are first transforaed to Pearsen type III 
standard <teviate and then to normal standard deviate. 



Transformation to Pearson type III deviate 't ' is 
performed by; 
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(3.S) 


(d) Transformation of Efearscn type III deviate 't' to normal 
standard deviate ‘K' is performed by; 


(e) After transforming the flows for all mcxiths in respect 
of all stations to normal deviates, the simple correla- 
tlcm coefficient 'R' between all pairs of stations^ for 
each current and preceding month are computed by using 
the formula; 


’l - 1 1 


N 
< S 
j=l 


X 


i,j ""i-.l,j 


\2 , 
) / 


N 
( S 
j=l 



N 

E 

j-1 


X 


2 

i-l, J 


)} 


CN - l)/(N-2)] 


{3.7) 


in which X * X - X. 

Cf) If the sets of historical record are not completed for 

all stations, first phase of computation involves comple- 
tion of all sets maintaining consistency of all corre- 
laticxi matrices. When the sets of flows are complete, 
generation of hypothetical stream flow is acccsnplished. 

The serial and cross-correlatican coefficients 
(correlation matrices) form the basis of relationship between 
dependent variable (relating to the missing month or month to 



be constituted) and the independent variables {for stations 
having flows for that month already generated and other stations 
with flows of the antecedent month). The regression equation 
from the above relations are then computed for each station 
and month by Crout method. 

(g) The regression equation is of the form; 


tt * 

^i.n 


= ^,K.' , + 8„K! - + B^K. ^ + ... + B , k; , 

3jl/ 3 rii ***x ****X 


^n ^i-l,n ^n+1 ^i-l^n+1 


1 m 2. „(1 - rJ 


( 3 . 8 ) 


where ^ 

K* 5= monthly flow logarithm expressed as a normal 
standard deviate; 

^ = beta coefficient computed from correlation 

matrix; 

i = subscript used for mcnth; 

n = subscript used for station; 

m = niirober of interrelated stations; 

R multiple correlation coefficient; 

Z K random ntmtoer from standard normal populatlcm. 

(h) By the above formula the normal standard deviates are 
computed for each station in t\arn for one racnth at a 
time. Ihe process is started with average valties (zero 
deviation) for all stations in the first month. 


(i) Normal standard variates are then converted to flows by 
back trans format icai using the following equations; 



(3.9) 
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and 


Q.^ . = - q^ + Antilog 


(3.11 ) 


For incomplete records the following measures are 
taken: 

(j) For each station and month with incomplete record# say 
station 1# a search is made among all the stations to 
find a station# say station 2# which has a longer record 
and good correlation with station 1. The length of equi' 
valent record (N|| ) to be used from statical 2 to correct 
the statistics of station 1 is given by; 
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(k) Using this equivalent length of record of station 2 and 
corresponding statistics# the statistics of station 1 
are corrected by; 


X£ - = CX| -X 2 ) R S^/SJ (3.13) 

* (S^ « S|p (3.14) 

The primes indicate tl^ lc«g period values and without 
prime are based on the same short period for both 
statlcms 1 and 2 • 



(l) When inccxnplete records or gaps are encotjntered the corr~ 
elation coefficients computed in step (e) are not complete. 
The irissing values of correlation coefficients are esti- 
matfidt- by the equation; 

^in ^ \i. ^ ~ ■“ (3.15) 

In order to be ccaisistent with the two related correla- 
tion coefficients, the correlation coefficient must lie 
between the limits given by equation (3.15), These are 
established for all related pairs, and the average of 
ttese two limits is ta3<en as the estimated correlation 
coefficient . 

(m) All the correlation matrices are tested for tl«ir consis- 
tency (determination coefficients should be less than 
one). The test of consistency for each complete matrix 
is made by recomputing the multiple correlaticai coeffi- 
cient, If this value is greater than one, an adjustment 
is made by introducing a coefficient, successively smaller 
by 0,2, on tl^ radical in equation (3.15). All triaat 
consistency tests are repeated until all matrices beccxne 
consistent , 

3.3, Streamflcwr Gauging Sites of the Narmada River 

The Narmada river basin extends over an area 38,114 
square miles and lies between east longitude 72‘’-32* to 
81 *-45* and north latitude 21°-20‘ to 23 "-45', Lying in the 
Northern extremity of the Deccan plateau, the basin covers 



large areas in the States of Madhya Pradesh and Gujarat and 
comparatively smaller area in Maharashtra. The State-wise 
distribution of the drainage area is as under: 


1 ■" 

State i 

1 

... - - ,. f 

Drainage area in square miles 

Madhya Pradesh 

33#123 

Maharashtra 

593 

Gujarat 

4 #398 

Total 

38#114 


The basin is bounded on the north by the VindhyaS / cxi 
the east by Maikala range# on the south by Satpuras and on 
the west by the Arabian sea. The basin has an elongated shape 
as shown in Figure 3,3# with a maximum length of 592 miles 
from east to west and a maximum width 145 miles from north to 
south. The basin is divided into three zones nan^ly upper zone# 
middle zone and lower zone. The upper zone extends frcmr 
source to Bargi dam# middle zone from Bargi dam to Nammdasagar 
and the lower zone from Narmadasagar to Sardar sarovar. 

The Narmacfe river is the fifth largest river in India, 

The total length of jriver from the head to its outfall in the 
sea is 815 miles. The river has 41 tributaries^ out of which 
22 are on the left bank and 19 on the right bank. There are 
three main gauging sites an the main stream namely 
Ja^mtafa# Mortakka and Garudeshwar. Daily discharge observaticns 
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are taken at each gauging site. Monthly discharge data are 
available for 28 years at each ga\:^ing site. 

There are 30 major projects proposed for the whole 
basin, exit of these# five are purely hydropewer projects while 
others are multipurpose projects. Five major projects to be 
developed can the main stem of the river Narmada are: Bargi# 
Narmadasagar / Omkareshwar# Maheshwar and Sardar Sarovar. 

3,4. Streamflow Data 

As already stated there are three gauging sites along 
the main stem of the river Narmada, Jamtara is situated 12 
miles below Bargi dam site. It is at a distance of 247 miles 
from source. It covers a catchment area of 6400 square miles. 
Observations are carried out daily with current weter. During 
rainy season, observations are taken from a bridge on tte 
South Eastern Railway lirte from Balaghat to Jabalpur. During 
dry weather# the site is shifted 1000 ft dewnstream where the 
water flews through a narrow channel. Observations froit 1948 
arc available. Monthly streamflow data are tabulated in 
Table 3.4.1. 

The Mortakka gauging site is situated 37 miles below 
Narmadasagar dam site. It is located at a distance of 562 
miles from the source, Ihe catchrrent area up to this gauge 
site is 25,942 square miles. Daily discharge readings are 
taken with floats at a section 3000 ft downstream of a Railway- 
bridge on the Western Railway from Khandwa to Indore. Obser- 
vations from 1948 are available. The mcMnthly streamflow data 
are tabulated in Table 3.4.2. 
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The Garudeshwar gauging site is situated 7 miles 
below Sardar Sarovar. It is located at a distance of 730 miles 
from the source. It covers a catchment area of 34,496 square 
miles. Prcxn 1948 to 1961, daily discharge observations were 
made by float method with the help of a boat plying across the 
river. In April 1961, a high level bridge was completed 
across the river for the Eastern State highway and the discharge 
observations are taken with current meter. Observations made 
with float have been correlated subsequently with the obser- 
vatiCMis made with current rreter. The monthly streamflow data 
are tabulated in Table 3.4.3. For 3-hour ly flood routing 
study, four floods namely floods of 1959 (September), 1961 
(September), 1970 (Septerrber) and 1975 (Septenber) are 
selected as standard floods observed at Mortakka. and 
Gartjdeshwar gauging sites. The daily discharge data for these 
floods are tabulated in Tables 3,4.4 to 3.4.7. 

The monthly as well as annual statistical paran^ters 
like mean and standard deviaticn of the observed historical 
streamflow records at each gauging site are tabulated in 
Tables 3.4,8 and 3.4.9. The maximum and minimum values are 
also shown in the same tables. 

In the present study, five reservoirs namely Bargi, 
Narmadasagar, Orokareshwar, Maheshwar and Sardar Sarovar are 
considered. The Bargi dam site is situated at a distance of 
235 miles fraon the source. The catchment area up to Bargi is 
6090 square miles. The Jamtara is the nearest gauging slte- 
The streamflow data observed at Jamtara are transferred cn 
the basis of drainage area ratio. 
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TABLE 3.4.4. AVERAGE tAILY DISCHARGE FOR SEPTEMBER 1959 FLOOD 

(Plows in cf s ) 


Date 

i 

1 

1 

Name of gauging sites 

f 

1 

f 

Mortakka 5 

- - - 1 ... 

Garudeshwar 

1 


287907 

494377 

2 


184372 

407793 

3 


814865 

283674 

4 


334715 

978000 

5 


278235 

601528 

6 


245229 

418529 

7 


261855 

408532 

8 


143000 

307856 

9 


161956 

278964 

10 


289742 

229000 

11 


414175 

380527 

12 


380710 

590560 

13 


265844 

463728 

14 


303368 

355350 

15 


971385 

1175890 

16 


255537 

759816 

17 


180277 

327511 

18 


145930 

307714 

19 


112007 

247921 

20 


94886 

199169 

21 


80978 

144484 

22 


73671 

131700 

23 


66364 

120776 

24 


61140 

93271 

25 


57892 

98765 

26 


52350 

80121 

27 


47373 

60698 

28 


47478 

80800 

29 


52597 

69649 

30 


52279 

75904 

sum 


6718117 

10111750 


Average 
monthly flow 


13.325 MAP 


20.0516RAP 






TABLE 3.4.5. AVERAGE DAILY 'DISCHARGE FOR SEPTEMBER 1961 FLOOD 

(Flows in cf s ) 


Date 

1 

1 

Name 

of gauging 

sites 

I 

« 

I 

1 

Mortakka 

■■ ■' 1 “ 

t 

1 

1 

Garudeshwar 

1 


167216 


122300 

2 


140247 


181075 

3 


124538 


183559 

4 


241276 


128870 

5 


386112 


154409 

6 


367932 


505373 

7 


281906 


510105 

8 


194856 


410570 

9 


630459 


283130 

10 


1191380 


1189166 

11 


888854 


1376788 

12 


405103 


916681 

13 


209859 


425313 

14 


137670 


299027 

15 


111760 


214393 

16 


998213 


169560 

17 


1191020 


1445455 

18 


862344 


1410000 

19 


237181 


520642 

20 


177277 


292592 

21 


177277 


238289 

22 


137882 


218179 

23 


118396 


173610 

24 


111760 


158225 

25 


137211 


124669 

26 


184266 


231367 

27 


314876 


521724 

28 


317418 


545235 

29 


153520 


380210 

30 


175335 


282740 

M\m 


10773143 


13739178 


Average 
monthly flow 


21.368 MAP 


27.245 MAF 
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TABLE 3.4.6. AVERAGE DAILY DISCHARGE FOR SEPTEMBER 1970 FLOOD 

(Flows in cf s ) 


Date 


Name of gauging 


Mortakka 


Garudeshwar 


1 

369244 

278941 

2 

398813 

370650 

3 

437141 

373474 

4 

1055901 

540090 

5 

1097028 

1242913 

6 

1366240 

2153300 

7 

812125 

2252140 

8 

354049 

360160 

9 

200322 

508320 

10 

161061 

367120 

11 

134173 

223802 

12 

118336 

68482 

13 

183143 

153555 

14 

161208 

218860 

15 

128634 

190832 

16 

104578 

148260 

17 

84265 

127080 

18 

76535 

113172 

19 

68982 

104629 

20 

61114 

93439 

21 

52658 

78401 

22 

51565 

66964 

23 

82303 

91074 

24 

84265 

78366 

25 

80924 

104135 

26 

80924 

102370 

27 

119984 

103429 

28 

79583 

114372 

29 

59274 

101664 

30 

53198 

75189 

sum 

8117570 

11131027 


Average 
mcmthly flew 


16,100 MAP 


22.075 MAP 
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TABLE 3.4.7. AVERAGE DAILY DISCWvRGE FOR SEPTEMBER 1975 FLOOD 

(Flows in c£s) 


Date 

t 

1 Name 

1 , . 

of gauging sites 


1 MortaXka 

f 

1 Garudeshwar 


i L 


1 

124397 

169864 

2 

155779 

169652 

3 

171381 

262350 

4 

161992 

264503 

5 

222108 

208552 

6 

177382 

312758 

7 

157579 

236828 

8 

136399 

208446 

9 

136823 

169758 

10 

155955 

180277 

11 

328537 

369979 

12 

211870 

840271 

13 

1029350 

1189610 

14 

611750 

430413 

15 

272304 

297050 

1 6 

168698 

214800 

17 

131563 

156061 

18 

102546 

132234 

19 

88285 

91590 

20 

76035 

77766 

21 

67600 

64281 

22 

60998 

56586 

23 

59339 

50479 

24 

61316 

52986 

25 

76601 

60504 

26 

61916 

75274 

27 

59339 

61493 

28 

56091 

52279 

29 

57433 

46772 

30 

64600 

48114 

sum 

5309143 

6377803 


Average 
monthly flew 


10 J.34 MAP 


12.653 MAP 
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TABLE 3.4,8. STATISTICS OF STREAMPLOWS AT THREE GAUGE SITES 

(Plcws in cfs) 


r- 

1 

Gauge | Month 

site No. ! 
and name J 

t 

1 i . - . 1. .j.™ 

* 1 

; Flow 1 

J max 1 

t t 

i 1 

1 • 

i - - , 

"" - ' 

Plow 

min 

Plow 

mean 

T 

1 

; Standard 

; deviation 

1 

I 

1 

1 

P_ 

75% 

dependable 

flow 

July 

77169 

1480 

25999 

17335 


Axig 

87984 

8457 

49076 

17585 


Sept 

73506 

3479 

27254 

15663 


Oct 

21240 

846 

7375 

5726 


Nov 

9159 

370 

2052 

1731 


11 Dec 

4521 

309 

1031 

813 


Jamtara Jan 

2862 

211 

713 

815 


Feb 

1170 

1 44 

449 

256 


March 

846 

65 

314 

219 


April 

773 

17 

187 

171 


May 

211 

16 

61 

48 


June 

16805 

17 

1831 

3335 


Annual 

statistics 

242885 

37723 

116341 

44906 

81749 
(4,95 MAP) 



July 

239427 

24606 

78998 

51068 



Aug 

278345 

29469 

164555 

60629 



Sept 

362255 

25494 

140060 

84255 



Oct 

96847 

4960 

33284 

23792 



Nov 

25964 

2706 

10757 

6105 


12 

Dec 

10604 

2017 

6010 

2496 


Mortakka 

Jan 

7985 

1594 

4300 

1648 



Feb 

5492 

1224 

3268 

1029 



March 

4277 

911 

2392 

860 



April 

4201 

689 

1807 

829 



May 

2017 

585 

1176 

432 



June 

36131 

739 

6007 

7597 


Annual 


870938 

142398 

452613 

167573 

354130 

statistics 





(21.43 MAP) 


Note : The monthly and annual flews are average discharge 
rates. 
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TABLE 3,4.9. STATISTICS OF STREAMFLCWS AT THREE GAUGE SITES 

(Plws in cfs) 


1 

f 

Gauge | Month 

site No, ; 

and name * 

1 

f 

Plow 

max 

Flow 

min 

Flow 

mean 

Standard 

deviation 

75% 

dependable 

flew 

July 

284720 

25273 

82240 

55750 


Aug 

339348 

41374 

183691 

70456 


Sept 

467172 

33779 

180503 

116684 


Oct 

154874 

6570 

42613 

35229 


13 

Garude- 

29561 

2706 

12808 

7069 


14751 

1887 

6825 

2920 


shwar if* 

Peb 

13027 

1545 

4799 

2283 


9831 

1278 

3681 

1664 


March 

6310 

878 

2564 

1162 


April 

4353 

555 

1849 

887 


May 

3301 

358 

1212 

613 


Jiine 

60348 

471 

9952 

13798 


Annual 

statistics 

1047601 

165490 

532737 

213372 

379180 
(22.94 MAP) 


Note; The monthly and annual flows are average discharge rates. 


Marma das agar is located at a distance of 525 miles 
from the soxurce. The catchment area up Vo Narmadasagar is 
23,800 square miles. Ihe nearest gauging site is MortaKka. 

The inflow series between Jamtara and Mortakka are determined. 
The streamflows are transferred to Narmadasagar and Omkareshwar 
dam sites on the basis of drainage area ratio. A schematic 
sketch as shown in Figure 3.4.1 shows the location of gauging 
sites, dam sites and drainage areas covered by them. 

Omkareshwar site is 34 miles below Narmadasagar i.e. 

559 miles from the source. The catchment area up to Omkareshwar 
dam site is 25,050 square miles. The nearest gauging site is 
Mortakka, Maheshwar is located at a distance of 585 miles 
from the source. It is 23 miles below Mortakka gauge site. 

The catchment area up to Maheshwar is 26,712 square miles. 

The inflow series between Mortakka and Garudeshwar gauging 
sites are evaluated. The streamflcws are transferred to 
Maheshwar and Sardar Sarovar dam sites on the basis of drainage 
area ratio. The transferred streamflows data at each dam site 
are tabulated in Tables 3.4,10 to 3,4.14, The maximum, minimum, 
mean and standard deviation for each month are computed at 
each dam site for 28 years of record of streamflows and 
tabulated in Tables 3.4.15 to 3*4.17, 

It is found that the time lag observed between tte 
peaks of flood observed at Mortakka and Garudeshwar is 12 hours, 
A travel time of 15 hoxirs is assumed between Narmadasagar and 
Sardar Sarovar. The flood hydrographs for four selected floods 
observed at Mortakka and Garudeshwar gauge sites are plotted 
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Garudeshwar gauging site 

Note Numbgrs ore in sq. miles 


3 4 1 SCHEMATIC SKETCH OF DAM SITES , GAUGING 
SITES AND DRAINAGE AREAS 
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TABLE 3 . 4 . 11 . INFLOW BE'T/JSEN BARGI AND NARMADASAGAR BASED ON OBSERVED PLOWS AT JAMTARA AND 

NORTAKKA GAUG1!«3 SITES 
(Flews in cfs) 
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TAB1£ 3.4.12. INFLOW BETWEEN NARMADASAGAR AND OMKARESHWAR BASED ON OBSERVED FLOWS AT JAMTARA 
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monthly fl<»?s shown are average discharge rates during the month 


TABI£ 3.4,14. INFLOW BETWEEN MAHBSHWAR AND SARDAR SAROVAR B?^SED ON OBSERVED PLCWS AT MOiTAKKA 

AND GARUDESHWAR GAUGING SITES 
{Plows In cfs) 
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monthly flews shown are average discharge rates diiring the month 



TABLE 3.4.15 


STATISTICS OP TRANSFERRED STREAMPLOWS 
(Plows in cfs ) 


■f 

f 


station J Month 

; Plow 

: Flow 

{ Plow 

J Standard ' 

75% 

No. and J 

I max 

1 min 

; mean 

1 deviaticai 

dependable 

name | 

1 

1 ... - 

t 

t 

1 

-1 - 


t 

t 

1 

1 

flow 

July 

73431 

1408 

24740 

16496 


Aug 

83722 

8047 

46699 

16733 


Sept 

69946 

3310 

27595 

15360 


Oct 

20211 

805 

7018 

5449 


Nov 

8715 

352 

1953 

1648 


101 Dec 

4302 

294 

981 

774 


Bargi Jan 

2723 

201 

678 

490 


Feb 

1113 

137 

427 

243 


March 

805 

62 

298 

208 


April 

736 

16 

178 

163 


May 

201 

15 

58 

46 


J\ane 

15991 

16 

1742 

3173 


Annual 

231120 

35895 

112367 

43061 

777.88 

statistics 





(4.71 MAP) 



July 

188591 

11314 

48449 

37082 

i 


Aug 

195060 

19119 

105198 

45035 



Sept 

260660 

19770 

100292 

63115 



Oct 

70928 

3704 

23426 

17103 


102 

Narmada- 

Nov 

19425 

2061 

7850 

4750 


Dec 

8234 

1536 

4484 

1793 


■J sn 

5877 

1242 

3229 

1231 


SSQS'jC 

Peb 

4299 

969 

2534 

807 



March 

3223 

387 

1866 

657 



April 

3330 

599 

1451 

631 



May 

1700 

466 

995 

362 



June 

18199 

221 

3920 

4591 


Annual 

statistics 

641561 

98195 

303694 

120933 

230034 
(13.92 MAP) 


Note: The iTKjnthly and annual flews are average discharge 
rates . 




TABLE 3.4.16. 


STATISTICS OF TRANSFERRED STREAMFLOWS 


(Plows in cfs) 


Station I 
No. and ! 

name J 

1 

f 

Month 

Plow 

max 

i 

Flow j 

min 1 

1 

1 

t 

Flew 

mean 

; Standard 

I deviation 

I 

75% 

dependable 

flow 




July 

93447 

763 

3390 

2638 



Aug 

13803 

1344 

7387 

3205 



Sept 

18470 

1408 

7104 

4486 



Oct 

5046 

263 

1657 

1213 


103 

Nov 

1380 

146 

557 

339 


JL,. W 

CXfikars - 

Dec 

587 

109 

318 

127 



Jan 

417 

88 

230 

88 


'53 it j|W(f 

Feb 

306 

69 

180 

58 



March 

230 

26 

133 

47 



April 

238 

43 

104 

45 



May 

122 

33 

71 

26 



June 

1294 

13 

276 

322 


Annual 


45536 

6933 

21406 

8583 

16201 

statistics 





(0.980 MAP) 


July 

13447 

763 

2718 

2541 


Aug 

13803 

1344 

6994 

3443 


Sept 

18470 

1408 

8710 

6295 


Oct 

5046 

263 

2036 

2051 


104 

Qmkare- 

shwar 

3 380 

146 

584 

421 


587 

109 

302 

191 


417 

88 

210 

149 


306 

69 

167 

120 


March 

230 

26 

111 

76 


■A^^r i 1 

238 

43 

78 

56 


May 

122 

33 

55 

41 


June 

1294 

13 

480 

818 


Annual 

48399 

7009 

22446 

10501 

14501 

statistics 





(0.877 MAP) 


Note: The monthly and annual flows are average discharge 
rates. 
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TABLE 3.4.17. STATISTICS OF TRANSFERRED STREAMFLOWS 

(Flows in cfs) 


i 

t 

Station 1 Mc^th 
No. and I 

name J 

1 

- - - - - 1 

T 

Plow ; 

max 

L - . f 

Flew 

min 

......i 

Flow 

mean 

- -J 

Standard 

deviation 

75% 

dependable 

flow 

-...„ , *,... 1 . ,r- ....1 

July 

281835 

25340 

82037 

55346 


Aug 

335597 

40642 

182514 

69689 


Sept 

451335 

33270 

178016 

114594 


Oct 

151068 

6471 

42039 

34411 


, fvc Nov 

28624 

2714 

12682 

6945 


„ , Dec 

Sardar 

Sarovar 

14283 

1898 

6775 

2854 


12671^ 

1557 

4768 

2221 


9501 

1288 

3656 

1605 


March 

6117 

881 

2553 

1132 


April 

4222 

563 

1846 

872 

1 

May 

3201 

397 

1210 

596 


June 

58114 

487 

9709 

13344 


Annual 

statistics 

1036638 

164069 

527806 

210408 

377640 
(22.85 MAF)i 


Note: The monthly and annual flows are average discharge- 
rates. 
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using daily discharge data and they are shown in Figxjres 3.4.2 
to 3.4.9. The 3-hourly flood ordinates are obtained and tabu- 
lated in Tables 3.4,18 to 3,4.25. 

The flood ordinates are moved forward in time 3-hours 
for Narmadasagar . They are modified on the basis of mean 
monthly flow ratio. For Qmkareshwar^ the travel time between 
Omkareshwar and Mortakka is small and hence no travel time is 
alleged. The ordinates are however moderated cn the basis of 
mean monthly flow ratio. For Maheshwar the flood ordinates 
are first lagged by 3 -hours and they are moderated on the basis 
of mean mcwithly flews ratio. For Sardar Sarovar, the travel 
time between Sardar Sarovar and Garudeshwar is assumed as zero. 
The flood ordinates are however modified on the basis of the 
mean monthly flews ratio. 

The 75% dependable flow at Sardar Sarovar is 27.0 MAP 
which is decided on the basis of the record of annual flows of 
79 years. They are tabulated in Table 3.4.26,. It is menticned 
in the NWDT Volume 3 that the 75% dependable flow at Mortakka 
and Garudeshwar gauging sites is 22.0 MAP and 27,22 MAP 
respectively on the basis of 79 years. It is found that with 
28 years of record, the 75% dependable flow at Jamtara, 

Mortakka and Garudeshwar gauging sites is 4.95 MAP, 21,43 
MAP, and 22.94 MAP respectively. It is mentioned in the NWDT 
Volume 1 that with 22 years of record the 75% dependable flow 
at Jamtara, Mortakka and Garudeshwar gauging sites is 5.15 
MAP, 20.30 MAP, and, 22,60 MAP respectively. 
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TABI£ 3.4.21. 3-HOURLY STRSAMFLOW DATA AT MORTAKKA, 1975 (SEPTEMBER) 

(Plows In c£s) 
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M3 
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uo 

o 

03 

O 

uo 

cn 


M3 

cn 

O 

00 

M> 

rH 

M3 

xf 


o 

ro 

M> 

M3 

rH 


xt 

ro 

r- 

03 

o 

in 

in 

UO 

rH 

00 


M3 

:U0, 

M3 


m 

UO 

M3 

M3 


\f \ 

rH 

rH 

rH 

rH 

CM 

rH 

rH 

rH 

rH 

CM 

o 

rH 

M3 

CM 

rH 

rH 













TABLE 3,4.22. 3-HOURLY STREAMFLCM DATA AT GARUDESHWAR, 1959 (SEPTEMBER) 

( Plows in Cf S ) 


O o tn o l> o O O tH o OD 

OroOrH^^yicnOtnoin 
r^in'^OoiOv0CNic4or-r^ 
vD ^ rO CNJ in ro tH 


o 

00 

a 

a\ 

o 

o 

o 

m 

o 

o 

tH 

a 

lO 

o 

o 

o 

M> 

o 

o 

o 

o 

rHf 



tn 


o 

o 

00 

o 

a\ 

5 


ro 

m 

CM 


o 

r- 

a\ 



a\ 

a\ 

.r- 

lO 


CM 

rHI 

1-1 

iH 





^o O 


m 

in 


OOOQOOOOOOOOOOOOOOOOOOOO 

oooooooooooooooooooooooo 

oooooooooooooooooasvooonoo 

moOOOtnOOOOOOOOOxOoanroovowcM-'d* 

oovoc0'=i'c>jc'i\0ro0f^o0tr)crvcococj\-5j<cn(ncrir'vor~r' 


r'O^ooxooooOiHO'iJoOv^ 
I>Ov0OCStOOOCMOv£)Or-(OtHOww 
nOT-HOinO<Nr-LnOi>o cx)0':f0'=i’>H 

■<^incnoc0OC0rHOO<»^OcriC * 

0'^'C'3 0^^^tnLnoocolno^lnln\c 
'<#nfoci>ritfnoicNroinni>r'rr 


o 


o 

O 

o 

iHi 

o 

m 

O 

o 

CX) 

o 

in 

Q 

CM 

o 

m 

o 

o 


fH 

as 

o 


fH 

O 

o 

o 


ro 

ro 

m 

O 


CO 

ro 

o 


ro 

a\ 

as 

CO 

lO 



CM 

r-i 

tH 








O O O O 
O OO O 
O O O O 
CD o O O 

CO 0\ T-l ri* 


o o o o o 
O O o O Q 
OO ooo 
CM O o ID O 


O 

o 

O 

o 

o 

O 

o 

o 

O 

O 

O 

o 

o 

o 

o 

8 

O 

O 

O 

8 

8 

O 

m 

O 

o 

ro 


o 

CM 

r-i 

ro 

m 

CO 


a\ 

o 

vD 

ro 

o 

Os 

r- 

OD 

ro 

CM 

CM 

r-J 


r-l 




ooo 

ooo 

a\oo 

O ^ Oi 


Q 

tH 

O 

00 

o 

o 

o 

in 

o 

Q 

o 

lO 

o 

rH 

Q 

CO 

o 

o 

o 


o 

O 

Q 

as 

O 

m 

o 

I> 

o 

t> 

o 


o 

i> 

a 

\0 

O 

i> 

CO 

00 

o 

tH' 

00 

iH 

ro 

O 

vO 

tH 

9 

CM 

a\ 

o 


ro 

o 

as 


vO 


i> 

O 

ro 

CM 

CM 

tH 

rH 

rH 







OC^OOOCOOcMOO w w w '-u 

OOOOo\vDOOOrocDcr>fninocMcrvO 

■<a<'M*r0rH'M''^C'l(NCNLn'5jtr0OOCMcM 


OOOOOOOOOOOOOOOOOOOOOOOO 

oooooooooooooooooc-oooooO 

OOOOOOOOOOOOOOOOO'«DOOinOOC0 

OmOOOOOOOOinCMOOCMOOrM-^diCMCMtHOOO 

^~.^c^f^cDT^c3^■^■5}'vO'^tncMCoOcMCDfOt-|a^r'vo^'^^ 

■^•^CMrH'Cl'Ti<CM<MCMin'«l'Cn«H0OmCMrHrMT-l 


rni ’Sj’ w* V-/ 1 t 

fDOt^OroOiDOcM 

t-iO'OO'^OOOOLn 

coOfDOD'unr'OtH 


O 

00 

O 

as 

O 


O 

O 

O 

SO 

o 

in 

Q 

CM 

O 

V0 

Q 

rH 

O 

ro 

ro 

r- 

o 

CM 

o 

l> 

o 

lO 

O 

r- 

O 

as 

vD 


in 

as 

o 

ro 

ro 

o 

o 


O 

so 

rH 

O 

CM 
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ro 

lO 

m 
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O 

ro 

00 

rO 

CM 

as 

MO 
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CM 


ro 

CM 

rH 

rH' 
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O 

9 

9 
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o 

O 

9 

O 

9 

Q 

O 

o 

9 

o 

9 

O 

O 

9 

'O 

o 

o 

o 

'O 

O 

O" 

9 

O 

o 

O 

O' 

o 

5 

O 

o 

■Q 

o 

■o 

O 
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o 
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o 

o 

00 
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m 
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ro 

CM 

rH 

rH 
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o 
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TABI£ 3,4,23» 3--H0URLY STRBAMFLCW DATA AT GARUDESHWAR, 1961 (SEPTEMBER) 

(Plows in cfs) 


Q£?S?C>Qi/> 0000000-00000000000 
OCOOr-rO^Ot^OcnoOOO^OcrvOroOOOOOO 
0^0'=^ Q COO rHOiHO-^O roOCOC'i-'-tO r'0<^0-^ 
tr><TiO'»:l*OfHtn o\OcDO'4‘OOOc'ic'iooin^OCT»Lncvi 

ror^'cj'tn coo\*-i oOt^ntfitH^o r^'^cr\'<:iiOvO cvr^iOinco 

r-ti--»rHT-l^^ni-HC>l|>nOItnLfyOOOJ(NCNli-!i-|CN|'^'^<N 


OOOOOOQOOOOOOOOOOOOOOOOO 

OOOOOOOOQOQOOOOOOOOOOOOO 

Oi>OOOOOOOOOOOOOOoJOOOOOOO 

ocnxofooiiioooomooCTvoooococrvoooooo 

cn CT\ fSi ovcnro-^cMCTvc^r-ivoOOtHrj'Ovo c'imr~a\C7v 
T-lT-^r-lr^•^'^^OT^C^^^nOInlr)OC^C^3CMr^T-|C'l-!a^^CM 


OOOOiTYOOOOOOOOOOOOOOOOOOO 
ooi>0c~-oooco0'^)'0vo0ir)oa\0r-j0c^or'0 
n^v00ro0rt*0r^0cnoui0o\0 cMO^OOfOOi^O 
r'i0<Nc«u*0v£>0^0 C0O0\C0cNif0a3rHr0(MrHin'=3'O 

C'3T-tiAiooO'^ot>oncDC'i\OCO\o^cOrHr~foocr\r'JO 

THiHr-ftHiAinrOTHfOCOroc'JiH^T-fro r^oJ tH tH c> 3 ro 


QOOOOOOOOOOOOOOOOOOOOOOO 

OOOOOOQOOOOOOOOOOOtnoOOOO 

OjiDOOOOOQOOOOOOOOOOcnoOOOO 

rHOinoinintnO'd'Otnooor'OCpLninfOvDOOOO 

C<^C0v0lf)00O'X)'MO^~0^^0^0^0f0vDTj^r^C^nf^r^r^^^ 

fHl r*i 


Ou^QQOQ^OOOOOOu^OOOOOOQ^nOLn 

Ot^OCOOOOOOOOOcooinOnOcoocTiOf'aoco 

OOO<NOr^Q0^Ov0OC0O^O^0Ow'<i<,^Ot>5^^ 

0^r^^00^ooOl>cM\Olr)^0l>tnO(nOcDnOO^^unO 

rHCDr''^r',-iCOvocM.HOrovO'sa<o\crivn,Hi>-^oJcMfOr-f 

t-)rH*-)tH'«l<inrOCJra<Cr\'?i«CSi<H'<^C'JnCMCMrHr-<r'lininrO 


gggooooooooooooooooooooo 

RRSS9RSSSS929§oeoooooooo 

ogoooooQoooooooooomooooo 

r'OinnininLnOOOtntn\ooo‘Jiu>von\Oinooin 

r^^^c^Ti^^/l0^c^c^rl<a\^lnvo■<;J‘Lnc^3Lncll^'■<;J•T^vO•cjt^o 

«-<T-l,-(,-|'ct'^rOC'J'^Cr»-«i*Cl.-<fOCO>s}'CMC4rHrHCM'4'Un(r) 


ogooogooooooooooooooootnoo 

ngvog'<4*oog.HOrHOa\ooo-:j<oooLnonoo 

t-40inOrHOir)OCOQroOnooOtr>0>ooi>OcsioO 

unoncot^c'jOuncootnoLnoooocMfnoOOinoo 

T-4'OCOfOcooOt-ir'invDc^'i?>i)C'j’HvD'£iroi>ijr)*H*-i'^inr^ 

r-<<HrH>H'^'<:l<-«a'CS'«a'0'^C4i-(OJ'«a'<j'C'J<NrH»-|CN^inoOCN 

r-t iH r-i mi 


OOQQOOOOOOOOQOQOOOOOOOOOO 

ggggoggggogoogggooooooQOO 
ogoooooQogQooooooomoooooo 
oog'^oooor^oomcogoo oocnr^ryooomo 
r^cnCDoOv£>C~•^OCO•-^fO(^Jr^C^O■<;CC7^^)fOCOtnO\vOC^JvO^~' 


S cnOOOLOOO'^OOOOmoOOOOinoOOOO 
OOOi>OroQoovOOQfOOl>OsJ'0'<#OoiO OvOt-io 
OvOOCOO C0Or-^c0r^OOO COOvOOxj.moiOT-iOC'JO 
tntn(NCDOOOOOO<^)COOC^O'<J'OO^J^\OOOCDOC^OOu^ 
O-^foooicriininoocomaiCRmcor^cir-TfOOtnr^c-iOcor^ 

tHI rHI «H 


g OOOOOOgOOOOOOQOOOOOOOOOO 
oggggogoogoogogoooooooQo 
ggoogSgQoggoogggoMtnoooogo 
OO^'ii'OOOOOOOirjtNgOO Tji tA H o o o o o 
Ot»«oOfoc'4t^i>QaKnvoc4O®Ococo««#(T»\0inav 03 cn 00 
»-( rH CM Xt nP in CM M vO fO CJ r' m vO C4 (N <-l iH Cl -.S' rt <M 



TABLE 3,4.24, 3 -HOURLY STRSAMFLOW DATA AT GARUDESHWAR, 1970 (SEPTEMBER) 

(Flows in cfs) 


OOQroOOOOO^OOOOQOOOQinOvOOO 
O-^O-r-fO -^O C^^O'OO^0O vDO -^ocOOfOO COOOv 
incMOCJ\ 0 ' 3 *Or-)OcaOCDO< 3 ' 0 ':tO' 5 i'LnTHO'^ 0 «-« 
rHr-iOiNO»-»Or'Ln<MinO'<d'Cj^nr'iC'>c»^'5j«c>icr\Ln 

CDOo^■^OOo<N\Or-(^-ooo^'^T^Oc^^'^0r^ooc^lOr~ 

CNjrO's^J'Olfnin^DnCNir-irH.Hi.HfHTH r-lr-< r^ 

t-l OJ tH 


oooo o o o o o o o o o o o o o o o o o o o 
oooooooooooooooooooooooo 

OOOOOOOOOOOOOOOOOvDOOOOOO 

Oo\oooooooOincnvDO(X)^inpocr\,Hnunfno 

oor'r'ir~Ocor^r~'Oiinvoc^’= 3 *cNO<T'r'\Dr'OOa\r-ioo 

^J^ 0 ■ 5 d*O<N^>v^)fnC^^r^r^r^r^^^r^ iHrH r^ 

T-t <N r^ 


OOOOOOOOMOOOOOOOOOunOOOvOO 
'?i’ 0 ' 5 f 00 oooooooOvi)OvoooocMOi>oa\o 
CT'QcaOcoOtHOCOO cMOcNOtnO^r^-^OrOLniHin 
oov 0 ,HOnot^Lnncoc^vOCDcMOcocO'^avc 3 \c^r'r~cM 
^'^~CDC'JtnmMr'<>^o^l'^OTl'CMr^o^^-vDr^a^oa^r^oo 

(N rH 


oooooooooooooooooooooooo 

oooooooooooooooooooooooo 

OOOOOOOOOOOOOOOOOCDOOtnmOO 

OfocriooooooooOroO'^<HOOi^<Nr~oior~in 

(Njr^r~r^Ooococoino\ir>T-)vOcMrHOCD\ooo<T\OOHOO 

nr 0 r 0 r~rHOr^r 0 C'I iHOJi-tT-HiHr-i tHi— trH 

(N CM 


8 oooqoooooooooo'=i' 0 ^ooooo'^ 

LnocoOM'or'OcoooooooovoOvDOrHOno'^ 

OvDOfMOr-fOr'OM'OcnooocrvOcriOfncOM'Oin 
Oo^'r^Oc^JO•^OoOQOOO^'r^<Nf^^oLnLnc^lnvD^~■ 
tnr~r^r»Oin'd'COi>vOincNi>cM.-<OCDvOCOCP\OOrHCO 
rOrorOOOCMCOrOCM r-4>-lr-t 


OOOOOOOOOOOOOOOOOOOOOOOO 

oooooooooooooooooooooooo 

OOOOOOOOOOOr'OOOOOCOOOOOOO 
OOOOOOOOOOOOOOcMfOOOOOOfnroinc^ 
CO':J'^'OOC^^Or^C 3 ^ 00 cMCOnr-^OCO^OCOC^OOr^OO 
^OfOnOO^C'JC 7 VM*CM>^r-^C^^^^T^r^^^ rHi-lrH 

rH <N 


OO-^l'OOOOOOOvOO^OOOOO^OOQOOO 
oOr'O ooOvoO'i'OinOrHOr'OvDOr'OroOr'OO 
oO'^OoOT-iOvoOLnunr'Ot-iinOcnoO'snonoo 
csOfnoooooOOfooooonfnChOrHvOcn'gtrtror' 
f-jCMOCOr4-d<vD'<i'T-lOinCMCO<n,-lOCOvO<TiCOOQi-lCr>l> 

■ 5 j'fnrOOl>rOaV'^fOt-tr- 4 rMr-trHiHTH 
1-4 CM 


OOOOOOOOOOOQOOOOOOOQOOOOO 

QQOOOQOOQOOOOOOOOOOOOOQOO 

OOOOOOOOOOOcnOOOOOcnOO^OOOO 

OOoOOooOOcMaDOooon 3 '»-i\oncM(ninoi >'0 

,--lOt~-V 0 OCr>rOI><nOOinCMC»'cd«rHOCJ\vjOCOCOOOrHa»r> 

r-l N t-4 


OcmOOOOOOOOOOOOOOOOOOO'OQ'^O 

Or-4OCr\OcMOcMOCnOv0OOOcnOOOr'OOOvDO 

O-^OOOCOOrooOOOOocsicnvOOCnOfncOvDOOO 

in-<d'C^OOa\OOOOr^MCOCOH 01 '«l'CMvOvOCOfOv 0 ^t-(vO 

Ott>r'''d»O'^®O-^>-i' 0 iHC 0 - 5 i'«-to<TV' 0 r~’r^ooOoi> 


W I t. w w ^ ^ w 1 I > - 

ca 01 tn ^ 04 in CO Oi fH «H T-l H T--i 
^ CM' fH' 


' ^ir|' 


8 comoOOOOOOg:>ioOcMcMm<o 0 \oicx) 3 tnaounir 
a)r^OOOcnir^airHv 0 OO^*<tr-iOcn\£)vor'OOa\or 

^ ro 01 01 ro CM -CM r-4 H rH'. . r-l f-4 t-4 

fH ' CM tHI 



TABIE 3,4.25. 3-HOURLY STREAMFLOW DATA AT GARUDESHWAR, 1975 (SEPTEMBER) 

(Plows in cfs ) 


O'a’OOOOOOOOOOOOQOOOOOOOOO 
O COO^OOtr)OQOOl>OmOOOLnOtH 

c^Jnou^lOr^OC'JO'5^«OOO0\O^'O0^OLnO'ia'C0■-f 
c5%votn ooonooO'5}*o^'O■<^tf•l>a)(sooo^^ln^0CO 
v0un^0OOOC0C0for^'?^0\a\'<4tor~mmtAv0r^ir)':j*'«a» 

.-^lHO^C>Jnolr^r^nOvOOJt- 4 »-|»^ 


OOOOOOOOOOOOOOOOOOOOQOOO 

OOOOOQOOOQOOOOOOQOOOQOOO 

^ocNjoooooooooooooOinooointno 

chOtnooomcDooooor^rHCOtHnrHOf^vor'Ch 

vOvOvOt-ir-^t-fCor^m 

tHiHOJdfOOJtHT-inOr^Ol OlTHt-t 


OOOOOOOOOOOOOOO OOOOinOmO 
voocMOvoO\ooooo<sooo,HOooocr»or-OvoO 
coo-s'Ot^or'OcrvOvOOOoOLnOC'JLnroOcsioooin 
o\rotnir)CvioovDcr\0<c0^0r~cM^'^'^i*a\tjr)a)t>o\ 
'0'OvO'-i.-tc4a\r'vX)vO'd'niHint-i cDvOinvOint^in-^-^i' 
^^r^C^lr'l^OCMT-^r^^OO^00OOC'lT-lT^I> 


OQOOOOOOOOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOQOOOOOOOOO 
OOOOOOOOOOOOOOOmOLnOOOLfiOO 
m^omooO'O'^OOOOOmO CDr-inocDcJc^cocn 
caoncnr'OJO<nur> r'JLncNir~\DtnvOinc^Ln'<i»«;i« 


O'OOiDOOOQOO <MOO\QOO COOmOCOO'^OCM 
ocr\^nQ\00(NinoO\OOvOir>coo\'i>r'vOO\i-(crvo^ 
cn'iJvo (N©fnoi>i>'5l*<-ii>fOincsjt^vOinyiLnvO\o^i<^ 
^-^r^C^JC^^CMC^^C^Jr-^'NC0Or0C^lr-^r^ 


qoooqoooooooooqoooqqoooo 

oooooooooooooooooooooooo 

O'S'OOOOOOOOOOQOOOOOOviDOOOO 

mo\fnOOO':i'rHOOomoc'JOca»-iLn^in\onOcri 

OvOvorovOunoc^f^t-tcrvCOTjivonoor'inLninvDvom'-^ 

CN,HCMCM<NCMCMtH<Nt~-Of0O3THT-1 


vOOOOOOOOinOOOinoOOinOOOOQOOO 
COomO'^OinOoJO<HOMO<nOoiOCOOOQOOOO 
OC'if0O<NO'ttr)0\O'DQr0O<NiO'<;fm'<d'C0inOC'JOO 
r'aiC'JO':i*incooooocr»0'«3’COCMtn<nroO'<^<(n'<3*csicoo 
cs[vOvO';l'xi<v£)Or'CX)r'OOOinvD<ncor'Uiir>invOvomxj<in 
cN»-4r'irsJr'>r'JOJiHrHLn<H'«j' ^Jr^r^ 


00'OOOc4 00voOOir»o->^»noovn«NrM'<#<ntfjfor>a’ 
r^^O'^3•tn^ncDo^'COro\o,^^or^n cor'inintnioxom-^'^ 

CNrHCMCMlN<Nlr'J*-<r-('^t-|rJ<CSlrHT-t 


OOOOOOOOOOOOOOOOOOOOOOOO 
OrHOOOino^oOr^Oi-iOtno axocDoaiOfoor' 
Ocno‘nO'?i'Or'OinO'^0{MOinO'^mo\OonOr' 
O CDO'<J'OvDT-)a»'it'5l'OOOmOO«-<\OOr-lCMoiyD>«a'<£> 

Ov0r-<'0CNicT\O'DC0Otn<nr^r^r><Tir>-Lnininv0'om-<^ 

CM,-^C^^<SCNI01C^^r^.Hnr^•<#Mr^H 


ooqoqo qqoq 


Uiw IXUI ^NV«.,»WWWWr-l VMWI n WIUI vu W 

os'X>oO'0<-iO'Or*oo«HO'ti' cooo«<»a»r'tnunin'Ovoin^^ 



TABLE 3.4.26. 79 YEARS ANNUAL SERIES AT SARDAR SAROVAR DAM SI' 

ARRANGED IN ASCENDING ORDER 


SI. ; 
No. : 

t 

L. 

Water 

year 

- 

r 

Inf lew 1 

MAP { 

1 

1 

SI. 

No. 

I 

1 

f 

t 

1 

i 

1 

t 



Water ; 

year 1 

f 
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CHAPTER 4 


OPTIMIZATION MODEL 


4.1, Introduction 

Various optimizaticai models have been described in 
Chapter 2. As mentioned therein optimization models can be 
used to screen alternatives. Decision makers can examine a 
wide range of alternatives and select the most desirable caie. 
The discontinuous linear programming technique is selected for 
application to the present problem for screening the alter- 
natives. 

4.2. Formulation of the Model 
4.2.1. Physical system 

For laroper utilisation of water resources of the 
Najonada river basin a master plan was prepared by Madhya Pradesh 
in 1972, in which 32 reservoir sites were proposed as shown in 
the report of the Narmada Water Dispute Tribunal (NWDT), Voltome Vj 
Out of these , 31 reservoir sites were in Madhya Pradesh and one | 
was in Gujarat. The tribxmal has awarded 18.25 MAF of water to j 
Madhya Pradssh, 9,0 MAF of water to Gujarat, 0,5 MAF of wciter to ; 
Rajasthan and 0.25 MAF of water to Maharashtra. Gujarat has 
planned to utilize its share of water by building a reservoir, 
namely Sardar Sarovar at Navagam. This reservoir will also 
providtei Rajasthan's share of water. Madhya Pradesh has planned i 
to build several reservoirs across the main stream, and across 
the trlbxitaries , A few of those are listed belcw; : 




Reservoirs on tributaries 


Reservoirs on main river 

(a) 

Bama 

(a) 

Upper Narmada 

(b) 

Tawa 

(b) 

Ragavapur 

(c) 

Chhota tawa 

(c) 

Rosra 

(d) 

Dudhi 

(d) 

Basania 

(e ) 

Kolar 

(e) 

Bargi 

(f) 

Man 

(f) 

Chinki 

(g) 

Jobat 

(g) 

N armadas agar 



(h) 

Qmkareshwar 



(i) 

Maheshwar 



(j) 

Sardar Sarovar 


There are 10 major reservoirs planned «along the 

main stem of river Narmada. Out of these, 4 are purely hydro- 
power projects and remaining six are multipurpose projects* 
As already stated earlier, the Narmada river basin is 
divided into three zones. The upper zone extends up to Bargi, 
middle zone from Bargi to Narmadasagar , and lower zone below 
Narmadasagar, All the utilisaticxi upstream of the Bargi reser 
voir are assumed to be clubbed into the Bargi. The water 
requirements at and above Narmadasagar are to be TOt by the 
Narmadasagar reservoir. The water requirements to be met by 
QmJcareshwar and water requdrements up to Sardar Sarovar are 
clubbed at Ond<areshwar. Maheshwar is planned as hydropower 
3 ?ro ject , 

It is very difficult to acccmmodate all the planned 
reservoirs in the formulaticxi of linear programming model, 
because of the large size involved* It requires not cxily a 


large amount of ccmputer time but also adds canplexity in the 
formulation. Therefore it is desirable to incorporate a small 
number of major reservoirs into the model, and to club several 
smaller reservoir into a larger one. 

From the hydrological and agricultural considerations, 
the water year is considered from 1st of July to 30th of June 
of next calender year. Keeping in view the above mentioied 
facts five reservoirs and four seasons in a year are xised in 
the formulation of LPD model for the Narmada river basin londer 
study. 


4.2.2. System representaticn 

The schematic diagram for the LPD model is shown in 
Figure 4.2. The descripticai of the model is as follows; 

(a) Numbers 1, 2, 3, 4, and 5 represent the reservoir sites j 

(b) f t inflow due to the independent catchment 

area of reservoir site i, during time t; 

(c ) is the annual irrigatican target at reservoir site i? 
<d) Kj *. I, is the proporticn of irrigation demand I at 

J« I* v» 

reservoir site i, during time t? 

(e) ^ is the return flow from irrigaticai utilisation 
at reservoir site i, during time t? 

(f) the release to meet irrigaticn requirements at 
reservoir site i, during time t; 

ig) c, . indicates release to meet energy requirements at 
reservoir site i, during time t; 

(h) shows the firm eneifgy target at reservoir site i. 




FIG. 4.2 SCHEMATIC DIAGRAM FOR LPD MODEL 
FOR FINAL STAGE OF DEVELOPMENT 






4.3, Demands 


4.3.1. Irrigation 

The term ’water reqxiirement ’ and 'water demand* will 
be used interchangeably to mean those quantities that a given 
water resource project is expected to deliver iinder specified 
condition of allowable shortages. The types of major water 
demands most ccanmonly found in developing countries are demand 
for irrigation and power generation, A deterministic approach 
has been used for their estimation. 

In order to estimate irrigation water reqtiirements ^ it 
is necessary to loicw the proposed cropping pattern in the 
ccKnmand area. After taking into consideration the existing 
and the proposed cropping pattern, the NWDT award has decided 
for major, medi\»n, minor, microminor and pumping schemes 
the following deltas: 2.56 ft for major, 2*07 ft for medium, 
1.89 ft for minor, 1.50 ft for micrcxainor and 2.56 ft for 
pumping schemes at the head regulator. The areas to be covere<i 
\jnder each are taken fran the data available frcxn the report 
of the NWDT Volume I. The water requirements for irrigaticai, 
and industrial and domestic water supply at each reservoir 
site are given in Tables 4.3.1 to 4,3,4. 

4.3.2. Hydropc3wer 

The energy requirements are calculated on the basis 
of proposed installed capacities of turbines and load factor. 
The load factor is assumed as 0,60* The actual energy prod\ice< 
is given by below: 


TABLE 4.3.1. MCMTOLY DISTRI BUTI CW OF WATER REQUIREMENTS AT AND ABOVE BARGI DAM 
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TABI£ 4.3.2. MC»ITOLY DISTRIBUTION OF WATER REQUIREMENTS AT AND ABOVE NARMADASAGAR (N.S, 
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TABLE 4,3.3, MCWTHLY DISTRIBUTION OP WATER REQUIREMENTS aT AND BELOW OMKARESHWAR BUT UP TO 

SARDAR SAROVAR 
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TABIiE 4,3.4. MCOTMLY IXtSTRIBUTIC^ OP WATER R5X}UIREMENT AT SARDAR SAROVAR 
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4.3.3. Industrial and domestic water 

The quantity of water requirement for industrial and 
domestic water supply is taken frcxn the NWDT, Volume I, Nearly 
uniform distributicai throughout the year is assumed and is inc- 
luded in Tables 4.3.1 to 4.3.4. 

4.4. Objective Functicxi 

Many theoretical and practical approaches have been 
proposed for identifying and quantifying objectives in water 
resources planning. It is a tedioxis and frustrating process to 
incorporate them into multiobjective planning models. It is a 
challenging job for the system analyst to reccxnmsnd a single 
planning alternative for adoption (Loucks, et al./ 1981). 

Marglin states, "In view of the three dimensional nature of 
national welfare — the size of the economic pie, its divisictfi, 
and the method of slicing — we believe it xinwise to attempt to 
define a single index for the broad objective; instead we shall 
develop alternative objectives for the most important ways in 
which water developnrient can contribute to national welfare. 

These are broadly two: efficiency, which expresses the objec- 
tive of optimizaticsi of national income, and income redistri- 
butics:!^ which expresses the objective of achieving a desiired 
slicing of the ecoriomic pie by a knife that suit community 
values" (Maass, et al,, 1962). 

In the present study, the efficiency objective 
has been used. This objective requires an evaluati<»i of costs 
incurred in building and operating the system, and benefits 
reaped from the system. 



4 •4*1, Construction of cost and benefit functions 


Evaluation of benefits from a system is a tough task 
for the systems analyst. For the estimation of benefit from 
irrigation, the proposed cropping pattern, command area and 
price of irrigation water (Rs./acre foot) are the basic require- 
ments. Similarly projected energy demands and unit price/kwhr 
are necessary for determination, of gross benefits from energy. 
Du« to lack of economic data of the projects, the following 
procedure has been adopted to develop cost and benefit func- 
tions. Bargi project is chosen to illustrate the procedure. 

The Tribunal has assumed 45 years for the full deve- 
lopment of the Narmada basin. It is suggested by the Tribxinal 
that the award be reviewed after 45 years. The basin is planned 
to be developed in three stages. The first stage extends up to 
first 10 years. The second stage is from the 11th year to 30th 
year and the third stage from. 31st year to 45 year. 

A time horizon of 70 years is assToroed in the analysis, 
including construction period. A discount rate of 5% is 
assumed* The operation, maintenance and replacement (OMR) 
costs for irrigation are assumed on a unit basis, i.e., R\^es 
per unit area to be irrigated at different reservoir sites. 
Similarly OMR costs for power are asstamed on a unit basis, i.e.. 
Rupees per tanit Installed capacity of turbines at varioixs 
reservoir sites* 

The rate of development of irrigation is represented 
by Figure 4.4.1. Since the costs and benefits are non-uniform 
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over the time horizon^ it is necessary to use one of the three 
commonly used methods namely the present worth method, annual- 
cost method, or the internal rate of rettirn method, for ccOTpa- 
rison of costs and benefits. In the present analysis, the 
annual-cost method is adopted. The details are given below: 

(a) Annual benefit frcan irrigation; See Figure 4,4*1. 

ABF « [j3,2(P/A, 5, 10) X (P/F, 5, 6d) + 0*03(P/G, 5, 10) 

X (P/F, 5, 6>) + 0.5 (P/A, 5, 20) x (P/F, 5, 16) + 

0*02 (P/G, 5, 20) X (P/F, 5, 16) + 0.9)P/A, 5, 15) 

X (P/F, 5, 36) + 0.0066 (P/G, 5, 15) X (P/F, 5, 36) 

+ 1(P/A, 5, 19) X (P/F, 5, 5lJ X (A/P, 5, 70) | 

where ABF stands for annual benefit factor for irrigation, 
P/A for series present-worth factor, P/G for uniform- 
gradient- series present-worth factor, P/P for single- 
payment present-worth factor, and, A/P for capital- 
recovery factor. Using 5% discount rate and the rate of 
develcjpment as shown in Fig\ire 4,4,1, ABF works out to 
be 0,46>* A delta of 2*57 ft at the head jregulator and an 
intensity of irrigation of 135% are assumed in calculating 
irrigation water requirements. The long term irrigation 
benefit is talcen as Rs, 418/acre irrigated. Therefore 
the benefit per acre ft of water will be 


418 

i:"5T‘xt:is' 


m 


RS, 120,4a/acre ft or 
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6 

Rs* 120.48.x 10 /Haf. Annual benefit fronrt irrigation will 
be Rs. 120.48 X 10^ x ABP per MAF, i.e., Rs. 10 x 10*^ x 
0.46 a Rs. 5.546 X IoVmAF. 

(b) Annual cost of irrigation works: 

Irrigation cost for continuing projects is taken from the 
Joxirnal of Irrigation and Power, Volume 3, 197.8, p* 98, 
as Rs. 3360/liectare for Madhya Pradesh where Bargi is 
situated. Assuming a \inif orm. ejqpenditure over 20 year 
development period and a time horizon of another 50 years, 
the annual coital cost recovery works out to be 
RS. 1.2672 X 10*^/HAF. 

(g) Operation, maintenance and replacement (OMR) cost: 

OMR cost is assumed as Rs. 6/acre, i.e., Rs, 1.734/acre ft, 
or Rs. 0.1734 x lo'^/MAF. 

(d) Reservoir cost: 

For Bargi reservoir, the capital cost is 112.5 crores. 

The gross storage is 3.28 MAP. 

Capital cost/MAF » ^ Rs. 34.30 x l0*^/lyiAF 

Assupne a construction period of 6 years, discoxjint rate 
5%, and, total time horizon of 7/0 years. Using these. 


the annxjal capital cost recovery works out to be 
1.51 X IoVmaf. 

(£,) Annual cost of energy; 

At present there is a provision for producing energy from 
Bargi dam by providing canal head power house (CHPH ) . 

Proposed capacity of left side CHPH; 2x7. 5 = 15.0 MW 

Capital cost chargeable to power = 37.5 crores 

Cost/kw installed = ^o65 ' ~ 2500/kw 

Cost/lO^ kw = Rs. 250 X 10^/10^ kw. 

Assuming a uniform expenditure over 6 years of installation 
period, and, a time horizon of another 64 years, the annual 
capital cost recovery for hydropower works out to be 
0,001296 X loVlO^ kwhr. 

ijf) Operation, maintenance and replacement (OMR) cost; 

Ca® cost is assumed as Rs. 25/kw installed, i.e., 

Rs. 2.5 X 10^/10^ kw or 0,0002854 x loVlO^ kwhr. 

(|) Energy benefit; 

The rate of cte^relopment of energy for Bargi dam is shewn ! 

in Pigxire 4.4.2. Using a discount rate of 5%, and, a i 

envelopment period as shown in the figxare, the annual 

benefit from sale of energy at a rate of Rs, 0.50/kwhr 

7 '6 

comes to Rs. 0.0328 x 10 /lO kwhr, f 

Thus for Bargi project, it is seen that the annual | 

7 ' i 

gross benefit from irrigation is Rs. 5,5460 x 10 /MAP, anniial 

■ - ! 
7 ^ ' 

cost for irrigation Rs, 1,2672 x 10 /MAP, OMR cost for irrigaticai 




7 

R8m 0.1734 X 10 /MAP and annualized reservoir cost Rs. 1.51 x 
7 

10 /Map. Also gross annual benefits frcxn power is Rs. 0.0328 x 

7 & 

10 /lO kwhr/ anniial cost of power generation Rs. 0.001296 x 
loVlO^ kwhr and OMR cost for power Rs, 0.0002854 x loVlO^ Icwhr. 
Using each of these xinit costs and benefits as a guide# cost- 
capacity curves, gross benefit curves, OMR cost curves, etc. are 
prepared and shown in Pigvires 4,4,3 to 4,4,9 for Bargi project. 
Cost and benefit functions for each project are tabulated in 
Tables 4,4,1 to 4,4,5. The net benefit function for irrigation 
for Bargi project is developed and shown in Piguire 4,4.10, 
Similarly for net benefit function for energy for Bargi project 
is developed and is shown in Fig, 4,4.11, For other projects, 
detailed calculations, and curves sucdi as cost-capacity, gross 
benefit curves, OMR cost cinrves and net benefit curves are given 
in -^pendix, 

4,5, Mathematical Formulation of LPD Model 

The objective of the problem is to maximize the net 
annual benefits fromi irrigation and energy for the whole system 

5 

Maximize R « £: INB, ^ + ENB, . - , (4.2) 

The first subscript, 1, 2 and 3 represent irrigation, energy and j 

reservoir respectively. The second subscript i, represents a 

site in the system. The terms used in the objective fxinction are; 

■ ' ' ^ I 

INBi ^ *■ Net annual irrigation benefit at site i and is equal to [ 

- =1,1 - °i,i 

IB^ . » Gross annual irrigation benefit at site i; 

1 If JL " ■ ' ' 

« Annual cfi^ital cost of irrigation works at site If 



Annual cost of reservoir in li’ rupeii^ 



FIG. A.4.3 ANNUAL COST OF BARGI RESERVOIR 
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OMR cost of irrigation in 10 rupees 



Annual irrigation in MAP 
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FIG.4.4.9 OMR COST OF ENERGY FOR BARG! RESER- 
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TABLE 4-4.1. COST AND BENEFIT FUNCTIONS 


BARGI 


SI, 

No. 

1 

1 

; Independent 

I variable 

t 

- - - 

1 

t 

J Dependent 

J variable 

1 

- 1 ... - - 

1 

1 Functional relaticnship 

1 

1 

1 

- 

1 

Reservoir 
capacity (MAP) 

Capital cost 
in crores 

C 3 ^ = 1.615 - 0.073 yJ 

2 

Reservoir 
capacity (MAP) 

OMR cost in 
crores 

1 ~ 0*000 

3 

Annual target 
output for 
irrigation 
(MAP) 

Gross irriga- 
tion benefit 
in crores 

IBj^ ^ = 6,01^ - O. 421 J 

4 

Annual target . 
for irrigation 
(MAP) 

Capital cost 
of irrigation 
in crores 

C^ ^ = 1.341^ - 0.0561^ 

5 

Annual target 
for irrigation 
(MAP) 

OMR cost in 
crores 

0 ^ ^ = 0.17341^ 

6 

Annxjal target 
for energy 
( 10 ^ kwhr) 

Gross energy 
benefits in 
crores 

EB 2 jL = 0.03433Ej_ - 0,000029 eJ 

7 

Annual target 
for energy 
( 10 ^ Kwhr) 

Capital cost 
in crores 

!i 

C 2 ~ 0.001296E^ 

8 

Annual target 
for energy 
( 10 ^ Icwhr) 

ca*iR cost in 
crores 

©2 1 = 0,0002854Ej^ ; 


7 

Note ; 1 cxore = 10 . 
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TABLE 4.4.2, COST AND BENEFIT FUNCTIONS 


NARMADASAGAR 


r 

SI. I 

Independent 

"HI j 

1 f 

; Dependent J 

Functional relationship 

No. : 

1 

variable 

1 variable J 

i • 

-J L. 



Reservoir 
capacity (MAF) 

Capital cost 
in crore s 

C^ 2 = 0.486 - 0.0138 

Reservoir 
capacity (MAF) 

OMR cost in 
crore s 

0_ - = 0.000 

Annual irriga- 
tion target 
output (MAF) 

Gross benefit 
frc»n irriga- 
tion in 
crores 

IB^ 2 ” 0*1461^ 

Annual irriga- 
ticm target 
output (MAF) 

Capital cost 
in crores 

Cj^ 2 = 1-1912 “ °*°26l2 

Annual irriga- 
tion target 
output (MAF) 

CMl cost in 
crores 

^,2 ^ O. 2 OI 2 

Annual energy 
target outputs 
(106 iwhr) 

Gross benefit 
frcTO energy 
in crores 

EB 2 2 ^ 0,03075E2 - 0.0000039E 

Annual energy 
target output 
ilO& Jcwhr) 

Capital cost 
in crores 


Annual energy 
target output 
(lOo kwhr) 

CMR cost in 
crores 

02^2 “ O.OOO 39 E 2 


7 

Note : 1 crore = 10 . 
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TABLE 4.4,3. COST AND BENEFIT FUNCTIONS 


OMKARESHWAR 


SI, 1 

Independent 

t 

1 

I Dependent 

Fxinctional relationship 

No. ! 
1 

1 

variable 

; variable 

I 

L.. — — J 

L- - - - - 


1 

Reservoir 

Capital cost 


capacity (MAF) 

in crores 

2 

Reservoir 

OMR cost in 


capacity (MAF) 

crores 

3 

Annual irriga- 

Gross benefit 


tion target 

from irriga- 


output (MAF) 

tion in 
crores 

4 

Annual irriga- 

Capital cost 


tion target 
output (MAF) 

in crores 

5 

Annual irriga- 

CMR cost in 


tion target 
output (MAF) 

crores 

6 

Annual energy 

Gross benefit 


target output 

from energy 


( 10 ° Jcwhr) 

in crores 

7 

Annual energy 

Capital cost 


target output 
( 10 ® kwhr) 

in crores 

8 

Annual energy 

OMR cost in 


target output 
( 10 ^ Kwhr) 

crores 


7 

Notes 1 crore= 10 . 


C3 2 = 2.O83Y2 - O.SSSY^ 

© 33 = 0.000 

IB^ 3 = 2,73513 - O.O8I4I3 

© 13 = O.955I3 - O.O3O4I3 
°1,3 = 

EB„ - = 0.021 15E- - 0,0000047E 

^ i ^ O 

C^ - - 0.00099E„ 

°2 3 “ O-OOO 4 OE 3 


cxj CO 



TABLE 4.4.4. 


COST AND BENEFIT 



3 Annual energy Gross benefit 

target output from energy 

( 10 ° Icwhr) in crores 


4 Annual energy Capital cost 

target oxitput in crores 

( 10 ° Icwhr) 


5 Annual energy OMR cost in 

target output crores 

(10° Icwhr) 


Note: 1 crore = 



100 


Functional relationship 

I 

S,4 = 

4 = 0.000 

EB 2 4 = 0.02437E^ - 0.0000089e| 
S 4 ^ 0-0014 55E^ 

0.00059E. 

J2,4 4 
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TABLE 4.4,5. COST AND BENEFIT FUNCTIONS 


SARDAR SAROVAR 


T 


SI. 

No. 

; Independent 
! variable 

I Dependent 

I variable 

1 

- 1 - - - - 

1 

Reservoir 

Capital cost 


capac it y ( MAF ) 

in crore s 

2 

Reservoir 

OMR cost in 


capacity (MAF) 

crore s 

3 

Annual irriga- 

Gross benefit 


tion target 

from irriga- 


output (MAF) 

ticxi in 
crores 

4 

Annual irriga- 

Capital cost 


tion target 
OXJtput (MAF) 

in crores 

5 

Annual irriga- 

OMR cost in 


tica:i target 
output (MAF) 

crores 

6 

Annual energy 

Gross benefit 


target output 
( 10 ° ]<whr) 

from eneirgy 


in crores 

7 

Annual energy 

Capital cost 


taraet output 
( 10 ° kwhr) 

in crores 

8 

Annual energy 

CMR cost in 


target output 
( 10 ® Iwhr) 

crores 


J Fxonctional relaticnship 

C_ c = 1.074Y_ - 0.04Yc 

0 = 0,000 

IBi 5 = ^*9641^ - 0.262815 

5 = I. 22 OI 5 - O.OSOlg 

EB 2 5 = 0.0193E^ - O.OOOOOaiEg 
C 2 5 * O.OOIOSE^ 

0_ c = 0.00026E_ 


Note : 1 crore = 10 ‘ 
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= Annual OMR cost of irrigation work at site i; 

ENB„ . = Net annioal energy benefit at site i and is equal to 
'Z / X 

*^^ 2 ^ - = 2,1 - ° 2,1 

EB,, . = Gross annual energy benefit at site i^ 

C,. . => Annual capital cost of energy at site i; 

£.iX 

^2 L ~ Annual OMR cost of energy generation at site i; 

C,j . = Annual capital cost of reservoir at site i; 

Yj^ => Gross storage capacity of reservoir at site i. 

Sx±)stituting for the variables in the objective function/ 
the following expression is obtair^: 

B = 4.48 - 0.364 + 0.03288 E^ - 0.000029 E^ 

- (1.615 - 0.07/3 Y^) + 2.51 - 0.12 + 0.0288 

- 0.0000039 e| - (0,48& Y^ - 0.0138 Y^) + 1.43 I 3 

- 0.051 I 3 + 0.01976. E^ - 0.0000047 E^ - (2.083 Y^ 

- 0,555 Y 3 ) + 0.02232 E^ - 0.0000089 eJ - (4.528 Y^ 

« 5.16; yJ) + 4.445 I 3 - 0.2128 I 5 + 0.018 E^ 

- 0.0000021 Eg - (1.074 Yg - 0.040 Yg) (4.3) 

The optimization of the objective function is subject 
to the following constraints; 



(a) The volxoitie of water released from reservoir must be suffi- 
cient to meet irrigation demands of that period from that 
reservoir ^ i . e . , 


i^t 


K. ^ I, 
1/t i 


V. 

lyt 


(4.4 ) 


where a. . is the water released frcxn the reservoir i in 

lyt 

time ty and K. is the proportion of irrigation demand I . 

jL / w jL 

to be diverted for irrigation in time ty from site i; 

(b) The volume of water released during any period cannot 

exceed the contents of the reservoir at the beginning of 
the period plus flow into the reservoir during the period y 
i.e. y 


a, , "4*c« . ^ STj . ■{"f* , “^a. - , — K- ^ , X > 

i i jfi *!jL 0 "t# JL "““JL f 'ti jL JL ^ TI> «L 


^i-lyt ^i ^ *^i-lyt 


iyt 


(4.5) 


where c. is the volume of water released thro\jgh river- 


bed turbine in time t at reservoir site i; 

t)S 

(c) The contents of the reservoir at, any period cannot exceed 

A 

the quantity left over from the previous period y i.e.; 


ST. . < ST. . + f . ^ + a. - -1 y. - K. ^ I. 

iyt-+-l lyt lyt i — 1 y t i -1 y t i— lyt i 


K! - ^ I . - a, ^ - G. ^ 

1-1 yt 1 iyt lyt 


iyt 


(4,6) 


(d) The contents of reservoir at any period cannot exceed the 
capacity of reservoir 


STy . < Y. 

iyt — i 


iyt 


(4-7) 


where ST^ ^ is storage in reservoir i at the beginning of period t; 



(e) The flow through the turbines should meet energy generation 
demand. Variation in generation efficiency at different 
sites is neglected at this stage; 

<0 V. (for canal 

*“ power house) 

<0 V (for river bed 

' power house ) 

(4,8) 

6 

where 1.025 is the conversion factor of MAP to 10 }<whr, 
e is the turbine and generator efficiency, and, h^, is 
the average head on tiurbine at site i. 

A few more constraints were added on the basis of some 
design criteria as listed below: 

(f) It has been reccxnmended by the Tribxanal that the total 
water used by Madhya Pradesh should be atleast 1.947 times 

the amount used by Gujarat, i.e., 

4 

E I, > 1.947 Ic (4.9) 

i=l ^ ^ 

(g) The total energy produced should be atleast equal to 
41(X) X 10^ kwhr, i.e, 

5 

S E, > 4100 (4,10) 

i=l ^ 


6.E, - (1.025) (e)(h. )(a, . ) 

t 1 1 X/t: 

6^E^ - (1.025)(e)(h^)(C^^^) 


(h) The dead storage of a reservoir puts a lower limit cxi the 
• reservoir storage in the case of reservoir producing power, 
i . e . , 


ST 


l,t 


> D, 


(4.11 ) 



where D. is the dead storage capacity of reservoir i- 
(i) The Full Supply Reservoir Level (FRL) is fixed at +45^ for 
the terminal reservoir (Sardar Sarovar) by the Tribunal. 

The gross storage at this level is 7.75 MAP, i.e., 

Y. = 7.75 (4.12) 

The return flow is taben as 10% of the annual water xised 
in the upstream. Evaporaticn losses and flood storages are 
not considered in the LPD model. Deterministic streamflows 
are used in this model. 

An illustration of application of the above constraints 
to, say, Sardar Sarovar is given by the following equations: 


- a^ ^ + 0.1934 < 0 

”^5 2 ^ 5^0 

- a. o + 0.2931 < 0 

- a^ ^ + 0,1645 < 0 
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^■^5,2 
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> 2. 
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S%,4 

> 2. 
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7,75 


(4.13) 
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4.6* Method of Solution 

The objective is to find the values of variables 

Y/ E for each reservoir ^ that maximize the objec- 
tive function siibject to the constraints mentioned above. The 
constraint equation (4.8) is nonlinear becaxise head is the 
function of storage which is variable quantity. Energy produc- 
tion is the function of the product of the head on the t\arbine 
and discharge passing through turbine. It is made linear by 
considering average head on turbine. 

The objective function is also nonlinear. However it 
is seen that this nonlinearity is of special kind namely it is 
sum of several nonlinear functions each of which is a function 
of one variable only. Such nonlinearities can be handled by the 
separable linear prograimning technique. In this technique each 
nonlinear function is replaced by a piece-wise linear function. 
The resulting estpression will be linear and can be solved by 
the familiar linear programming technique. 

It is well known that linear programming can be used 
generally to maximize concave functions or minimize convex 
ftinctions. If the nonlinear concave or convex functions are 
relati\rely smooth the global optimum solution resxilts. If not 
there is a possibility of getting multiple optima (Loucks^ 

1981, page 58). 



CHAPTER 5 


SIMULATION 


5,1 . Introduction 

Simulation attempts to stxjdy the performance of a 
system and examines conseqxoent economic effect of the system 
by a ccxnputerized algorithm. Simulation models are more popiular 
due to their simplicity and versatility. They are widely used 
in the discipline of science and technology. Simulation studies 
a complex system which cannot be analyzed by the vise of analy- 
tical methodologies. Simulation studies have been greatly 
facilitated by the advent of high speed ccanputers , The economic 
response of the various alternatives can be examined^ if eccxnomic 
data for each project aind other activity in the basin can be 
identified well in advance, Ncaillneary dynamic, and stochastic 
response of the system can be incorporated in the simulation 
program. 

A major difficulty with simulation technique is, to 
develop a simulation program. Debugging and testing of the 
program are the other difficulties associatacl with it, A large 
amount of input and output data to be handled during the execu- 
tion is also a difficulty^ one has to face. 

A slmvilation model is time sequenced. Weekly^ ten- 
daily,, or monthly slmvilation stucjy is adopted depending on the 

f ■ 

type of problem and information required from it. For flood 
aimulatlcn, the time is reduced to 6-hours or 3-hours. Many 
assxjmptions are made while developing a simulation program. 



For larger time interval some of the approximations are invalid 
and thxis the simulation results become unreliable. 

There are three main limitations to simulaticai appro- 
ach. It does not yield a xinigiie optimal solution. It gives 
a large set of alternative design variables. It is a tedious 
job to examine such large sets. To overccxne these difficulties 
partly^ sampling methods for efficient selection of combinatlcai 
of variables are used. A second limitation of simulation 
results from the lack of flexibility in operating procedure of 
the system. This can be overcame partly^ by using flexible 
operating rules. The third limitation of simulation derives 
from the use of historical streamflow record. This can be 
overcome by using synthetic streamflows. 

In the present study, the sampling method used to 
explore the response surface is the steepest ascent method. 

The flexibility Introduced in the operating policy is in the 
form of modification in the release policy according to the 
hydrological classif icatian of the basin as very wet, wet, 
average, and, dry year. This is discussed later in Section 5.4 

As reported in Chapter 4 the system under study 
consists of five reservoirs. Out of these, four reservoirs 
(Bargl, Narmadasagar, Qmkareshwar, Sardar Sarovar) have three 
design variables namely irrigation quantity, reservoir capacity, 
and energy. The remaining one (Maheshwar) has two, namely 
reservoir capacity and energy. Thus there are 14 variables. 

The net benefits to simulating combinations of these 14 variables 
in the Narmada basin under stxjdy, creates a 14-dimensional 
response surface, that includes atleast cxie point with the 



highest net benefits. Even if we assume a coarse grid of 
5 divisions for each variable it would be required to examine 
5 6 ^ combinations. The only practical solution to 

find the optimal combination out of these 14 variables is to 
use a suitable method of sampling the variables to explore the 
net benefit response surface. 

5.2. Sampling Methods 
5.2.1. Types of methods 

There are five methods for exploring the response 
surface (Maass^ et al., 1962): 

(a) Uniform grid method 

(b) Randcan sampling 

(c) Single factor 

(d) Marginal analysis 

(e) Steepest ascent method. 

(a) Uniform grid method: In this method each of the n vari- 
ables is divided into m equal increments. Net benefits are 
simulated for each of m^ points on the resulting uniform grid. 

The point on the surface with the highest net benefit is selected. 
The advantage of this method is its ability to map the entire 
response surface of the system with small nuirher of variables. 

The main disadvantage of this method is that it spends more time 
in examining areas of low net benefit compared to search the 
higher areas for the peak. 

(b) Random sampling: Random values of variables are picked 

up from its population « and net benefits for the resulting point 



on the response sxarface are selected. This method is capable 
of handling large number of decision variables. 

(c) Single factor method: Initial values are selected for 
the variables. The values of all variables except one, say the 
first one, aire held constant, and the value of the first vari- 
able is altered \jntil the maximum net benefit is attained. 

This process is repeated with a second variable and so cffi. 

This method works best when the variables are independent of 
one another. The decislcMi variables are interrelated particu- 
larly for water resource system, hence this method is inherently 

ctfote, 

unsuit'^ to water resource systems, 

(d) Sampling by incremental or marginal analysis: In this 
method two variables are altered at a time, while remaining 
variables are held constant. The relationship between the 
variables is isolated. It provides the information regarding 
the direction in which the change should be made. This method 
is most suitable when there are few variables. The important 
pairs are identified by prior analysis. 

(e) Steepest ascent method: This is an Iterative method. 
Sampling moves sequentially from lower to higher elevation on 
the response surface. The starting points are selected on 
the basis of information available from preliminary design. 

Each variable is incremented/ keeping other variables constant, 
by preselected amount Ax^^. The change in benefit value AB^^ 
is calculated for each variable. The associated change in net 
benefits and its rate of change, ccmputed, 

A move is calculated from the starting values to the first 



revised base valxaes It requires that the change in each 

variable be made propoartional to the associated rate of impro- 
vement A B^/ A . The move to the new base is syiribolized by 
the equation; 

X^ - X° = AB^/AX^) (5.1) 

where is the proportionality constant. The relationship 
between the distance ’d* from the original base to the revised 
base is symbolized by the equation; 

% = S ( AB_^/A (5.2) 

Alteimatively can be determined as below (Maass, et al., 
1962). Cg can be of such magnitude that the maximtam absolute 
value of change for any variable will equal a specified magni- 
tude Kg or max ~ ^ 5 ' fran equation (5.2) 

that 

Cg max lABj^/AX^l = Kg 

Kg/max IAB^/AX^I < Kg/ IAb^/ A x^ 1 

^ ^5 

Cg = Kg min | AX^^/ AB^, 1 (5.3) 

The routine of calculations is repeated \intil there is 
no fxirther improvement. To xise this method/ the variables must 
have same units. This is not possible in water resoxirce system* 
To overcome this difficxxlty/ the \inits of the system eompcnents 
as well as sampling parameters d and I^are made dimensicnless by 
expressing them as ratios of their permissible range. 



The major advantages of steepest ascent method is that 
it measxxres cc«nbined influence of variables on the net benefits 
of the system in each lift- The major disadvantage is that it 
may miss trxie surnnit on a response surface- 

5,2.2. Application of steepest ascent method to the system 

under stxjdy 

A computer program is developed for the steepest ascent 
method to suit the reqxoirements of the present system. The 
range of variables in natural units are selected on the basis 
of the information obtained frcan the results of LPD model. 

The size of the step for each variable is set at about cane 
fifth of its range- The objective fxonction developed for the 
LPD model is •used to represent response s\irface. A set of 
design variables is selected arbitrarily as starting base 
values. The net benefit for this starting set is calculated. 
Each variable is then incremented^ keeping other variables 
constant. The change in the benefit val\:ie AB^ dxje to incre- 
menting each variable is calculated. The step valxie A in 
natural unit is converted into range \init for each variable. 

The ratio, AX^/A is calculated for each variable. An 
absolute minimtjm valxae out of these ratios is found. The 
value of Cg is calculated as the product of absolute minimiam 
value X20. The move from original base to new base is equal to 
AX^/ ABj^), which is then converted to natural units fear 
each variable. A second base = original base + move. The 
revised net benefit is calculated. It is then canpared with the 
net benefit already calculated for the starting base. If the 
valtje is positive then the procedtire is repeated, considering 



the second base as the original base. The procedure is repea- 
ted until no improvement is found in the benefit value. It is 
fovind that it takes about 7 iterations for this phase. A 
second sample is taken and similar computations made. Results 
of 9 such samples are given in the next chapter. Tables 6.3.1 
to 6.3.9. Based can these tabular values the final selected 
set is shown in Table 6.3.10. 

5.3. Data for Simxilation 

Nairmada river basin as modeled has five reservoirs: 

(i) Bargi (ii) Narmadasagar , (iii) Omkareshwar. (iv) Maheshwar 
and (v) Sardar Sarovar. Pour serve two purposes namely irri- 
gation and power. One reservoir namely Maheshwar is a hydro- 
power project. Canal head power hoToses (CHPH) are provided at 
reservoir sites (i), and (v) . River bed power houses (RBPH) are 
provided at the remaining reservoir sites. 

In the present study two aspects are considered in the 
simulation study: (a) monthly simulation, (b) 3-hourly flood 
simulation. Monthly simulation study is carried out for conser- 
vation ptirposes, i.e,, reservoirs are filled up during monsoon 
period (filling period) and stored water is used to meet the 
demand d’urii^ nonmonsoon period (depletion period). Monthly 
discharge data at each reservoir site are given in Tables 
3.4.10 to 3. 4. 14. 3-hourly flow data for the four selected 
floods for selected reservoir sites are given in Tables 3.4.18 
to 3.4.25. The types of flood classification on the basis of 
monthly mean flow and trigger value are given in Tables 5.3.1 
to 5.3.4. Ihe release policy diiring flood month based cn 



TABLE 5.3,1. TYPE OF FLOODS AND RANGE FOR SELECTION 
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TABLE 5.3,2. TYPE OF FLOCKS AND RANGE PCR SELECTION 
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TABLE 5.3.3. TYPE Cff* FLOODS AND RANGE FOR SELECTION 


r-| CO E 
fO CO Q) < 

POPS 

+J O H 
0X0 4-) 

< o > (d 



>rt ^ 
O X CO 

O cd CM 


o 

+j cd 
Cd CM o 
Q O ft 



3-ho\irly mean flews is given in T^le 5.3.5. Other data liJce 
elevati on-area-capacity ^ monthly evaporaticn rate^, monthly 
release policy, and maximum and minimxan permissible storage 
at the end of each month for each reservoir, are 
given in Tables 5.3.6 to 5.3.19. 

The net inflow to any reservoir is taken as the inflow 
due to indeiendent catchment area - (u/s water requirements for 
major projects, medium, minor and pumping schemes) + Regenerated 
flews + Release frean the u/s reservoir. 

5.4. Simulation Procedure 

A simulation program developed by the HBC (1966) 
namely 'Reservoir yield' is used in the present study after 
making suitable modification in the program. The design vari- 
ables obtained frean the steepest ascent method are used as 
input data. Other data needed for simulation are described in 
Section 5.3. Ihe system performance is tested using historical 
streamflcw data of 28 years and generated streamflcw data of 
50 years for each simulation rtua. 

The water requirements and releases are taken from 
the NWDT award (see Appendix). A sfudy is also carried out 
\asing modified water requirements and releases modified accor- 
ding to the classification of the basin as very wet, wet, 
average and dry year. The 75% dependable flew at Sardar 
Sarovar denoted by, say, w is fixed as 27 MAP cn the basis of 
79 years of annxaal flow series available at Sardar Sarovar. 

If the flow at Sardar Sarovar in any year is more than 1.5 w, 
the year is termed as very wet, and the water raqiiirements and 
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TABIE 5.3.5, FLOOD RELEASE POLICY FOR FLOCD MOITH 


Total 3-ho\arly mean flow 
= Flood ordinate 

+ Release from upstream 
reservoir 

(cfs) 

- - - - - 1 

Flocd release 
(cfs) 

0 - 50000 

15000 

50001 - 75000 

20000 

75001 - 100000 

25000 

100001 - 200000 

30000 

200001 ~ 300000 

35000 

300001 - 400000 

40000 

400001 - 500000 

45000 

500001 - 600000 

50000 

600001 - 750000 

60000 


y 750001 


70000 





X4X 


TABLE 5.3.6. ELEVATION -AREA-CAPACITY DATA 

BARGI 


• 

1 

SI. i 

No. 1 

1 

1 

1 

. - - - i - 

Elevaticffi in 
feet 

1 

I 

1 

1 

f 

f 

1 

1 

1 

t 

1- 

Area in 
acres 

Storage capacity in 
MAP 

L. - - - - 

1 

1300 


6225 

0.210 

2 

1320 


14993 

0.548 

3 

1324 


18278 

0.600 

4 

1340 


25565 

1.100 

5 

1360 


41250 

1.800 

6 

1380 


60000 

2.749 

7 

1387 


66750 

3.180 

8 

1390 


69750 

3,406 

9 

1400 


83500 

4.055 

10 

1410 


97500 

4.750 

11 

1420 


115000 

5.650 

12 

1426 


130200 

6.250 

13 

1430 


143500 

6.600 

14 

1435 


160000 

7.500 





TABLE 5.3,7 


ELEVATIC»J-AREA-CAPACITy DATA 


NARMADASAGAR 


f 

SI. i 

No, ; 

• 

• 

Elevation in J 

feet J 

, - - 

Area in 
acres 

1 

1 

; Storage capacity in 

1 MAF 

t 

. 1 - — 

1 

690 

000 

0,000 

2 

700 

2700 

0.120 

3 

710 

3900 

0.160 

4 

720 

4240 

0.193 

5 

730 

5000 

0.210 

6 

740 

6960 

0.304 

7 

750 

11000 

0,385 

8 

760 

18040 

0.505 

9 

770 

27000 

0.770 

10 

780 

38040 

1,100 

11 

790 

52000 

1.500 

12 

798 

63800 

2.000 

13 

800 

66040 

2.150 

14 

810 

82000 

2.870 

15 

820 

99800 

3,760 

16 

830 

121000 

4.700 

17 

840 

147300 

6.210 

18 

850 

180000. 

7.800 

19 

860 

224840 

9.900 

20 

864 

238000 

10.800 

21 

870 

266000 

12.500 

22 

880 

308000 

16.000 

23 

890 

360000 

20.000 

24 

900 ■ 

430000 

25.500 







ELEVATION-AREA-CAPACITY DATA 


TABLE 5.3.8. 


OMKARESHWAR 


1 

SI. i 
No. ; 

- ,.-i- 

t 

1 

Elevation in | 

feet I 

I 

I 

Area in 
acres 

f 

I Storage capacity in 
; MAP 

i- - - - - - 

1 

510 

0 

0.0000 

2 

520 

75 

0.0025 

3 

530 

300 

0.0075 

4 

540 

650 

0.0130 

5 

550 

1126. 

0.0225 

6 

560 

1700 

0.0410 

7 

570 

2400 

0.0630 

8 

580 

3300 

0.0910 

9 

590 

4540 

0.1280 

10 

600 

6100 

0.1800 

11 

610 

8050 

0.2470 

12 

620 

11100 

0.3470 

13 

630 

15100 

0.4800 

14 

635 

17950 

0.5600 

15 

640 

21100 

0 , 6560 

16 

650 

28360 

0.9000 

17 

660 

34976 

1.2160 

18 

665 

41500 

1.4375 

19 

670 

45952 

1.6190 

20 

675 

55000 

1.8750 

21 

680 

75000 

2.5000 






TABLE 5.3.9. ELEVATION-AREA-CAPACITY DATA 

MAKES HWAR 


SI. 

No. 

t 

1 

1 

t 

Elevation in 
feet 

j 

1 

t 

1 

1 , 

Area in 
acres 

f 

; storage capacity in 
; MAF 

L.-. - - - - - : - 

1 


456 


0 

0.00000 

2 


460 


500 

0,00067 

3 


470 


1300 

0.00935 

4 


480 


2400 

0.02757 

5 


490 


3700 

0.05741 

6 


500 


5400 

0.10307 

7 


510 


7000 

0.16490 

8 


520 


9350 

0.24537 

9 


530 


11000 

0.34800 

10 


532 


11800 

0.37500 

11 


534 


12000 

0.39800 

12 


540 


13800 

0.47175 

13 


545 


16000 

0.55000 

14 


550 


18150 

0.63100 

15 


555 


21000 

0.75000 

16 


560 


24000 

0.85000 

17 


565 


28000 

0.95000 

18 


570 


34000 

1.15000 

19 


575 


38000 

1.30000 

20 


580 


44000 

1.50000 

21 


585 


52000 

1.70000 

22 


590 


64000 

2.00000 

23 


595 


80000 

2.30000 
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TABLE 5.3,10. EIEVATION-AREA-GAPACITY DATA 

SARDAR SAROVAR 


T r 


SI. ; 
No. ! 

I 

Elevaticffi in 
feet 

1 

i 

1 

-- A 

Area in 
acres 

J Storage capacity 
; MAP 

t 

1 

100 


1208 

0.000 

2 

120 


1998 

0.080 

3 

140 


2761 

0.100 

4 

160 


4155 

0.150 

5 

180 


5162 

0.240 

6 

200 


6468 

0.360 

7 

220 


8260 

0.500 

8 

240 


9968 

0*700 

9 

260 


12340 

0.920 

10 

280 


14340 

1.120 

11 

300 


17688 

1.540 

12 

320 


21188 

1.940 

13 

340 


22852 

2.380 

14 

360 


27224 

2.880 

15 

363 


28200 

2.970 

16 

380 


33372 

3,480 

17 

400 


41532 

4,280 

18 

420 


55148 

5.190 

19 

440 


70340 

6.440 

20 

455 


91340 

7.700 

21 

460 


99969 

8.140 

22 

465 


107500 

9,000 

23 

470 


117500 

9,500 

24 

475 


125000 

10.150 

25 

480 


135000 

11,000 

26 

492 


150000 

13,230 
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TABLE 5.3.11. MOtWHLY EVAPORATION RATE 

BARGI NARMADASAGAR 


SI. 

No. 

f 

I Month 

f 

t 

t 

1 

1 1 

t f 

1 Evaporation rate J 

I in inches I 

1 t 

1 - - - . 1 

SI. 

No. 

- - 

I Month 

J - 

I Evaporaticai rate 
; in inches 

1 

July 

2.5 

1 

July 

5.44 

2 

Aug 

2.3 

2 

Aug 

5.48 

3 

Sept 

3.1 

3 

Sept 

5.10 

4 

Oct 

4.0 

4 

Oct 

5.32 

5 

Nov 

3.1 

5 

Nov 

4.20 

6 

Dec 

2.4 

6 

Dec 

3.60 

7 

Jan 

2.8 

7 

Jan 

4.30 

8 

Feb 

3.7 

8 

Feb 

5.10 

9 

March 

7.3 

9 

March 

8.76 

10 

April 

11.6 

10 

April 

13.20 

11 

May 

15.6 

11 

May 

15,20 

12 

Jxine 

7.7 

12 

June 

10.00 
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TABLE 5.3.12. MONTHLY EVAPORATIOJ RATE 


OMKARESHWAR MAHESHNAR 


SI. 

No, 

1 

1 

I Health 

1 

( 

1 

1 - 

t I 

; Evaporation rate J 

; in inches ; 

1 • 

..- 1 - - - - . ■ • 

SI. 

No. 

1 

I Month 

I 

1 

t 

1 

t 

; Evaporaticai rate 

I in inches 

1 

1 

July 

5.1 

1 

July 

5.1 

2 

Aug 

4.8 

2 

Axjg 

4,8 

3 

Sept 

5.0 

3 

Sept 

5.0 

4 

Oct 

5.2 

4 

Oct 

5.2 

5 

Nov 

3.8 

5 

Nov 

3.8 

6 

Dec 

3.3 

6 

Dec 

3.3 

7 

Jan 

3.8 

7 

Jan 

3.8 

8 

Feb 

4.5 

8 

Feb 

4.5 

9 

March 

6.6 

9 

March 

6.6 

10 

April 

8.1 

10 

April 

8.1 

11 

May 

10.6 

11 

May 

10.6 

12 

June 

8.1 

12 

June 

8.1 
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TABLE 5.3.13. MC5NTHLY EVAPORATION RATE 


SARDAR SAROVAR 


1 — 

1 

SI. No* ; 

• 

f 

... - I 

Month 

1 ■' • 

t 

I Evaporation rate in inches 

1 

f 

1 

July 

3.30 

2 

Aug 

3.00 

3 

Sept 

4.44 

4 

Oct 

6.46 

5 

Nov 

5.87 

6 

Dec 

4.73 

7 

Jan 

5.50 

8 

Feb 

6,39 

9 

March 

9.80 

10 

April 

12.00 

11 

May 

12.90 

12 

Jme 

8.00 
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TABLE 5.3.14, MOSITHLY RELEASE POLICY 

BARGI NARMADASAGAR 


SI. 

No. 

T ! — 

{ Month { 

1 t 

1 1 

1 t 

Release in 
cfs 

1 1 

1 1 

1 ! SI . 

I 1 No. 

t 1 

i-J - 

t t 

1 Month 1 

f f 

1 t 

i 1 

Release in 
cfs 

1 

July 

260 

1 

July 

13224 

2 

Aug 

3306 

2 

Aug 

24795 

3 

Sept 

3306 

3 

Sept 

56203 

4 

Oct 

104 

4 

Oct 

3306 

5 

Nov 

104 

5 

Nov 

3306 

6 

Dec 

52 

6 

Dec 

3306 

7 

Jan 

52 

7 

Jan 

3306 

8 

Feb 

52 

8 

Feb 

3306 

9 

March 

52 

9 

March 

1653 

10 

April 

52 

10 

April 

1653 

11 

May 

52 

11 

May 

1653 

12 

June 

52 

12 

June 

1653 


Total 


7444 

(0,45 MAP) 


Total 


117364 
(7.10 MAP) 
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TABLE 5.3.15. MONTHLY RELEASE POLICY 

OMKARESHWAR MAHESHWAR 


SI. 

No. 

t f 

J Month J 

1 1 

I t 

1 f 

» I 

Release in 
cfs 

1 1 

1 1 

1 ! SI . 

1 ! No. 

1 1 

X 1 

“T r 

I Month ; 

t 1 

r f 

1 1 

. .. i - - . - -X 

Release in 
cfs 

1 

July 

13224 

1 

July 

11569 

2 

Avig 

24795 

2 

Aug 

24795 

3 

Sept 

56202 

3 

sept 

59508 

4 

Oct 

4133 

4 

Oct 

4133 

5 

Nov 

4133 

5 

Nov 

4133 

6 

Dec 

4133 

6 

Dec 

4133 

7 

Jan 

4133 

7 

Jan 

4133 

8 

Feb 

4133 

8 

Feb 

4133 

9 

March 

2479 

9 

March 

2479 

10 

April 

2479 

10 

April 

2479 

11 

May 

2479 

11 

May 

2479 

12 

June 

2479 

12 

June 

2479 


Total 

124803 
(7,55 HIAF) 


Total 

126453 
(7.65 MAF) 
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TABLE 5.3.16. MONTHLY RELEASE POLICY 


SARDAR SAROVAR 


1 

SI, No. i 

- ■ .1 - 

Month 

I Release in c£s 

1 - ...... - - 




1 

July 

16 

2 

Aug 

99 

3 

Sept 

1& 

4 

Oct 

16 

5 

Nov 

16 

6 

Dec 

16 

7 

Jan 

16 

8 

Feb 

16 

9 

March 

16 

10 

April 

16 

11 

May 

16 

12 

June 

16 


Total 

275 

(0.0ia7 MAF) 
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TABLE 5.3.17. MINIMUM AND MAXIMUM PERMISSIBLE STORAGE AT 

THE END OF EACH MONTH 

BARGI NARMADASAGAR 


SI. 

No. 

Monthi 

Minimum 

MAP 

Maximum 

MAF 

f I 
f i 

; ! SI. 

! 1 No. 

t 

I 

; Month 

t 

1 

Minimum 

MAP 

; Maximum 
; MAP 


\ 

- - 


t 1 

II 

1 

.1 



1 

July 


1.8750 

1 

July 


4.5000 

2 

Aug 


4.3750 

2 

Aug 


10.5000 

3 

Sept 


6.2500 

3 

Sept 


15.0000 

4 

Oct 


5.6250 

4 

Oct 


13.5000 

5 

Nov 


4.6875 

5 

Nov 


11.2500 

6 

Dec 

0.6000 

4.0625 

6 

Dec 

2.0000 

9.7500 

7 

Jan 


3.7500 

7 

Jan 


9.0000 

8 

Feb 


3.4375 

8 

Feb 


8.2500 

9 

March 


3.1250 

9 

March 


7.5000 

10 

April 


2.8125 

10 

April 


6.7500 

11 

May 


2.5000 

11 

May 


6.0000 

12 

June 


2.1875 

12 

June 


5.2500 
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TABIE 5.3.18. MINIMUM AND MAXIMUM PERMlfSSIBIE STORAGE AT 

THE END OF EACH MOlSftH 

OMKARESHWAR MAHESHWAR 


SI. 

No. 

t 

1 

{ Mcaith 

1 

1 

1 

» 

Minimxm | 

MAP ; 

1 

. . -■ 

Maximviin 

MAF 

TT 

i i SI. 

I I No. 

-1 r 

; Month { 

f 1 

1 f 

1 

Minimum 

MAF 

Maximum 

MAP 

1 

July 


0.6400 

1 

Jxily 


0.3750 

2 

Aug 


1.1200 

2 

Aug 


0.3800 

3 

Sept 


1.6000 

3 

Sept 


0.4000 

4 

Oct 


1.4400 

4 

Oct 


0.3950 

5 

Nov 


1 . 2800 

5 

Nov 


0.3900 

6 

Dec 


1.2000 

6 

Dec 


0.3850 



0.5600 




0.3750 


7 

Jan 


1.0400 

7 

Jan 


0.3800 

8 

Feb 


0.9600 

8 

Feb 


0.3750 

9 

March 


0.8800 

9 

March 


0.3750 

10 

April 


0.8000 

10 

April 


0.3750 

11 

May 


0.7200 

11 

May 


0.3750 

12 

June 


0.6400 

12 

Jxme 


0.3750 
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TABLE 5.3.19. MINIMUM J\ND MAXIMUM STORAGE PERMISSIBLE 

STORAGE AT THE END OF EACH MOOTH 

SARDAR SAROVAR 


SI. No. 

Month 

Minimuin 

MAF 

1. ... - - - 

Maximum 

MAF 

, - - ' 

1 

-T - - . , - „ 

July 


3.1000 

2 

Aug 


5.4250 

3 

Sept 


7.7500 

4 

Oct 


^ 7.3625 

5 

Nov 


7.1688 

6 

Dec 

2.9700 

6.9750 

7 

Jan 


6.7813 

8 

Feb 


6.5875 

9 

March 


6.3937 

10 

April 


6.2000 

11 

May 


/ 

6,0000 

12 

June 


5.7800 






releases are increased by 5%. If the flew at Sardar Sarovar in 
any year is greater than 1.1 but less than 1.5 W, the year is 
defined as wet. and water requirements and releases are increased 
by 2.5%. If the flow at Sairdar Sarovar in any year is less than 
0.9 W, the year is classified as dry and water requirements and 
releases are decreased by 7.5%. If the flow at Sardar Sarovar 
in any year is between 0.9 W and 1.1 W, the year is termed as 
average year, and water requirements and releases are taken same 
as normal demands and releases respectively. 

It is necessary to use seme type of loss functions to 
evaluate the economic consequences of failure to neet target out- 
puts. This failiire is due to uncertainties caused by stochastic 
natxire of both supply and demand. The 'shortage index as defined 
by Beard (1970) is used in the simulation program and is discu- 
ssed in the next section. In addition to Beard's shortage index 
concept^ a second method based on classification of shortages 
into different types ^ is also used. The different types are shewn! 
in a tabular form below. 


Annual irrigation shortages 


Particulars 

t 

1 

1 

- 1 

Shortage 

ratio 

Type 

1 


0.00 - 

0.05 

Type 

2 


0.06 - 

0.10 

Type 

3 


0.11 - 

0.15 

Type 

4 


0.16 - 

0.20 

Type 

5 


0.21 - 

0.25 

Type 

6 


0.26 - 

0.30 

Type 

7 


0.31 - 

0.35 

Type 

8 


0.36 - 

0,40 

Type 

9 


0.41 - 

0.45 

Type 

10 


0.46 - 

0.50 

Type 

11 


>0.51 





Monthly power shortages for each season 


1 

Particxilar 1 

^ L 

Shortage 

ratio 

Type 

1 

0.00 - 

0.50 

Type 

2 

0.06 - 

0.10 

Type 

3 

0.11 - 

0.15 

Type 

4 

0.16; — 

0.20 

Type 

5 

0.21 - 

0.25 

Type 

6 

> 0.26 



Shortage = Reqiaired quantity - Actual quantity supplied 


Shortage ratio 


Shortage 

Required quantity 


5.5. Concept of Shortage Index 

Piering (1969) gave the following expression for loss 
function; 


T 

L = E L, 
t=l 


where^ 


P K a. 

( E b. . Sh/. ) N.-" 

Xjrw V 


in which ^ 


(5,4) 


(5.5) 


L a= total loss in T periods'; 

= loss in period t; 

p = number of stobperiod i in period t; 

SH^ t ~ shortage in subperiod i (positive or otherwise 
zero) in period t; 

® number of consecutive shortages in period t; 





b. ^ = coefficient (derived from the relationship 

between economic loss and water shortage); 

” coefficients. 

Beard (1970) '.h . defined shortage index as the sum of 
the square of annual shortage ratios over a 100-year period. 
This shortage ratio is the difference between quantity of water 
req-uired and supplied, divided by the qxoantity reqviired. He 
gives the following expression; 


SHI = 


T 

E 

t=l 


(SH^) X 


100 

T 


where , 



in which, 

= quantity of water required in period t; 
Xj. = quantity of water supplied in period t; 
T = nxmber of years considered. 


(5.6) 


(5.7) 


One can get the same expression as (5.6) by substituting = 0; 
p = 1 ; =2 and assuming b as constant in equation (5.5). 

There are certain advantages claimed by shortage index concept 
in the economic evaluation of water resource system as listed 
below; 

(a) Shortage index is a single representative value for all 
shortages occurring during the economic life of a project 
for a given sequence of supply and demand. 

(b) The magnitude and frequency of shortages iure'^pcedict^. 

If shortage index is 0.25^ it may mean 25 annxjal shortages 



of 10% each or 10 annual shortages of about 16% each/ or 
100 annual shortages of 5% each/ over a period of 100 
years. Losses that occur near the beginning of the 
project life have severe effect than those occur near 
the end of the economic life of project. Hence average 
economic effect is considered for each project. 


5.5.1. Evaluation of constant b 

One can easily find the loss fxanction if the valxjfi of 
b^^^ is kncwn in equation (5.5). Assviming an average square 
shortage ratio, and b as a constant, one can write the loss as; 

L = b(SH^) (5,8) 

in which , 

2 2 
(SH ) = average of squared shortage ratio, i.e., E SH./T 

t=l 

2 

b(SH ) = annual average losses as a ratio to gross anniial 
benefits . 

2 100 

We know, shortage index SHI = E (SH ) x --r- 

t=l ^ ^ 


. SHI = E (SH^/T) X 100 
t=l 


. • . SHI = 100 X SH^ 


b X SHI 
100 


(5.9) 

(5.10) 


For the whole system, the loss function for irrigation is 



(5.11) 


^RR 


5 


in which/ 


^RR “ loss for irrigation in whole system in 

7 

10 rupees; 

7 

IR = annual benefits in 10 rupees/MAF; 

IRR = actual irrigaticffi target met in MAP; 

SHIj^ = shortage index for irrigation , 

The same loss function is assumed for energy. 


^NG 


5 

E 

i=l 


"2/i 


SHI 


2/i 


PE. EACT. 


(5.12) 


in which, 

= annual loss for energy in whole system in 

7 

10 rupees; 

PE = energy price/kwhr; 

EACT = actual energy produced in kwhr; 

SHI 2 = shortage index for energy. 

Now the problem is to determine the value of and b 2 

for each project. If the actual data of annual losses and 

annual shortages were available, one can evaliiate the values of 

b^ and b 2 by fitting a quadratic function of the type y = b 2 ^(SHj^) 

^ ~*2 

for irrigation and y = b 2 (SH 2 ) for power. 

Irrigation is given first priority in the present study, 
hence irrigation shortages may be more serioxis than pcwer 
shortages. As no data are available to construct qxaadratic 
function, the value of b^ for irrigaticn is taken as. 10.0 for 
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each project. The power shortages are considered to be less 
effective^ the value of b 2 is assumed as 8.0 (Martinez^ 1971). 

5.6. Simul at i on Rxans 

The following simulation runs were conducted; 

(1) Monthly simulation runs are condiocted by considering 
water reqxiirements and releases according to (i) tribunal 
awards (ii) classification of year as very wet, wet, 
average and dry. 

(2) In order to know 3-hoxarly maximxam spill from the terminal 
reservoir, the flood subroutine is activated and 3 -hourly 
flood roxiting is carried out for the flood month having 
monthly meanflow greater than the specified trigger valiae 
at ireservoir sites. 

For the first study, foiar simulation rxans for each 
case using historical flow and synthetic flow are carried o\it 
while for latter study two sim\xlation runs are carried oxjt. 

The results of the simulation runs are tabulated in Tables 


6.4,1 to 6.4,60 



CHAPTER 6 


RESULTS AND DISCUSSIONS 

6.1, Streamf low Generation 

Twenty eight years of historical streamflow record at 
the three gauging sites are available. Using this, two hundred 
and fifty years of synthetic flows are generated at eacdi reser- 
voir site using HEC-4 program. This program occupies about 
50 K memoary on DEC-1090 computer system. It requiires about 
80 seconds to generate two hundared and fifty years of data. 
Tables 6.1.1 to 6.1,5 show monthly and annual means and 
standard deviations of flow at all the sites. It may be seen 
that the mean and standaard deviation of generated flows are 
close to historical stareamflcws at all stations, 

6.2. LPD Model 

The following configurations are studied with LPD 
model (Table 6.2,1). The model is run for different cases 
such as linear and nonlinear objective function, and different 
flow conditions like 1.25 times average flow, average flow, 
and 0,9 times average flow. Another configxiraticai in which 
release is allowed to pass through txirbine for power generation 
and then to irrigate the land is also studied. Fourteen alter- 
natives out of the possible twenty fo\ir are studied. Each run 
takes about two minutes CPU time. The restilts ai?e tabulated in 
Tables 6,2.2 to 6.2,5. Based cai these results, limits on 
variotis design variables are selected as shewn in Table 6.2.6 
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TABLE 6.1.1. COMPARISON OF STATISTICAL PARAMETERS AND OTHER 

INFORMATION FOR HISTORICAL AND GENERATED 
STREAMFLOWS 





— 

— 




— 

Station 

Month 


Flow 

Flow 

Flow 

Flow 

Standard 

No, and 



(sTjm) 

(max) 

(min) 

(mean) 

deviation 

r^ame of 
station 

- T j 




cfs 

. - 

cfs 

cfs 

cfs 

cfs 


Jxily 

H* 

692713 

73431 

1408 

24740 

16496 


G* 

6254277 

61575 

0 

25017 

14617 


Aug 

H 

1307568 

83722 

8047 

46699 

16733 


G 

11617173 

78317 

4227 

46469 

17654 


Sept 

H 

772670 

69946 

3310 

27595 

15360 


G 

6752786 

84036 

2978 

27011 

14724 


Oct 

H 

196507 

20211 

805 

7018 

5449 


G 

1369017 

40374 

572 

6676 

5392 


Nov 

H 

54675 

8715 

352 

1953 

1648 


G 

440976 

5990 

227 

1764 

1114 


Dec 

H 

27457 

4302 

294 

981 

774 


G 

242725 

8969 

231 

97.1 

816 

101 

Bargi 

Jan 

H 

18987 

2723 

201 

678 

490 


G 

180467 

7945 

142 

722 

680 


Feb 

H 

11956 

1113 

137"’ 

427 

243 


G 

113807 

2840 

75 

455 

304 


March 

H 

8357 

805 

62 

298 

208 


G 

81078 

3501 

20 

324 

302 


April 

H 

4991 

736 

16 

178 

163 


G 

48405 

1975 

13 

194 

211 


May 

H 

1621 

201 

15 

58 

46 


G 

13667 

269 

6 

55 

48 


June 

H 

48786 

15991 

16 

1742 

3173 


G 

339758 

15960 

0 

1369 

1929 

Anmjal 

H 

3146288 

231120 

35895 

112367 

43061, 

Statistics 

G 

27754136 

232076 

26178 

111017 

36708 


*H stands for historical and G for generated flows. 








TABIjE 6.1.2 


COMPARISON OF STATISTICAL PARAMETERS AND OTHER 
INFORMATION FOR HISTORICAL AND GENERATED 
STREAMFLOWS 


— r~~-" 

1 

Station ; Mooth 
No . and { 
name of { 

station J 

1 

- , ..1. 


Flow 

(sum) 

cfs 

L-- 

Flow 

(max) 

cfs 

Flow 

(min) 

cfs 

L...,. 

Flow 

(mean) 

cfs 

. - - -J 

Standard 

deviation 

cfs 

July 

H* 

1356572 

188591 

11314 

48449 

37082 

G* 

10057156 

145813 

5770 

40229 

24100 

Aug 

H 

2945551 

195060 

19119 

105198 

45035 

G 

23946667 

195623 

13194 

95787 

40338 

Sept 

H 

2808165 

260660 

19770 

100292 

63115 

G 

24406818 

661820 

17059 

97627 

77641 

Oct 

H 

655923 

70928 

3704 

23426 

17103 

G 

6227084 

145280 

3049 

24908 

20392 

Nov 

H 

219812 

19425 

2061 

7850 

4750 

G 

2106158 

35090 

1625 

8425 

5707 

102 

H 

125550 

8234 

1536 

4484 

1793 

G 

1125604 

12418 

1232 

4502 

1882 

Narm'itd’^ 

Jan 

H 

90414 

5877 

1242 

3229 

1231 

G 

811943 

8330 

1084 

3248 

1307 

Feb 

H 

70948 

4299 

969 

2534 

807 

G 

617566 

5299 

808 

2470 

845 

March 

H 

52243 

3223 

387 

1866 

657 

G 

472963 

3271 

386 

1892 

664 

April 

H 

40622 

3330 

599 

1451 

631 

G 

367468 

5280 

516 

1470 

646 

May 

H 

27873 

1700 

466 

995 

362 

G 

252009 

2839 

362 

1008 

363 

Jxme 

H 

109749 

18199 

221 

3920 

4591 

G 

1014309 

32029 

88 

4057 

5092 

Annual 

H 

8503422 

641561 

98195 

303694 

120933 

Statistics 

G 

71405745 

912508 

97889 

285623 

112076 


*H stands for historical and G for generated flows 








TABLE 6.1.3 


COMPARISON OF STATISTICAL PARAMETERS ANB OTHER 
INFORMATION FOR HISTORICAL AND GENERATED 
STREAMFLOWS 




■ 

station ; Month 
No. and I 
name of \ 
station | 

t 

- .. - - 


Flow 

(s\jm) 

cfs 

-.1 

Flow 

(max) 

cfs 

L- 

Flow 

(min) 

cfs 

J 

Flow 

(mean) 

cfs 

- - - ^ 

standard 

deviaticai 

cfs 

L.,_ 

July 

H* 

94924 

13447 

763 

3390 

2638 

G* 

698066 

9873 

371 

2792 

1681 

Aug 

H 

206823 

13803 

1344 

7387* 

3205 

G 

1677367 

13963 

922 

6709 

2856 

Sept 

H 

198911 

18470 

1408 

7104 

4486 

G 

1729820 

47294 

1186 

6919 

5534 

Oct 

H 

46400 

5046 

263 

1657 

1213 

G 

441323 

10314 

216 

1765 

1448 

Nov 

H 

15591 

1380 

146 

557 

339 

G 

149673 

2499 

115 

599 

406 

103 

Cinkar- 

H 

8916 

587 

109 

318 

127 

G 

79924 

879 

88 

320 

134 

Jan 

H 

6426 

417 

88 

230 

88 

G 

57692 

591 

78 

231 

93 

Feb 

H 

5048 

306. 

69 

180 

58 

G 

43938 

381 

58 

176 

60 

March 

H 

3722 

230 

26 

133 

47 

G 

34508 

250 

20 

138 

54 

April 

H 

2899 

238 

43 

104 

45 

G 

34206 

1318 

14 

137/ 

141 

May 

H 

1997 

122 

33 

71 

26 

G 

21099 

482 

13 

84 

58 

J\ine 

H 

7714 

1294 

13 

276 

322 

G 

62665 

1662 

07 

251 

276 

Annual 

H 

599371 

45536 

6933 

21406 

8583 

statistics 

G 

5030281 

65149 

6860 

20121 

7993 


*H stands for historical aid G for generated flows. 










TABLE 6.1.4 


CCMPARISON OF STATISTICAL PARAMETERS AND OTHER 
INFORMATION FOR HISTORICAL AND GENERATED 
STREAMFLOWS 


1 

ff 

Station 1 Month 
No. and I 
name of ! 

station ; 

1 

-..1 - - 


— 

Flow 

(sum) 

cfs 

L- -- 

Flow 

(max) 

cfs 

- 

Flow 

(min) 

cfs 

Flow 

(mean) 

cfs 

- - -J 

Standa 

deviat. 

cfs 

July 

H* 

76114 

13673 

0 

2718 

2541 

G* 

550816 

3752 

0 

2203 

1207 

Aug 

H 

195818 

14933 

2031 

6994 

3443 

G 

1689339 

20970 

1296 

6757 

3553 

Sept 

H 

243877 

21724 

1751 

8710 

6295 

G 

2165054 

55836 

857 

8660 

7877 

Oct 

H 

57012 

9173 

0 

2036 

2051 

G 

517532 

7590 

0 

2070 

1706 

Nov 

H 

16357 

1938 

0 

584 

421 

G 

130629 

1048 

1 

523 

320 

104 

Mahe- 

H 

8462 

969 

0 

302 

191 

G 

7046t8 

481 

0 

282 

153 

Jan 

H 

5875 

805 

0 

210 

149 

G 

52844 

336 

0 

211 

104 

Feb 

H 

4684 

663 

0 

167 

120 

G 

39884 

251 

0 

160 

76 

March 

H 

3106 

413 

0 

111 

76 

G 

28022 

183 

0 

112 

58 

April 

H 

2187 

280 

0 

78 

56 

G 

19742 

158 

0 

79 

48 

May 

H 

1532 

219 

0 

55 

41 

G 

11091 

77 

0 

' 44 

21 

June 

H 

13451 

4194 

0 

480 

818 

G 

133857 

6067 

0 

535 

867 

Annual 

H 

628475 

48399 

7009 

22446 

10501 

Statistics 

G 

5409278 

73638 

4993 

21637 

9859 


*H stands for historical and G for generated flows. 


TABLE 6.1.5, COMPARISON OF STATISTICAL PARAMETERS AND OTHER 

INFORMATION FOR HISTORICAL AND GENERATED 
STREAMFLOWS 


Station 

Mcffith 


Flow 

Flew 

Flow 

Flew 

Standard 

No, and 



(sxjin) 

(max) 

(min) 

(mean) 

deviatican 

nanie of 
station 


- -i 

cf s 

- - 

cf s 

cfs 

- J 

cfs 

- 

cfs 


July 

H* 

2297025 

281835 

25340 

82037 

55346 


G* 

19574078 

344957 

17484 

78296 

41791 


Aug 

H 

5110391 

335597 

40642 

182514 

69689 


G 

44154058 

351609 

26872 

176616 

69621 


Sept 

H 

4984461 

451335 

33270 

178016 

114594 


G 

43302035 

999898 

31137 

173208 

132985 


Oct 

H 

1177105 

151068 

6471 

42039 

34411 


G 

11628423 

381146 

6378 

46514 

44577/ 


Nov 

H 

355097 

28624 

2714 

12682 

6945 


G 

3329294 

53689 

2157 

13317 

8596 


Dec 

H 

189706 

14283 

1895 

6775 

2654 

105 

G 

1753853 

22746 

1967 

7015 

3091 

Sardar 

Sarovar 

Jan 

H 

133506 

12571 

1557 

4768 

2221 


G 

1301783 

18086 

1614 

5207 

2408 


Feb 

H 

102368 

9501 

1288 

3656 

1605 


G 

984602 

12665 

1105 

3938 

1742 


March 

H 

71493 

6117 

881 

2553 

1132 


G 

708915 

8647 

582 

2836 

1258 


April 

H 

51690 

4222 

563 

1846 

872 


G 

518609 

8730 

673 

2074 

1026 


May 

H 

33879 

3201 

397 

1210 

596 


G 

329730 

7704 

429 

1319 

663 


June 

H 

271858 

58114 

487 

9709 

13344 


G 

2509253 

98978 

652 

10037 

12630 

Annual 

H 

14778579 

1036638 

164069 

527806 

210408 

Statistics 

G 

130094630 

1431669 

182185 

520379 

194600 


*H stands for historical and G for generated flows. 
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TABLE 6.2.1. TYPES OF CONFIGURATIONS 


C onf igxar ati on 

Name of I Purposes 

number 

project ; ; EnPTOv 

I Irriaation 

1 I f 

■ 1 

; ; RBPH ; chph 

1 1 t 



Bargi 

Yes 

No 

No 


N armadas ageu: 

Yes 

Yes 

No 

1 

Orokareshwar 

Yes 

Yes 

No 


Maheshwar 

No 

Yes 

No 


Sardar Sarovar 

Yes 

No 

YWS 



Bargi 

Yes 

Yes 

No 


Narmadasagar 

Yes 

Yes 

No 

2 

Orrikareshwar 

Yes 

Yes 

No 


Maheshwar 

No 

Yes 

No 


Sardar Sarovar 

Yes 

No 

Yes 



Bargi 

Yes 

No 

Yes 


Narmadas agar 

Yes 

Yes 

No 

3 

Omkareshwar 

Yes 

Yes 

No 


. Maheshwar 

No 

Yes 

No 


Sardar Sarovar 

Yes 

No 

Yes 



Bargi 

Yes 

Yes 


Narmadasagar 

Yes 

Yes 

4* 

Orrikareshwar 

Yes 

Yes 


Maheshwar 

Yes 

Yes 


Sardar Sarovar 

Yes 

Yes 


*In configuration 4 release is allowed to pass through txarbines 
for pcwer generaticaa and then to irrigate the land. | 

[ 

I 

" 

RBPH - River bed power house i 

CHE^ - Canal head pcwer house. 


! 
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TABLE 6.2.2. RESULTS OF LPD MODEL FOR CONFIGURATION 1 


Variable 


Objective 

function 



Linear j 

^ ^ 1 

Ncailinear 

i. ^ 


1.25 

times 

average 

flow 

I 

1 Average 

; flow 

1 

1 

1 

I 

-1---'.- 

I 0.9 

I times 

I average 

1 flow 
• 

1.25 

times 

average 

flow 

1 

5 Average 
; flow 

f 

1 

1 

1 

t - - 

i 0.9 

I times 

1 average 

I flow 

1 - - 

11 

8.930 

7-146 

6,431 

4.500 

4.000 

3.500 

Y1 

6.070 

4.980 

4.540 

2.880 

2.710 

2.430 

El 

0,000 

0.000 

0.000 

0.000 

0.000 

0.000 

12 

12.720 

6.660 

7.120 

11.000 

10.350 

9.250 

Y2 

13.260 

11.600 

10.450 

14.720 

13.250 

12.260 

E2'’" 

1601.000 

1802.000 

1456.000 

2000.000 

1526.000 

1400.000 

13 

5.260 

6.270 

5.080 

4,750 

4.150 

3.750 

Y3 

1.020 

0.880 

0.830 

0.560 

1.000 

0.940 


787.000 

693.000 

584.000 

1000.000 

( 

800.000 

750.000 

Y4 

0.375 

y 0.375 

\ 

’ 550.000 

0.375 

0.375 

0,375 

0.375 

E4'*’ 

559.000 

464.000 

600,000 

600.000 

550.000 

15 

13.800 

10.300 

9.550 

10.400 

9.500 

8.500 

Y5 

6.850 

7.750 

5.580' 

5,640 

6.080 

4.950 

E5* 

1253.000 

1153,000 

996.000 

1200,000 

1150.000 

1000.000 

Net 







benefits 
in 10^ 

188.240 

171. ISCP" 156,370 

125,680 

113.750 

107.480 

Rupees 








I, Y are In MAF and E is In million Icwhr, 
* CHFH considered 


+ RBPH considered 



X** 5r 


TABLE 6.2.3. RESULTS OP LPD MODEL FOR CONFIGURATION 2 



Objective function 

Variable 

Nonlinear 


Average flow 

11 

3.950 

Yl 

5.000 

El"^ 

300.000 

12 

10.400 

Y2 

11.060 


16-74.000 

13 

4.150 

Y3 

0.650 


800.000 

Y4 

0.375 


600,000 

15 

9.500 

Y5 

6.350 

E5* 

1150.000 

— 

Net benefits in 10 Rupees 

149.660 


I^ Y are in MAP and E is in million kwhr. 

+ RBEH considered 
* CHPH considered 

Kjptej In this ccaifiguration caily average flow case is included. 







TABLE 6.2.4 


RESULTS OF LFD MODEL FOR CONFIGURATION 3 


Variable 


Objective 

function- 




Linear 



' Nonlinear 


1.25 

times 

average 

flow 

; Average 

I flC3W 

I . 

i 0.9 

! times 

I average 

I flow 

1 

1.25 

times 

average 

flow 

Average 

flew 

0.9 

times 

average 

flow 

11 

8.930 

7 . 146 

6.430 


3.800 


Y1 

6.070 

4.980 

4.540 

- 

5.220 

- 

El* 

173.000 

139.000 

125.000 

- 

150.000 

- 

12 

11.700 

5.470 

4.530 

- 

10.000 

- 

Y2 

13.430 

11.800 

10.880 


12.500 

- 

E2'*' 

1752.000 

1979.000 

1839,000 

- 

1609.000 

- 

13 

5.240 

6.140 

7.560 

- 

4.000 

- 

Y3 

1.010 

0.880 

1.240 

- 

0.760 

« 

E3‘*' 

874.000 

805.000 

564.000 

- 

758.000 

- 

Y4 

0.375 

0.375 

0.375 


0.375 

— 

E4'^ 

617,000 

624.000 

468.000 

- 

575,000 

- 

15 

13.270 

f 

9,620 

9.500 

- 

9,600 

- 

Y5 

6.700 

5.880 

5.610 

- 

6.150 

- 

E5* 

1383.000 

1253.000 

1004,000 

- 

1175.000 

- 

Net 

benefits 
in 107 
Rupees 

213.430 

181.130 

162.130 

- 

116,230 

- 


I, Y are ±n MAF and E Is In milliai kwhr. 

* CHPH considered 
+ RBPH ccxisidered 

Note: Only average flow case is studied in case of ncnlinear 
objective functicn. 


TABLE 6,2.5 


RESULTS OF LPD MODEL FOR CONFIGURATI CW 4 


Variable 

t 

1 

1 

L— 

Objective functicai 

f 

t 

t 


Linear 


t 

I 

« 

1 

1 

- - 1 - 

1.25 times 
average flow 

• 

• 

J Average flew 
• 

1 

f 

J- . - - . 

. ' ■ ■ 

1 

1 0,9 times 

; average flow 

1 

-.-1 - - - - 

11 


8.750 

7.150 

6.430 

Y1 


6.125 

4,980 

4.540 

El 


193.000 

139.000 

125.000 

12 


6.740 

12.340 

12.300 

Y2 


15.340 

12.200 

10.690 

E2 


3514.000 

2692.000 

2328.000 

13 


10.390 

3.030 

1.530 

Y3 


1,040 

1.870 

1.650 

E3 


1793.000 

840.000 

656 ,000 

Y4 


1.490 

0,550 

0.430 

E4 


727.000 

450.000 

391.000 

15 


13.270 

11.540 

10.390 

Y5 


7.750 

7.750 

5.860 

E5 


1383.000 

778.000 

700.000 

Net benefits 
107 Rupees 

in 

268.220 

207.400 

1 81.^940 


1/ Y are In MAP and E is In million kwhr. 


Notes In this configuration only linear objective fxjnction is 
considered. 
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TABLE 6.2.6. LIMITS ASSUMED FOR EACH VARIABLE 


Variable 

f 

f 

1 

f 

1 

t 

1 

- -1 

Limit® assumed 

11 


2 


5 

Y1 


3 

- 

7 

El 


50 

- 

150 

12 


5 

- 

13 

Y2 


10 

- 

16 

E2 


600 

- 

1600 

13 


2 

- 

5 

Y3 


0.5 

- 

2.0 

E3 


600 

- 

1000 

Y4 


0.3 

- 

0.7 

E4 


200 

- 

600 

15 


7 

- 

10 

Y5 


7 


8 

E5 


600 

— 

1100 


Y are in MAP and E is in million kwhr 





for use in steepest ascent analysis. It is found that for 
each configuration the benefits are larger in the case of line 2 ir 
objective fvinction than in the case of nonlinear objective 
functicn. 

6.3. Steepest Ascent Method 

The range for each variable is fixed on the basis of 
infornvation obtained from LPD results. A set of variables is 
assumed arbitrarily within this range as a starting base. The 
objective f\inction developed for LPD model is used to represent 
response surface. The net benefits for nine samples are 
obtained and tabulated in Tables 6.3.1 to 6.3.9 as already 
mentioned in the last chapter. A set of variable which produce* 
maximum net benefit is selected as input data for simulation 
study. This set is shewn in Table 6.3,10. 

6.4, Simulation 

The design variables obtained frcrni steepest ascent 
method of sampling are used in simulaticn analysis. Simulation 
can be conducted by adopting varioxxs operating policies. An 
operating policy is a statement which stipulates the amount of 
water to be kept in storage at the end of a period/ say a 
month/ and the amount of water to be released during the month. 
This decision depends upc«n the config\iration of reservoirs/ arxi 
the purpose for which the system has been designed. The purpose 
may be conflicting or complementary. For instance/ flood 
control requires that the reservoir should be kept as empty as 
possible to hold expected flood/ whereas water should be held 
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TABLE 6.3.1. RESULTS OF STEEPEST ASCENT METHOD FOR SAMPLE 1 


Variable 

— 1 1 

1 

J Assumed 
i limits of 

5 variable 

t 

L— 



Range 



Starting 
base value 

• 

1 

1 Final valxie 

1 

1 

I 

f 

1 

... t _. - - 

1 

• 

2 

3 

4 

! 5 . 

11 

2 

- 5 

3,0 

2.0 

3,600 

Y1 

3 

- 7 

4.0 

3.0 

1.500 

El 

50 

- 150 

100.0 

50.0 

27.000 

12 

5 

- 13 

8.0 

5.0 

8.600 

Y2 

10 

- 16. 

6.0 

10.0 

8,430 

E2 

600 

- 1600 

1000.0 

600.0 

2000.000 

13 

2 

- 5 

3.0 

2.0 

2,800 

Y3 

0,5 

- 2,0 

1.5 

0.5 

0,230 

E3 

600 

- 1000 

400.0 

600.0 

761.000 

Y4 

0.3 

- 0.7 

0.4 

0.4 

0.398 

E4 

200 

- 600 

400.0 

200.0 

403.000 

15 

7 

- 10 

3.0 

7.0 

7.840 

Y5 

7 

- 8 

VO 

7.0 

6.960 

E5 

600 

- 1100 

500.0 

600.0 

882.000 


Net benefit- R3. 113.93 x 10^ 


I, Y are in MAP and E is in million Icwhr 
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TABLE 6.3,2. RESULTS OF STEEPEST ASCENT METHOJ FOR SAMPUE 2 


Variable 

* I 

; Assumed I 

I limits of J 

J variable J 

1 1 

L 

Range 

Starting 
base value 

Final value 

i 

1 

1 - 

2 . i 

- 3 


w 

11 

2 

- 5 

3.0 

2.15 

3.730 

Y1 

3 

- 7 

4.0 

4.00 

2.680 

El 

50 

- 150 

100.0 

60.00 

37.200 

12 

5 

- 13 

8.0 

6.00 

8.890 

Y2 

10 

- 16 

6.0 

11.00 

10.510 

E2 

600 

- 1600 

1000.0 

625.00 

2025.000 

13 

2 

- 5 

3.0 

2.25 

3.040 

Y3 

0.5 

- 2.0 

1.5 

0.75 

0.530 

E3 

600 

- 1000 

400,0 

650.00 

807.000 

Y4 

0.3 

- 0.7 

0.4 

0.41 

0.409 

E4 

200 

- 600 

400.0 

250.00 

445.000 

15 

7 

- 10 

3.0 

7.50 

8,220 

Y5 

7 

- 8 

1.0 

7.25 

7.210 

E5 

600 

- 1100 

500.0 

650.00 

931.000 


Net benefit 

-Rs.115.05 

X 10*^ 



If Y ajfe in MAP E 1b in million l<wbri 











TABLE 6.3,3. 

RESULTS OP 

STEEPEST 

ASCENT METHOD 

FOR SAMPLE 3 

Variable 

t 

1 Assumed 
! limits of 

I variable 

Range 

I 

f 

; Staarting 
; base value 

r 

1 

t 

1 

; Final value 

1 

1 

t 

1 


-1 - 



f 

■f-..- 


1 

1 

-J- - . ... 

2 

3 

* 4 

f.. - - 

* 5 

11 

2 

- 5 

3.0 

2.30 

3.800 

Yl 

3 

-7 

4.0 

5.00 

3.860 

El 

50 

- 150 

100.0 

70.00 

47.400 

12 

5 

- 13 

8.0 

7.00 

9.190 

Y2 

10 

- 16 

6.0 

12.00 

11.590 

E2 

600 

- 1600 

1000,0 

650.00 

2050.000 

13 

2 

5 

3,0 

2.50 

3.280 

Y3 

0.5 

- 2,0 

1.5 

1.00 

0.840 

E3 

600 

- 1000 

400.0 

700.00 

853.000 

Y4 

0.3 

- 0.7 

0.4 

0.42 

0.418 

E4 

200 

- 600 

400.0 

300.00 

487.000 

IS 

7 

- 10 

3.0 

8 00 

8.600 

Y5 

7 

- 8 

1.0 

7.50 

7.460 

E5 

600 

- 1100 

500.0 

700.00 

980.000 


Net benefit* Rs.ll6. 

32 X lo"^ 



I# Y are in MAF and E is in millican kwhr. 









157 


TABLE 6.3.4* RESULTS OF STEEPEST ASCENT METHC3D FOR SAMPLE 4 


Variable 

I 

• 

; Assximed 
; limits of 

1 variable 

f 

J 

Range 

1 

t 

1 Starting 

1 base value 

1 

1 

1 

Final valxjfi 

i - 

1 

^ 

2 ‘ 

3 

: 4 _ 

: 

11 

2 

- 5 

3.0 

2.45 

3.910 

Y1 

3 

- 7 

4.0 

6.00 

5.040 

El 

50 

« 150 

100.0 

80.00 

57.610 

12 

5 

- 13 

o 

f 

CD 

8.00 

9.460 

Y2 

10 

- 16 

6.0 

13.00 

12.660 

E2 

600 

- 1600 

1000.0 

675.00 

2075.000 

13 

2 

- 5 

3.0 

2.75 

3.530 

Y3 

0.5 

-2.0 

1.5 

1.25 

1.140 

E3 

600 

- 1000 

400.0 

725.00 

877.000 

Y4 

0.3 

- 0.7 

0.4 . 

0.43 

0.428 

E4 

200 

- 600 

400.0 

325.00 

509.000 

15 

7 

^ 10 

3.0 

8.25 

8.790 

Y5 

7 

^ 8 

1.0 

7.60 

7,560 

E5 

600 

- 1100 

500.0 

750.00 

1029.000 


7 

Net benefit =Rs.ll7.18 x 10 


I, Y are in MAF and E is in million Xwhr 
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TABLE 6.3.5. RESULTS OF STEEPEST ASCENT METHOD FOR SAMPLE 5 


variable 

t ' ' ' ■ 

• 

! Assumed 
; limits of 
; variable 

t 

Range 

j 

1 

; Starting 
; base valxie 

t 

t 

I 

- 1- .. - 

Final value 

i 

-1^ 2 

3 

I 4 

5 

11 

2-5 

3.0 

2.60 

4.000 

yi 

3-7 

4,0 

6.50 

5.640 

El 

50 - 150 

100.0 

90.00 

. 68.000 

12 

5-13 

8,0 

8.50 

9.600 

Y2 

10 - 16 

6.0 

14.00 

13.740 

E2 

600 - 1000 

1000.0 

685.00 

2085.000 

13 

2-5 

3.0 

3.00 

3.760 

Y3 

0;5 - 2.0 

1.5 

1.50 

1.440 

E3 

600 - 1000 

400.0 

750.00 

900.000 

Y4 

0.3 - 0.7 

0.4 

0.42 

0.418 

E4 

200 - 600 

400.0 

350.00 

529.000 

15 

7-10 

3.0 

8 50 

8.970 

Y5 

7-8 

1.0 

7.70 

7.660 

E5 

600 1100 

500.0 

775.00 

1053,000 


Net benefit 

=sRs,ll8. 

16, X 10*^ 



y are in MAP. and E is in million Ksdu:. 
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TABLE 6.3.6. RESULTS OP STEEPEST ASCENT METHOD FOR SAMPLE 6 


Variable 

1 ' 
{ Asstmed 

I limits of 

J variable 

1 1 

Range 



r— — 

Starting 
base value 

Final value 

1 

^ ■ 

2 ! 

3 

! 4 , 

5 

11 

2 

- 5 

3.0 

2.60 

4.007 

Y1 

3 

- 7 

4.0 

6.75 

5.940 

El 

50 

- 150 

100.0 

95.00 

73.000 

12 

5 

- 13 

8.0 

8.75 

9.670 

Y2 

10 

- 16 

6.0 

14.50 

14.280 

E2 

600 

- 1600 

1000,0 

700.00 

2100.000 

13 

2 

- 5 

3.0 

3.25 

4.000 

Y3 

0.5 

- 2.0 

1,5 

1.65 

1.623 

E3 

600 

- 1000 

400.0 

775.00 

923.000 

Y4 

0.3 

- 0.7 

0.4 

0.41 

0.409 

E4 

200 

- 600 

400.0 

375.00 

550.000 

15 

7 

- 10 

3.0 

8.75 

9.150 

Y5 

7 

- 8 

1.0 

7.80 

7.760 % 

E5 

600 

- 1100 

500.0 

800.00 

1078.000 


Net benefit ='Rs,119.34 x lo"^ 


Y are in MAP awi E is in million kwhr 
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TABLE 6.3.7. RESULTS OF STEEPEST ASCENT METHOD FOR SAMPLE 7 


Variable 

Assumed 
limits of 
variable 



Range 

- _J 

Starting 
base valxie 

L__ 

Final value 



2 

3 ! 

4 

^ - 

11 

2-5 

3.0 

2.75 

4.09 

Y1 

3-7 

4,0 

o 

o 

6.24 

El 

50 - 150 

100.0 

100.00 

78.30 

12 

5-13 

8.0 

9.50 

9.86 

Y2 

10 - 16 

6.0 

15.00 

14.82 

E2 

600 - 1600 

1000.0 

700.00 

2100.00 

13 

2-5 

3.0 

3.40 

4.14 

Y3 

0.5 - 2.0 

1,5 

1.65 

1.62 

E3 

600 - 1000 

400.0 

800.00 

945.00 

Y4 

0.3 - 0.7 

0.4 

0.45 

0.44 

E4 

200 - 600 

400,0 

400.00 

570.30 

15 

7-10 

3.0 

9.25 

9.52 

Y5 

7-8 

1.0 

7.80 

7.76, 

E5 

600 - 1100 

500.0 

825.00 

1101.00 


Net benefit =Rs..l20.50 x 10^ 


Y are In MAF.. and E is in million Icwhr. 








TABLE 6.3.8 


RESULTS OP STEEPEST ASCENT METHOD FOR SAMPLE 8 


variable 

1 1 

; Assxamed I 

! limits of *, 

J variable ; 

1 1 

Range 

1 

J Starting 

J base value 

1 

1 

1 

Pinal value 

1 

-J- 2 : 

3 

! 4 

-L . 

11 

2-5 

3.0 

3.00 

4.24 

Y1 

3-7 

4.0 

7.00 

6.24 

El 

50 - 150 

100.0 

100.00 

78.30 

12 

5-13 

8.0 

9.00 

9.73 

Y2 

10 - 16 

6.0 

15.00 

14.82 

E2 

600 - 1600 

1000.0 

700,00 

2100.00 

13 

0.2 - 5 

3.0 

3.40 

4.14 

Y3 

0.5 - 2.0 

1.5 

1.75 

1.74 

E3 

600 - 1000 

400.0 

800.00 

945.00 

Y4 

0.3 - 0.7 

0.4 

0.40 

0.40 

E4 

200 - 600 

400.0 

400.00 

570,00 

15 

7-10 

3.0 

9.00 

9.34 

Y5 

7 - 8 

1.0 

7.70 

7.66 

E5 

600 - 1100 

500.0 

825.00 

llOl.OO 


Net benefit 

=RS. 120.64 

X 10 ^ 



I, Y are in MAP and E is in miliion kwhr 
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TABLE 6.3.9. RESULTS OP STEEPEST ASCENT METHOD FOR SAMPLE 9 


Variable 

— I j- 

■ 1 

I AssTomed 5 

{ limits of 1 

J variable J 

< I 

t- - . 1 ^ 

Range 

i 

; starting 

; base value 

1 

1 i 

• 1 

1 -J 

Final value 

1 



2 

T" 

- 1 

3 

! .4 _.J 

_ 

11 

2 

- 

5 

3.0 

2,60 

4.00 

Yl 

3 

- 

7 

4.0 

7.00 

6.24 

El 

SO 

- 

150 

100.0 

100.00 

78.30 

^2 

5 

- 

13 

8.0 

10.00 

10.29' 

Y2 

10 

- 

16 

6.0 

15.00 

14.82 

E2 

600 

- 

1600 

1000.0 

700.00 

2100.00 

13 

2 

- 

5 

3.0 

3.40 

4.14 

Y3 

O.S' 

- 

2.0 

1.5 

1.6S 

1.62 

E3 

600 

- 

1000 

400,0 

800.00 

945.00 

Y4 

0.3 

- 

0.7 

0.4 

0.41 

0.41 

E4 

200 

- 

600 

400.0 

425.00 

590.00 

IS 

7 

mm 

10 

3.0 

9 25 

9.52 

Y5 

7 

- 

8 

1.0 

7.80 

7.76 

E5 

600 

- 

1100 

500.0 

825.00 

1101.00 


Net 

ber^f it 

= Rs,120.€l5 

X lo"^ 



I, Y are In MAP and E is in millicn kwhr 
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TABUS 6,3.10. PINAL SELECTED SET OF VARIABLES 


t 

1 

SI. No. i 

« 
t 

— - - - 1 

variable 

t 

1 

1 Selected valve 

f 

1 

11 

4.00 

2 

Yl 

6.24 

3 

El 

78.30 (9 MW ) 

4 

12 

10,30 

5 

Y2 

14.82 

6 

E2 

2100,00 (240 MIf) 

7 

13 

4.14 

8 

Y3 

1.62 

9 

E3 

945,00 (108 m) 

10 

Y4 

0,41 

11 

E4 

590.00 (67 Mtf) 

12 

15 

9.52 

13 

Y5 

7.76 

14 

E5 

1100.00 (125 m) 

I^ Y are in MAP 

and E is 

in million Icwbr. 




in storage for use later for irrigation and other beneficial 
pxirposes. A well co-ordinated plan is generally, required to 
be developed by ccaisidering demands for the various pvirposes 
and distribution of releases amcaig varioias reseirvoirs. 

An operating policy showing releases for each reser- 
voir (Tables 5.3.14 to 5.3.16) and the minimum and maxiimm 
permissible storages at the end of each month for each reser- 
voir (Tables ,5,3.17 to 5.3.19) are adopted and simulation 
carried out. A separate flood release policy is assumed for 
flood month (Table 5.3,5). Ten simulation runs are carried out 
for various cases mentioned in Section 5,6. 

Tte results of simulation studies are given under the 
following groups , 

6.4.1. Annual irrigaticn shortages 

Under the operating policy adopted, the performance of 
the system is evaluated in terms of irrigation and power shor- 
tages. Tables 6.4.1 to 6.4,10 shews annual irrigation shortages 
and their frequencies in terms of "types" of shortages. The 
"types" of shortages have already been described in Section 5.4. 
Table 6,4.1 refers to the irrigation shortages for twenty eight 
years of historical streamflows data for modified NWDT award 
(modification as explained in Section 5.4). Tables 6.4.2 to 
6.4.4 refer to irrigation shortages using fifty years each of 
synthetic flows for modified NWDT award. Similar results for 
the cases incorpearating the NWDT award (see extract in Appendix) 
are given in Tables 6.4.5 to 6,4.8. The Table 6.4.9 and 6.4.10 
refer to flood studies . 
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6.4.2. Power shortages 

Two seasons (July to September^ and, October to June) 
are considered for finding the actual power production and 
shortages during each simxilation stxidy. The types of pcwer 
shortages have been described for each season in Section 5,4, 
Table 6,4.11 refers to pcwer shortages for each seasoi for 
twenty, eight years of historical streamflow data for modified 
NWDT award. Tables 6.4.12 to 6.4.14 refer to power shortages 
for fifty years each of generated data for modified NWDT award. 
Similar results for the cases incorporating the NWDT award are 
given in Tables 6.4.15 to 6.4.18, Tables 6,4.19 and 6,4,20 
refer to flood studies. 

6.4.3. Shortage index 

Another method of indicating the shortages by a 
"shortage index" as defined by Beard (see Section 5.5) is also 
Incorporated in the simulation program. The shortage indices 
for various simulation runs are given in Tables 6,4.21 to 6.4.3C) 
Shortage index for canal (irrigation) are found to be higher at i 
each reservoir site. Suppose canal (irrigation) shortage index 
5 occxars for a 50-year simulation period, one can say that 
there may be 10 annual shortages of 50% each, or 20 annual 
shortages of about 36% each, or 50 annual shortages of about 
22% each, etc. in a 50 year period. One cannot say whether 
shortages occxirred consecutively or not. 















TABLE 6.4.13. TYPES OF POWER SHORTAGES AND TOEIR FREQUENCIES FOR EACH SEASON FOR 50 YEARS OF 

SYNTHETIC FLOWS FOR MODIFIED NWDT AWARD 
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TABLE 6.4.16. TYPES OP PC5WER SHORTAGES AND THEIR FREQUENCIES FOR EACH SEASON FOR 50 YEARS OF 

SYNTHETIC FLOWS FOR NWDT AWARD 
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6.4.4, Averages 

For each simulation run^ average monthly results are 
obtained and shewn in Tables 6.4,21 to 6,4.30. For instance 
Table 6.4 •21 shews the average monthly results for twenty eigi 
historical flows. In this tables column (2) shows the 
regenerated flew due to upstream uses (other than ccanmand arec 
under a particular reservoir)^ column (3) the regenerated flow 
caused due to water uses under immediate upstiream reservoir/ 
colxjmn (4) the inflow due to independent catchment area/ | 

colxamn (5) the release from upstream reservoir/ column (6-), 
colxjmn (7)/ and column (8) the required release for irrigation 
actual release/ and shortages jcespectively* colximn (9) the 
upstjream water jrequirements above a particular jreservoir/ 
column (10), column (11). and column (12) jrequiaired release to : 
downstream, actual release and shortages respectively. The 
shortage mentioned in column (12) has not the same meaning as 
that for irrigation shortages in column (8) because in calcu- 
lating river shortages, negative shortages are set equal to 
zero. Similar results for power are shown in Tables 6,4,31 to 
6.4.40. 


years^ 


6.4.5, Rigid jrules 

I 

Another output of the siraulaticai consists of the end ! 

l 

of month storages which can be used for the development of rigi< 

■ ; , •• i: 

jTules for each reservoir. End of menth storages are obtained 

■f 

for each hydrological state of basin like very wet, wet, average 
and dry. The mean value for each month is calculated and are | 
tabulated in Tables 6.4,41 to 6.4,50, The graphs of mean end 
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TABI£ 6.4.22. AVERAGE MONTHLY RESULTS FOR 50 YEARS SYNTHETIC PLOWS FOR MODIFIED NWDT AWARD 

(Flews in cfs) 
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calculating shortages in colimn (12)^ negative shortages are set to zero 

















TABLE 6,4.23. AVERAGE MONTHLY RESULTS FOR 50 YEARS SYNTHETIC FLOWS FOR MODIFIED NWOT AWARD 

(Plows in c£s ) 
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calculating shortages negative shortages are set to zero in colun\n(.12,) 



TABLE 6,4.24. RESULTS FOR SO YEARS SYNTHETIC PLOWS FOR MODIFIED NWDT AWARD 
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shortages negative shortages are set to zero in column ^12.) 



TABLE 6.4.25. AVERAGE MONTHLY RESULTS FOR 28 YEARS OP HISTORICAL PLOWS FOR NWDT AWARD 

(Plows in cfs ) 
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calculating shortages negative shortages are set to zero in column jl'l 2) 


TABLE 6.4.26. AVERAGE MONTHLY RESULTS FOR 50 YEARS SYNTHETIC PLOWS FOR NWDT AWARD 

(Flews in cf s ) 
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*In calculating shortages negative shortages are set to zero in column (l^ 














TABLE 6.4.27. AVERAGE MONTHLY RESULTS FOR 50 YEARS SYNTHETIC PLOWS FOR NWDT AWARD 

(Plows In cfs) 
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calculating shortages negative shortages are set to zero in coluimi 


table 6.4,29. RESULTS FOR FLOOD SIMULATION FOR 28 YEARS OF HISTORICAL FLOWS 
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calculating shortages negative shortages are set to zero in column (12| 
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calculating shortages negative shortages are set to zero 
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In calculating shortages, negative shortages are set to zero. 
Load factor = 60%, 



32. AVERAGE SEASONAL ENERGY FOR SO YEARS OF SYNDETIC PLOWS FOR MODIFIED NWDT AWARD 
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* In calculating shortages^ negative :^crt ages are set to zero 
4-Load factor = 60%, 


TABLE 6.4.33. AVERAGE SEASONAL ENERGY FOR 50 YEARS OF SYNTHETIC FLOWS FOR MODIFIED NWDT AWARD 
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In calculating shortage, negative shortages are set to zero. 
Load factor = 60%. 
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in calculating shortages, negative shortages are set to zero 
Lo^ factor = 60%. 
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In calGulating shortages, negative shortages are set to zero 
Load factor =60%, 



TABLE 6.4.36. AVERAGE SEASONAL ENERGY FOR 50 YEARS OF SYNTHETIC FLOWS FOR NWDT AWARD 
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calculating shortages, negative shortages are set to zero 
id factor = 60%, 



TABIE 6.4,37. AVERAGE SEASWAL ENERGY FOR 50 YEARS OF SYNTHETIC FLOWS FOR NWDT AWARD 
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in calculating shortages, negative shortages are set to zero. 
Load factor *= 60%. 
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In calculating shortages, negative shortages are set to zero 
Load factor =60%, 
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In calculating shortages, negative shortages are set to zero 
load factor « 60%. 



TABLE 6 . 4 , 41 . RIGID RULES FOR DIFFERENT HYDROLOGICAL STATES OP BASIN FOR 28 YEARS OF 

HISTORICAL PLOWS FOR MODIFIED NWDT AWARD 
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Dry 0,289 0.385 0.154 0.181 0.175 0.247 0,228 0.266 0.248 0.219 0.184 0,165 

Ave 0.316 0.367 0.293 0.316 0.297 0.297 0.280 0.273 0.264 0.258 0.247 0.321 

Wet 0.367 0.385 0.368 0.358 0.366 0.362 0.349 0.336 0.326 0.322 0.329 0.344 

Very wet 0.375 0.385 0,400 0.381 0.377 0.371 0.364 0.356 0.352 0.350 0.349 0.375 



TABLE 6.4,41 (Contd. 
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TABLE 6.4.48. RIGID RUI^S FOR DIFFERENT HYDROLOGICAL STATES OF BASIN FOR 50 YEARS OF 

SYNTHETIC PLOWS FOR NWDT AV7ARD 
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TABIiE 6 . 4 . 50 . RIGID RULES FOR DIFFERENT HYDROLOGICAL STATES OF BASIN FOR FLOOD SIMULATION 
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Dry 0.300 0.385 0.149 0.175 0,170 0.250 0.207 0.287 0.269 0.243 0.212 0.192 

Ave 0.317 0.368 0.294 0.329 0.310 0.293 0.298 0.288 0.284 0.274 0.274 0.325 

wet 0.365 0.385 0.369 0.377 0.378 0,373 0.360 0.347 0.344 0.344 0.339 0.347 

very wet 0,375 0.385 0.405 0.380 0,377 0.371 0.364 0.356 0.352 0.352 0.350 0.375 


TABI^E 6 . 4 . 50 < Contd 





—it CO CM m 

CD 




O CM f-t r* 





CO CTv CO 

a> ftJ 



r-t 


> 




CM CM m lA 

-H >1 









<S vO 03 o 

vJ 

>1 




<«> a» r~ O 

H tJ 



ro 

C^ 00 

O.C 



fH 

t ♦ • • 

e (0 




CM CM CO ^ 

0} 




O in i> o 

3i 
+3 D> 




ov CO O m 

H fO 



CM 

r* o\ vD \0 




r-t 

• ♦ . • * 





CM CM CO 't:!* 

• > 

fxi 




*0 (TJ 

QJ • 

2 



CO O O O 

-r-t > C2 




c^ o\ o 

44 -P 3 



r-t 

03 o^ r~ o\ 

•H 0) 

-H 


r-t 

# ♦ • • 

O 3: ’ 




CM CM CO 

a -s 

Di 

O O O m o 



CO 4J -H 

m 


O vO in 

G) 4-> 

y 

O O vO pv 


CM CO CM tH 

CO s <0 

o 

so o m CO cn 

o 

<T» r-t rH CM . 

m r-t 

4J 

« * «^ « • 

rH 

• • • ♦ 

a 

CO 

O CMO O CM 


CM CO tn 

CO M H 

0) 0) H 

M > CQ 




T-l m vD 

U 

e 



CO r- in vo 

0 

-H 



CO CO o\ 

Cn CO 0 

C 



♦ • • • 

CO CO ^ ir> 

m 00 CO cj\ 

OJ G 0 
+3 0 H 
(0 H 44 

O +3 

S 



n) M 




00 in CO 00 

y 0 0 



CO 

O r* in 

3 -H MM 




• • • f 

0 44 




^ ^ in vD 

44 -H 73 





CO y 

C 




0 CO 10 

•ri 



^ CM CO 

43 to 36 




O CJV rHI 

C r4 rO 

CD 


>- 

1-1 o CO cn 

H 0 


m o o o m 


• % • «l 




1 tn in vo vO 

C S 

U 

CM O vO ^ l> 


' 

•H 3 5 

0 

• * •• '» • 



CO 0 z 

4D 

vx> m th o r- 


O in CM SD 

tO 44 


r-t 


O O CO CO 

£1 -d 



vD 

iH r- cri CM 

^ CD 

£ 



• • ♦ * , 

<D U -H 

3 



nO vo vo 

£3 03 44 

S 




43 > -H 1 

Tj i 



CO mo o 

0) -rt 1 

X j 



>1^ 2 

m 



f-t CO m m 

iM 0 e 

S j 


in 

^ m r- r> 

•H X 





CO 43 




r- t> r-* 

1 

(0 43 

to *44 

• 




H >1 5 

0 



GO tH CO vO 

O M 




CM ^ CO t-t 

■O C 




t> CO f-t 

43 0 

u 


] 

* ♦ ■ 

Odra 

•irl 


\ 

^ m 

C G -H 

0 

r-t CM CO Sf* m 


CO M 

> 



CO cp 

y 



O Oi 00 o 

(U ' Oi 

<D 



n vO o 

O 03 £ 

CO 


fO 

CJi O O r-« 

d O' o 

S 



^ * • ' , • 

(0 0 

0^ 



CM CO CO CO 

d M 





U 03 l4 





(3 > o 





S to 44 




<D 

to 





'd 





® 4J 0) 





£ ® ra 



CM 

>1 m 4> M 

P ^ 0 





•* 




m 

1 1 




970 


rr 

of month storage versijs period (month) for a £ew selected 
cases are shewn in Figures 6.4,1 to 6,4,16. The storages at 
the end of June generally give a pictxare of carry over storagi 

6.4.6, Spills 

The statistics of monthly spills fran the terminal 
reservoir namely Sardar Sarovar are tabulated in Tables 6.4,! 
to 6.4i58 for each simulation run. It is fomd that there is 
spill during the July month because of the restriction imposed 
on the level of storage (40% storage) by the end of July. If 
this limit is increased, the spill will be shifted to other 
months say August, September, and October. There is a possi- 
bility of spilling during the month of November in very wet 
years, 

■ 1 ! 

6.4.7, Flood simulation study 

3-ho\arly maximum spill from the terminal reservoir, 
Sardar Sarovar is the output of flood simulation study. Two 
simxilation rvins are taken and the observed maximum spills are 
shewn in Tables 6,4,59, and 6.4,60. In the flood simulatloii 
studies ireported here, the lag method of flood r cubing is used 
for channel routing. In this method the ordinates of flood 
hydrographs are lagged by the time of travel between the two 
reservoirs. The time of travel from odb reservoir to another ; 
has been given in Chapter 3 (Secticai 3.4) . 

From Tables6.4.59, and 6.4.60, it is seen that there 
is more attenuation of flood in the case of NWDT than in the 
case of modified NWDT. - | 
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TABLE 6.4.51. STATISTICS OF SPILL AT SARDAR SAROVAR FOR 

28 YEARS OF HISTORICAL FLOWS FOR MCHDIFIED NWDT 
AWARD 


t 

1 

Month 1 

t 

f 

1 

t 

Mean 

cfs 

; standard 
; deviation 

1 cfs 

1 

1 

t 

; Maximum spill 
; cfs 

1 

.-.-1 - - -1 

Minimum spil 
cfs 

' ' - 

July 

42535 

6074 

157017 

0 

Aug 

14021 

2413 

- 92246 

99 

Sept 

36457 

247 

200000 

16 

Oct 

15681 

2238 

162977 

16 

Nov 

6712 

342 

41311 

16 

Dec 

2056. 

321 

12856 

'i6 /: j 

Jan 

14 

0 

16 

0 ■ ^ 

Feb 

9 

1 

16. 

0 i 

March 

5 

2 

16 

O' i 

April 

5 

2 

16 

0 

i 

May 

5 . 

2 ’ 

16 

0 

June 

888 

158 

18342 

' ■ ; ■ . P 

0 ■ i, 
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TABLE 6.4.52. STATISTICS OF SPILL AT SARDAR SAROVAR FOR FIRS 

50 YEARS OF SYNTHETIC FLOWS FOR MODIFIED NWDT 
AWARD 


Month 

Mean 
cf s 

Standard 
deviation 
cf s 

L.,.,. 

1 

1 

t 

1 Maximum spill 

J cfs 

1 

1 

1 

Minimum spil. 
cfs 

July 

26659 

1293 

147288 

16 

Aug 

10289 

1456 

128219 

99 

Sept 

23970 

3422 

213939 

16 ! 

Oct 

6366 

907 

121456. 

16 

Nov 

2025 

287 

40490 

16 

Dec 

692 

97 

9751 

16 

Jan 

13 

0 

16 

0 

Feb 

7 

1 

16 

0 i 

March 

3 

0 

16 

0 ! 

April 

1 

0 

16 

0 

May 

1 

0 

16 

0 


Jvine 


497 


69 


20390 


0 
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TABLE 6.4.53, STATISTICS OF SPILLJ- AT SARDAR SAROVAR FOR SECOND 

50 YEARS SYNTHETIC FLOWS FOR MODIFIED NWOT AWARD 


f 

1 

Month J 

t 

f 

t 

1 

t 

Mean 

cfs 


1 

1 

1 

; Maxim'um spill 
; cfs 

t 

1 

• - - 

Mininnim spiJLl 
cfs 

July 

38418 

17664 

144966 

O 

Aug 

14670 

2804 

105817 

99 

Sept 

31726 

1775 

143591 

16 

Oct 

10517 

2021 

142299 

16. 

Nov 

4026 

772 

37540 

16 

Dec 

1393 

265 

10196 

16 

Jan 

14 

0 

16. 

O 

Feb 

10 

1 

16 

O 

March 

5 

1 

16 

O 

April 

4 

1 

16c. 

O 

May 

4 

1 

la 

O 

June 

394 

76.- 

19450 

0 


[ 
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TABLE 6.4.54. STATISTICS OF SPILL AT SARDAR SAROVAR FOR THIRD 

50 YEARS SYNTHETIC PLOWS FOR MODIFIED NWDT AWARD 


f 

ff 

Month 1 

1 

• 

1 

1 

t 

Mean 

cfs 

Standard 

deviation 

cfs 

- j - . - - J 

Maximum spill 
cfs 

f 

t 

; Minimum spill 
; cfs 

f 

f 

July 

41964 

2678 

125789 

0 

Axig 

11901 

1633 

103290 

99 

Sept 

47285 

6014 

667986 

16 

Oct 

15532 

2986 

152976 

16 

Nov 

8175 

142 

68107 

16 

Dec 

2110 

996 

14669 

16 

Jan 

14 

0 

16 

0 

Feb 

10 

1 

16 

0 

March 

7 

2 

16 

0 

April 

7 

2 

16 

0 

May 

5 

2 

la 

0 ^ 1 

June 

1512 

288 

36171 

0 
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TABLE 6.4,55. STATISTICS OF SPILL. AT SARDAR SAROVAR FOR 

28 YEARS OF HISTORICAL FLOWS FOR NWDT AWARD 


1 

» 

1 

Month I 

r 

t 

1 

I 

Mean 

cfs 

» 

1 

} 

; standard 

J deviation 

1 cfs 

1 

1 

1 

f 

I Maximxjm spill 

; cfs 

1 

1 

1 

1 

t 

1 

1 Minimum spill 

I cfs 

1 

...1 , - - - 

July 

28261 

5436 

144337 

16 

Aug 

10013 

1556 

67842 

99 

Sept 

33108 

1379 

179343 

16 

Oct 

7829 

1504 

131306 

16 

Nov 

1638 

312 

14653 

16 

Dec 

16 

0 

16 

16 

Jan 

14 

0 

Id. 

6 

Feb 

10 

2 

16 

0 

March 

6 

1 

Id 

0 

April 

3 

1 

Id 

0 

May 

2 

0 

16 

0 

June 

5 

2 

Id 

0 
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TABLE 6.4.56. STATISTICS OF SPILI* AT SARDAR SAROVAR FOR 

50 YEARS OF SYNTHETIC FLOWS FOR NWDT 
AWARD 


1 

1 

t 

Month ; 

i 

t 

1 

i 

1. 

Mean 

cfs 

1 

1 

t 

J Standard 

I deviation 

1 cfs 

t 

1 

1 Maximxm spill 

! cfs 

1 

1 

.... 1 ... , , . 

Minimtan spill 
cfs 

July 

18851 

3625 

128305 

0 

Aug 

10334 

1970 

113920 

99 

Sept 

22139 

4258 

179382 

16 

Oct 

4216 

808 

75420 

16 

Nov 

560 

105 

18284 

16 

Dec 

15 

0 

16 

0 

Jan 

13 

1 

16 

0 

Feb 

9 

1 

16. 

0 

March 

7 

2 

16 

0 

April 

4 

2 

16 

0 

May 

2 

3 

16 

0 

June 

4 

2 

16 

0 
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TABL£ 6.4.57. STATISTICS OF SPILL AT SARDAR SAROVAR FCR SECOND 



50 

YEARS SYNTHETIC 

FLOWS FOR NWDT 

AWARD 

t 

Month 1 

1 

1 

t 

1 

t 

Mean . 

cf s 

t 1 

; standard I 

I deviation I 

! c£s ; 

1 f 

1 

1 

1 

Maximum spill [ 

cfs 1 

1 

1 

- - i. 

Minimum spill 
cfs 

July 

22733 

9411 

192653 

0 

Aug 

12129 

2315 

82696 

99 

Sept 

27031 

1063 

127788 

16; 

Oct 

7355 

1412 

123265 

16 

Nov 

871 

165 

15913 

16 

Dec 

16 

0 

16: 

16 

Jan 

12 

1 

16; 

0 

Feb 

9 

2 

16 

0 

March 

6l 

1 

16: 

0 

April 

2 

0 

16 

0 

May 

2 

0 

1€d 

0 

Jione 

4 

1 

16 

0 
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TABLE 6.4.58. STATISTICS OF SPILL AT SARDAR SAROVAR FOR THIRD 

50 YEARS SYNTHETIC FLOWS FOR NWDT AWARD 


1 

1 

Month 1 

1 

S 

f 

1 

. - 

Mean 

cfs 

t 

1 Standard 

deviation 

• cfs 

1 

-t . 

1 

1 

J Maximum spill 

I cfs 

1 

f 

..... 1 . - - .. 

Minimtim spill . 
cfs 

- - - - 

July 

20032 

4030 

85664 

0 

Avg 

10395 

1982 

87758 

99 

Sept 

41058 

7898 

249514 

16 

Oct 

9795 

1882 

143188 

16 

Nov 

2598 

497 

29411 

16 

Dec 

16 

0 

16 

0 

Jan 

15 

0 

1& 

0 

Feb 

7 

2 

16- 

0 

March 

7 

1 

16 

0 

April 

5 

1 

16> 

0 

May 

3 

1 

16 

0 

June 

4 

1 

16 

0 
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TABLE 6.4.59. STATISTICS OP SPILL AT SARDAR SAROVAR FOR PLOCO 

SIMULATION FOR 28 YEARS OF HISTORICAL FLOWS FOR 
NWDT AWARD 


f 

1 

Month 1 

f 

i 

t 

■1. • ■ 

Mean 

cfs 

“T 

1 

1 

• 

f 

t 

f 

« 

i 

r 

1 

I 

Standard I 

deviation | 

cfs . 

i 

1 

" ' t 

I 

1 

Maximum spill 
cfs \ 

t 

- - 1- 

Minimum spill 
cfs 

July 

35667 


6861 

124029 

16 

Aug 

14836 


2852 

59007 

99 

Sept 

40650 


7820 

136102 

16 

Oct 

17086 


925 

136102 

16 

Nov 

8281 


203 

56549 

16 

Dec 

326 


60 

2770 

16 

Jan 

270 


52 

16 

0 

Feb 

9 


2 

16 

0 

March 

9 


2 

16 

0 

April 

5 


1 

16- 

0 

May 

9 


2 

16 

0 

June 

830 


157 

15082 

0 

3-ho\irly maxlnvum spill 

= 5,42,407 

Cfs. 





253 


TABLE 

6.4.60. : 

STATISTICS OF SPILL AT 
SIMULATION FOR 28 YEARS 
MCDIFIED NWDT AWARD 

SARDAR SAROVAR 
OF HISTORICAL 

FOR FLOOD 
FLOWS FOR 

Month 

t 

1 

t 

J Mean 

r 

1 

1 

1 

1 

j j 

1 I 

1 1 

Standard J Maximum spill J Minimum; spill 

' 

! cf s 

1 

1 

deviation ; 

cfs ; 

cfs 


1 

1 

1 

1 

i 

cf s « 

1 

.. - - -...1 

f 

t 


July 

46419 


8930 

171710 

0 

Aug 

20984 


5589 

104186 

99 

Sept 

43812 


3949 

160135 

16 

Oct 

18980 


806 

176216 

16. 

Nov 

9007 


2340 

41311 

16 

Dec 

2190 


407 

12856 

16 

Jan 

14 


0 

16 

0 

Feb 

9 


1 


0 

March 

6 


2 

16. 

0 

April 

5 


2 

16 

0 

May 

5 


2 

l€v 

Q 

June 

889 


213 

18342 

0 

3-hour ly maximum spill 

= 11,00,840 Cfs 

• 







254 


6.4.8. Econonic evaluation 

For determining the net annual benefits achieve<^ from 
each simulation run, loss fxinctlons for irrigation and power 
are incorporated in the objective fxinction (equation 4.3), Ihe 
loss functions for irrigation and power are expressed by 
equations 5,11, and 5.12, The net annual benefits for each 
simulation is tabulated in Table 6,4^61. 

TABLE 6.4.61. ECONOMIC EVALUATION OF THE SYSTEM 


Simulation 

I j 

t 

I Annual benefits 

1 

J Average 

i 

1 

} Remarks 

rxm 

I in crores 

1 

t 

1 

f 

1 

t 

f 

1 

. 1 - - - - 

1 

69.02 



2 

62.41 

66.37 

Modified NWDT award 

3 

67.80 



4 

66.25 



5 

77.87 



6 

66.84 

73.93 

NWDT award 

7 

73.85 



8 

77.17 



9 

i 

80.13 

80.13 

NWDT award 

10 

71.53 

71.53 

Modified NWDT award 
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The difference between NWDT and modified NWDT is 
e3<plained in Section 5.4. FrcMii Table 6*4,61/ it is seen that 
the average annual benefits, averaged over four simulation runs 
or 178 years of simulation, in the case of NWDT award is higher 
than that in the case of modified NWDT award. It is also noted 
that flood peak reduction in the case of NWDT at Sardar Sarovar 
is more than that in the case of modified NWDT. In the former 
case a flood magnitude of 24,50,000 cusecs gets attenuated to 
5,42,407 cusecs, while in the latter case it gets attenuated 
to 11,00,840 cusecs. It is not easy to calculate benefits due 
to reduction in flood damage . Hence no attempt has been made 
here to quantify these benefits. 

6.4,9i, Comparison of optimizing model results with 
simulation model results 

LPD model gives a wide range in the magnitiKie of the 

decision variables as it is run for various configurations and 

conditions mentioned in Section 6.2. If the design is cotrple- 

tely based on these results it may lead to either over optimal 

or under optimal design. The net benefit from the system 

7 7 

varies from 107.5 x 10 to 268.22 x 10 rupees. Therefore the 

LPD model is used for preliminary screening purpose. 

Pxarther screening of the variables is done by the 

steepest ascent method. It seen that the changes in the 

magnitude of each variable changes the optimal value because 

of the further screening. The net benefit from the system is 

'7 

reduced to 120,64 x 10 rupees. 
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Results of the simulation model give the performance 
of the system as the selected system is analysed with a time 
series inflows. The performance is measured in terms of shor- 
tages in irrigation and power. High penalty was introduced 
for not meeting the target for irrigation and power demand. 

Prom the results of simulation it is seen that power shortages 
are more which means that installed capacity of turbine at 
various sites required to be reduced. Prom the simulation 

7 

results the net benefit from the system varies from 66.37 x 10 
7 

to 80.13 X 10 rx^ees. 



CHAPTER 7 


SUMMARY AND RECOMMENDATIONS FOR FUTURE WORK 


7.1. Summary 

A combined use of optimization and simulation techniques 
proves to be efficient for planning and designing of large 
water resource systemjg . This technique is applied to the 
Narmada river basin in this thesis . The system consists of 
five serially linked reservoirs located on the main stem of 
river. A mathematical model is used to represent the system. 
Four periods in a year are selected for LPD model. This model 
is run for different cases such as different configurations , 
different flow conditions and linear and nonline^ objective 
function to obtain 14 design variables. Deterministic flows 
are used in the model. This model is planned for preliminary 
screening purpose. 

If the system analysis is based only on the LPD model, 
the design of system may be overoptimal or underoptimal . Hence 
the results of the LPD model are needed to be screened fiirther. 

A range for each variable is assumed on the basis of LPD results 
for use in the steepest method of sampling, Ibe results of the 
steepest ascent method are lased in simulation study. A set of 
design variables which gives maximum benefit based on steepest 
ascent method is selected for simulation study. 

Simulation attempts to study the performance of a 
system for specified operating policies. The operating policies 
include (i) release policy, (ii) the maximiim and minimum 
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permissible storage for each period and for each reservoir. 
Simulation study requires long sequences of streamflcws. Since 
cnly 28 years of streamflow record is available at the gauge 
sites on Jamtara Mortakka and Garudeshwary synthetic flows are 
generated \ising stochastic hydrology preserving important 
statistical parameters of the historical streamflcws. HEC-4 
computer program: 'mcanthly streamflcw simulation* is used to 
generate 250 years of synthetic flows at each reservoir site^ 
using transferred streamflow data at each reservoir site. 

Two types of simulation studies are made (i) monthly 
simulation and (ii) flood simulation. Four monthly simulation 
runs are taken for modified NWDT award and for NWDT awaord/ using 
historical and synthetic flows respectively. Flood simulation 
study is carried out to find the maximum 3 -hourly spill frcxn 
the terminal reservoir namely Sardar Sarovar. Economic evalu- 
ation of the system is also carried out using certain simplified 
loss functions. Rigid rules are developed for each reservoir 
using end of month storages. Shortage index and types of 
shortages for irrigation and power are evaluated at each reser- 
voir site and for each simulation run. The loss functions are 
necessary to take into account the contingency of not being, 
able to meet irrigation and power targets. Monthly statistics 
of spills from the terminal reservoir are evaluated for each 


simulation run 
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7 . 2 . 

(i) 


(i±) 


(ill) 


(iv) 


(v) 


Reconunendations for Future Work 

In the present study, the discontinuous LP (LPD) is 
used for optimization. It is possible to use stochastic 
LP technique in this phase. 

The stochastic nature of demand should receive more 
attention. Fxirther research work is needed to develop 
the technique for the generation of equally likely 
sequences of the different types of demands. 

Weekly or ten daily simulation is probably better 
suited for the system basically operated for conser- 
vation pvirposes. For stach simulation stxadies, adequate 
techniques of generation of weekly or ten daily seqxje- 
nces of streamflows should be developed. 

There is a lot of scope for improving the loss fmc- 
tions used to account for losses incxirred in not meeting 
target irrigation and power demands. 

To reduce the spill, it is reccsrunended to carry out 
further studies by changing the demand pattern, release 
policy and maximum permissible storage in each period. 
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APPENDIX 


COST BENEFIT FUNCTIONS OP PROJECTS 


NARMADASAGAR 


(a) Annual benefit front, irrigation; See Figure Al, 

ABFl = [c.l(P/A, 5, 10) X (P/F, 5/ 10) + 0.02(P/G, 5^ 10) 

X (P/F, 5# 10) + 0,3(P/A^ 5/ 20) X (P/P/ 5/ 20) 

+ 0.03(P/G/ 5/ 20) X (PA/ 5/ 20) + 0.9(P/A/ 5/ 15) 
X (P/F/ 5/ 40) + 0.0007 (P/G/ 5/ 15) (P/P/ 5/40) 

+ 1(P/A/ 5/ 15) (PA/ 5/ SS'Q X (AA# 5/ 70)/, 


where ABPl stands for annual benefit factor for irrigation/ 
P/A for series present-worth factor/ P/G for uniforn>> 
gradient- series present- worth factor/ P/F for single- 
payraent present-worth factor/ and AA for capital-recovery 
factor. Using 5% discovint rate and the rate of develop- 
ment as shown in Figure Al/ ABFl works out to be 0.2975. 

A delta of 2.57 ft at the head regulator and an intensity 
of irrigation of 135% are assTjmed in calcxilating irrig- 
ation water reguirements . The long term irrigation 
benefit is taken as Rs. 450/acre irrigated. Therefore 
the benefit per acre ft of water will be; 


450 

2,57/ X 1.35 


Rs* 129.7./acre ft or Rs. 129,7 X 10^/MAF, 


Annual benefit from irrigation will be Rs. 


129.7 X 10^ X 
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ABFl per MAP, i.e. Rs. 12.97 x 10*^ x 0.2975 « Rs. 

3,8564 X 10*^/MAF. 

(tt) Annual cost of irrigation works: 

Irrigation cost for continuing projects is taken from the 
Journal of Irrigation and Power, Volume 3, 1978, p. 95, 
as Rs. 3360/hectare for Madhya Pradesh where Narmadasagar 
is situated, Ass\iming a uniform esqpenditure over 25 year 
development period and a time horizon of another 45 years, 
the annual coital cost recovery works out to be Rs. 1.150 
X lo’^/MAP . 

(c) Operation, maintenance and replacement (OMR) cost: 

OMR cost is assumed as Rs, 7/acre, i.e., Rs. 2*0 23/acre 
or Rs. 0.2023 x lo’^/MAF. 

(d) Reservoir cost; 

For Narmadasagar, the capital cost is 120.0 crores. The 
gross storage is 9.9 MAF. 

Capital cost/MAF =» => Rs. 12.12 x loV^AF 

Assume a construction period of 10 years, discount rate 
5% and total time horizon of 70 years. Using these, the 
annual capital cost recovery works out to be 
Rs. 0.485 x lo’^/kAF. 
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(t^) Annual cost of energy: 

At present there is a provision to provide RBPH (river 
bed power house) at Narmadasagar dam. 

Proposed capacity of RBPH: 8 x 125 = 1000 MW. 

Capital cost chargeable to power = 340 crores. 

7 

Cost/kw installed = ibb5 x ' i§00 ~ 

.*. Cost/1 O^ kw = Rs, 340 X lo"^/10^ )<w. 

Assijming a uniform expenditxire over 10 years of installa- 
tion period and a time horizon of another 60 years, the 
annual capital cost recoverry for hydropewer works out to 
be Rs. 0.00155 x loVlO^ kwhr. 

(^) Operation, maintenance and are placement (OMR) cost: 

OMR cost is assumed as Rs. 34/kw installed, i.e., 

Rs. 3,4 X loVlO^ kw, or Rs. 0.00039 x 10^/10^ kwhr. 

(i|) Energy benefit: 

The rate of development of energy for Narmadasagar is 
shown in Figure A2. Using discoaont rate of 5% and a 
development period as shown in the figure, the annual 
benefit from sale of energy at a rate of Rs. 0.50/kwhr 

n £1 

comes to Rs. 0.02667 x 10 /10 kwhr. 

Thus, for Narmadasagar project, it is seen that the 

' ■ 

annual gross benefit from irrigaticn is Rs. 3,8564 x 10 /MAF, 

7 

annual cost of irrigation Rs, 1.150 x 10 /MAF, OMR cost of 

7 

irrigation Rs. 0.2023 x 10 /MAF and annualized reservoir cost 
Rs, 0.485 X 10^ /mF, Also gross annual benefits from power is 



NAHM A U 





SGS,SCP 



RS. 0.0266,7 X loW l=«hr- generatxon 

RS. 0.001296 X loVlo" l^hr ar^l OMR coat 6or power Rs. 0.00039 
X loVlO^ I»hr. using each of these unit coats and benefits 
as a guide, cost-capacity curves, benefit-yield curves, 
cost curves, etc, are prepared and shown in Figures AS and 

A12 to A17 . 


OMKARESHWAR 

(a) Annual benefit from, irrigations See Figure A3. 

X CPA/ 5, 20) + 0.75(P/A, 5, IS) CPA/ 5, 40) 

+ 0.0067CP/G/ S- 15) (PA/ 5/ 40) + KP/A/ 5, 15) 

X (pA*/ 5/ 55)J X (aA/ 70); 

where ABF2 stands for annual benefit factor for irrxgatron, 

P/A for aeries present-worth factor, PA £°p uniform- 
gradient-series present-worth factor, P/F for single 
payment present-worth factor a.^ aA £or capital-recovery 
factor, using 5^ discount rate a»l the rate of development 
as Shown in Figure A3, ABF2 worlts out to be 0.2158. 
delta of 2.57 ft at the head regulator and an intensity of 
irrigation of 135% are assuned in calculating irrigation 

, ta The long term irrigation benefit is 

water requirements. The i.ong 

n /acre irrigated. Therefore the benefit 

taken as Rs, 437/anre 

per acre ft of water will be 

437 . „ na. 125 . 95 /acre ft or RS. 135.95 x lO^A^F. 

TTSTTTr^B 
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Annual benefit from irrigation will be Rs, 125-95 x 10 
X ABF2 per MAF, i.e., Rs. 12.95 x 10*^ x 0.2158 « 

Rs, 2.7/ia X lo'^/MAP. 

(b) Annual cost of irrigation works: 

Irrigation cost for a continuing project is taken from 
the Journal of Irrigation and Power, Volume 3, 1978, 
p. 95 as Rs, 3360/hectare for Madhya Pradesh where 
Omkareshwar is situated. Assuming a xiniform eitpenditure 
over 15 year development period and total time horizon 
of 70 years, the annual coital cost recovery works out 
to be Rs. 0.87 x IoVmAF. 

(c) Operation, maintenance ar^ replacement (OMR) cost: 

OMR oost is assvimed as Rs, 12/acre, i.e,, Rs. 3,5/acre ft, 
or Rs. 0.35 X 10*^/MAF. 

(d) Reservoir cost: 

For Omkareshwar, the capital cost is Rs, 76 crores. The 
gross storage is Rs, 1.216 MAP. 

Capital cost/hAF = “ Rs, 62.5 x lo’^/MAF. 

Assume a construction period of 10 years , discoimt rate 
5% and total time horizon of 70 years. Using these. 



the annual capital cost recovery works out to be Rs. 1*530 

X ioVmaf. 

(^) Annual cost of energy: 

At present there is a provision to provide RBPH at 
Qnkareshwar dam. 

Proposed installed capacity of RBPH: 6 x 65 = 390 MW 
Capital cost chargeable to power =138 crores 

7 

.*. CostAw installed = 396 x ' i§Q0 "" 3540/kw 

Ass\jme a loniform expenditure over 10 years of installation 
period and total time horizon of 70 years. Using these# 
the annual capital cost recovery of hydropower works out 
to be Rs. 0.00099 x 10^/10^ )<whr. 

(J") Operation, maintenance and replacement (CMR) costs 

OMR cost is assumed as Rs. 35/kw installed, i.e,, Rs. 3,5 x 
loVlO^ kw or Rs. 0.00040 x loVlO^ kwhr. 

(||^) Energy benefit : 

The rate of development of energy for Qnnkareshwar is shown 
in Figure A4, Using discount rate of 5%, and, a develop— 
Trent period as shown in the figure, the annual benefit frcnm 
sale of energy at a rate of Rs. 0,50/kwhr comes to 
Rs. 0.01900 x ICp /IQ^ kwhr. 

Thus, for Qmkareshwar project, it is seen that the 

1 '"'' 

annual gross benefit from irrigation is Rs, 2.71,8 x 10 /MAF, 
annual cost of irrigation Rs. 0.87 x IoVmaP/ OMR cost of irri- 
gation Rs. 0.35 X IoVmaf and annualized reservoir cost 
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7 

Rs. 1.53 X 10 /MAP. Also gross annual benefits from power is 
V 6 

Rs, 0.0190 X 10 /10 l<whr, annual cost of power generation 

Rs. 0.00099 X lo'^/lO^ kwhr^ OMR cost of power Rs. 0.00040 x 

7 6 

10 /lO kwhr. Using each of these unit costs and benefits as 
guide, cost-capacity curves, benefit-yield curves, OMR cost 
curves, etc., are prepared and shown in Figures A9 and A12 to 
All. 

MAHESHWAR 

(a) Reservoir cost: 

For Maheshwar, the capital cost is 42 crores. The gross 
storage is 0,4 MAF 

AO 7 

Capital cost/MAF = = Rs. 105 x 10 /MAP 

Assume a construction period of 10 years, discount rate 
5%, and total time horizon of 70 years. Using these, the 
annual capital cost recovery works out to be Rs. 2.57 x 

7 

10 /MAF. 

(b) Annual cost of energy: 

At present there is a provision to provide RBPH at 
Maheshwar dam. 

Proposed installed capacity of RBPH: 6 x 40 = 240 MW 
Capital cost of power = 125 crores 

195 •jf lo"^ 

Cost/kw installed - 245 x ' lUo S ~ Rs. 5200/kw 

Assume a uniform expenditure over 10 years of installat lean 
period and total time horizon of 70 years , Using these / 
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the annual capital cost recovery of hydropower worlcs out 
to be Rs. 0.001455 x loVlO^ l<whr. 

(c ) Operation, maintenance and replacement (OMR) cost: 

OMR cost is assxomed Rs. 52/kw installed, i.e., Rs. 5.2 x 
loVlO^ Tcwhr, or, Rs. 0.000593& x loVlO^ l<whr. 

(d) Benefit from energy: 

The rate of development of energy for Maheshwar is shown 
in Figxjire A5 . Using discount rate of 5% and a develop- 
ment period as shown in the figure, the annual benefit 
from sale of energy at a rate of Rs. 0.50/kwhr comes to 
Rs. 0.002221 X lo’^AO^ ]<whr. 

Thus, for Maheshwar project, it is seen that the annual 

7 6 

gross benefit from energy is Rs, 0,002221 x 10 /lO kwhr, 

7 6 i 

annual cost of power generation is Rs, 0.001455 x 10 /10 Iwhr# ; 

7 6 

OMR cost for power Rs. 0.00059 x 10 /lO kwhr and the annual- 

7 ' 

ized reservoir cost is Rs. 2.57 x 10 /MAF. Using each of these 
unit costs and benefits as guide, cost-capacity cxorves, benefit- 
yield curves, OMR cost curves, etc., are prepared and shewn in j 
Figures AlO and A15 to A17. ! 

SARDAR SAROVAR 

.'■t 

. I 

(a) Annual benefit from irrigation: See Pigxore A6. i 

ABF3 = ’[o,13(P/A, 5, 10) x (PA# 5, 10) + 0.017(P/G, 5, 10) 

X (P/F, 5, 10) + 0.3(P/A, 5, 20) X (P/F, 5, 20) 

+ 0.035(P/G, 5, 20)(PA/ 5/ 20) + 1(P/A, 5, 30) | 

X (P/F, 5, 4oR X (A/P, 5, 70)j 



MAHE5HWAR 
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where ABP3 s-tands for annual benefit factor for irrigation/ 
P/A for series present-worth factor# P/G for uniform- 
gradient-series present -worth factor# P/F for single- 
payment present-worth factor# and A/P for capital-recoveiy 
factor. Using 5 % discoiint rate and rate of development as 
shown in Figure A6.> ABF3 works out to be 0,32. A delta of 
2.57 ft at the head regxilator and an intensity of irriga- 
tion of 135% are assimed in calculating irrigation water 
requirements. The long term irrigation benefit is taken 
as Rs. 625/acre irrigated. Therefore the benefit per 
acre ft will be: 


625 

2.57 X 1.35 


Rs, 180/acre ft# or# 
Rs, 180 X 10^/MAF. 


Annual benefit from irrigation will be Rs, 180 x 10^ x 
ABF3 per MAF# i.e,# Rs, 18 x 10 ^ x 0.32 = Rs. 5,763 x 
10*^/hAF. 

(b.) Annual cost of irrigation works; 

Irrigation cost for continuing projects is taken f rcxn the 
Journal of Irrigation and Power# Volume 3# 1978# p. 95 as 
Rs, 2790/hectare for Gujarat where Sardar Sarovar 
situated. Assuming a uniform expenditxire over 20 years 
development period and a time horizon of another 50 years# 



the annual capital cost recovery works out to be 
7 

RS. 1.052 X 10 /MAP. 

(<J) Operation, maintenance and replacement (OMR) cost: 

6 

OMR cost is assumed Rs. 10/acre, i.e., Rs. 2*89 x 10 /MAP, 

7 

or, Rs. 0.2890 x 10 /MAP. 

(4) Reservoir cost: 

For Sardar Sarovar, the capital cost is 200 crores. The 
gross storage is 7.75 MAP. 

. . Capital cost/MAF = = Rs, 25.80 x 10 '/MAP. 

Assume a construction period of 10 years, discount rate 5%, 
and total time horizon of 70 years. Using these, the 
annual capital cost recovery works out to be Rs. 1.03 x 

7 

10 /MAP. 

(g) Annual cost of energy; 

At present there are provisions for RBPH and CHPH at 
Sardar Sarovar, In final stage the energy frcxti RBWi will 
become nil. 

Proposed installed capacity of RBPH and CHPH: 1200 MW 
(RBPH 6 X 150 = 750 MW; CHPH 6. x 75 = 450 MW) 

Capital cost chargeable to power = 276 crores 

7 

Cost/kw installed = i§o ' 5 ' "x ' I § o 5 ~ Rs, 2300/0cw. 

Assuming a unrif orm expenditure over 10 years of install- 
ation period and a time horizcai of another 60 years, tte 
annual capital cost recovery for hydropower works o\Jt to 
be Rs. 0.00105 X loVlO^ kwhr. 



(i^) Operation, maintenance and replacement {OMR) cost: 

OMR cost is ass\amed as Rs. 23/kw installed, i.e. 2.3 x 
loVlO^ Kw, or, Rs. 0.00026 x 10^/10^ Icwhr- 

(^) Energy benefit: 

The rate of developT>ent for energy for Sardar Sarovar is 
shown in Figure A7. Using discount rate of 5% and a 
development period as shown in the figure, the annToal 
benefit from sale of energy at the rate of 0.35/kwhr ccxnes 
to Rs. 0,01867 x‘ loVlO^ Iwhr. 

Thus, for Sardar Sarovar project, it is seen that the 

7 

annual gross benefit from irrigation is Rs, 5.763 x 10 /MAF, 

7 

annual cost of irrigation Rs. 1,053 x 10 /MAF, OMR cost 

7 

Rs. 0,2890 X 10 /MAP, and annualized reservoir cost Rs, 1.03 x 

7 

10 /MAF. Also gross annual benefit from power is Rs, 0.01867 x 
7 6 

10 /lO kwhr, annual cost of power generation Rs. 0.00105 x 

7 6 7 6 

10 /lO kwhr and CWR cost for power Rs. 0.00026 x 10 /lO kwhr. 

Using each of these unit costs and benefits as guide, cost- 

capacity curves, benefit -yield curves, OMR cost curves etc., 

are prepared and shown in Figures All to A17. 

EXTRACT FROM THE NWDT AWARD 

(1) In the event of the available utilisable waters for allo- 
cation in any water year from 1st of July to 3Qth of June 
of the next calender year falling short of 28 MAF, the 
shortage should be shared between variotJS States in tiie 
ratio of 73 for Madhya Pradesh, 36 for Gujarat, 1 for 
Maharashtra, and 2 for Rajasthan, 
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(2) The utilisable flow of Narmada in excess of 28 MAP in any 
water year, i,e., from 1st of July to 30th of June of 
next calender year is apportioned in the following ratios 
of allocation^ i.e., 73 for Madhya Pradesh, 36 for Gujarat, 
1 for Maharashtra, and 2 for Rajasthan. 

(3) The available utilisable water in a water year will include 
the waters carried over from the previous water year as 
assessed on the 1st of July on the basis of stored waters 
available on that date. 

The available utilisable waters on any date will be 
inclxisive of return flows and exclrosive of losses due to 
evaporation from various reservoirs. 


( 4 ) 
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^2^ 0.00 }55E2( Narmada sagar ) 
^2^^ 0-00099E3{Omkareshwar} 

C2 4= 0-001455E4( Maheshwar) 





ruf * = 
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• ENB2,2= 0.0288E2- 0.0000033E2 {Narmada sagar | 
O ENB 2 4 = 0-0 18 E 4 - 0.0000021 E 4 (Sardar Sarovar ) 


Annual energy in 10 kwhr 

FIG.A18 NET ENERGY BENEFIT FUNCTION FOR NARMODA 
SAGAR AND SARDAR SAROVAR 
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